A GE ORY SUD Ss) KA FO R EN ING BIN 


VOLUME 5, NUMBER 2 Te l | uU S MAY 1953 


A QUARTERLY JOURNAL OF GEOPHYSICS 


Comparison of Microseisms in Greenland, Iceland, 
and Scandinavia 


By MARKUS BATH, Meteorological Institute, University of Uppsala 


(Manuscript received February 16, 1953) 


Abstract 


Microseismic amplitudes and periods recorded at Scoresby-Sund, Reykjavik, Bergen, and 
Uppsala in seven different cases from the years 1949—1950 are studied. At all stations the polar 
air is of essential importance for the generation of microseisms, whereas there is in general no 
close connection with the cyclone centres themselves. A coast effect is of importance for Scan- 
dinavia, whereas the source for microseisms recorded at Scoresby-Sund is located over the 
open ocean within the polar air. Standing ocean waves may be of importance at the Norwegian 
coast but in many cases not on the open ocean. The microseismic waves propagate much farther 
over the continent than along the ocean bottom. The microseismic periods in Scandinavia vary 
generally in parallelism with the amplitudes, in Iceland and Greenland generally not. Period 
minima and rapid amplitude increases are observed in Scandinavia when cold fronts cross the 
Norwegian coast. There is no sign of microseismic barriers in the Atlantic outside Scandinavia. 
The microseisms at Scoresby-Sund have a regular, group character; at the other stations they 


are generally continuous. 


Introduction 


In a previous paper (BATH, 1949, p. 145) 
I expressed the idea that the main source o 
microseismic energy may be different in 
different parts of the world. This result was 
based on the fact that in Scandinavia some 
kind of coast effect was the dominant source, 
whereas in many other parts of the world a 
cyclone effect seems to be of greatest impor- 
tance. Later I discussed these problems more 
thoroughly (BATH, 1951c) and inferred a 
difference with regard to microseisms be- 
tweén east and west coasts of the continents. 
The reason would be differences in the prop- 
erties of air masses blowing on-shore in 
cyclones—warm air blowing on-shore on 
an east coast, but cold air on a west coast. 
With a view to investigate these and related 
problems more closely as well as to get some 
basic data for possible judgment between 
various theories, a comparison of microseisms 
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in a number of cases has been made between 
Scoresby-Sund (Greenland), Reykjavik (Ice- 
land), Bergen, and Uppsala (Scandinavia). 
The microseisms on Greenland have been 
studied by LEHMANN (1949, 1951); the mi- 
croseisms at Reykjavik have not earlier been 
studied in detail, but some data are given by 
GUTENBERG (1932); the microseisms in Scandi- 
navia have earlier been studied by the present 
author. Copenhagen and Helsinki were not 
included in this study, as it is already known 
that the microseisms at these stations are very 
closely related to those at Bergen and Uppsala 


(see my papers 1952). 


Methods and materials used 


The microseisms studied are those in the 
usual period range of 4—8 sec. 
Data about the stations and instruments used: 
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Station Long. 


Height Ground Instrument 


G0 EEE a EL NE u Kae AL EU RE Fe ee ea na 


Scoresby-Sund... 70° 29’N 21° 57’W 69 m Gneiss : Galitzin NZ 
Reykjavik....... 64° 08’ N 21° 54° W 44 m Doleritic rock | Mainka E, N 
Fos BR f Bur. 60° 24’N ie) atte 28) 20 m Gneiss Wiechert E, N 
Uppsala ned etekche 59° 51’N 1723808 I4 m Granite Wiechert E, N 
The mean values of the seismograph constants 
for the seven cases studied from the years 
1949 and 1950 are as follows: 
pee OSPR) eS eee 
Reykjavik Bergen Uppsala 
Constant (Mainka) (Wiechert) (Wiechert) 
E N is N E N 
Wee oe Oe nl ee PP EN PO er EN 
Free period of pendulum............. 5.8 7.6 8.4 9.2 9.3 9.1 sec 
Static magnification..........+s..+-:- 92 58 IgI 180 188 193 
Damping ratio....................... 4-5 6.5 2.3 202 3.9 3.8 
Déviation due to friction rer 0.3 0.3 1.0 222 0.9 0.7 mm 


Scoresby-Sund 


Constant (Galitzin) 
13 N Z 

Seismometer period.....| 11.8 11.7 9.4 sec 
Galvanometer period...| 12.0 11.9 10.0 sec 
Seismometer damping...| appr. aperiodic 
Galvanometer damping.. aperiodic 
Transference factor..... 1070, LOOM OO m= secgl 
Reduced pendulumlength| 12.0 12.0 14.9 cm 
Distance galvanometer 

TORTECOLÄ Sere nee sys 100 I00 I00 cm 


In computing the ground amplitudes naturally 
the values of the constants for each separate 
case were used. The Z-records of Bergen were 
often faulty (obviously too large friction) and 
were therefore not measured. 

Of the seven cases studied each comprises 
in general four days. Representative values of 
the maximum amplitudes and the corre- 
sponding periods were measured on all records 
mentioned for every full hour + about 15 
minutes. 

There are various methods of measuring 
microseisms. The amplitude given should be 
a measure of the microseismic activity at the 
given time. It is not sufficient to give only the 
absolute maximum amplitude. In one case 
there may be a practically continuous train 
of waves with a certain amplitude, and in 
another case there may be only one or two 
waves within a 30-minute interval with the 
same amplitude as in the former case. If the 
absolute maxima were given, clearly the same 


values would be obtained in these two cases, 
but obviously the microseismic activity is 
larger in the former case than in the latter. The 
author has here used the same method as 
earlier, ie. a representative maximum is 
obtained by eye inspection of a 30-minute 
interval (method I). Another way is to measure 
for instance the five largest amplitudes within 


30 minutes and compute their mean value - 


(method II). Theoretically this would be 
better than method I, but practically there is 
no difference in reliability between these 
two methods. By increasing the number of 
measurements used in calculating a mean 
value, we could get a more representative 
measure of the microseismic activity. The 
author has therefore tried a third method 
(method III), measuring the maximum ampli- 
tudes for every minute within + 20 minutes 
of a full hour. Successive mean values, % 
(Au + A_1), 14 (Au + A9 + Ay + A), eétc’are 
formed (A; is the maximum amplitude for 
the i-th minute after the full hour, and A_; for 
the i-th minute before). The successive mean 
values are plotted against i. The mean ampli- 
tudes converge rapidly towards a final value. 
In general + 7 minutes are sufficient to get a 
mean value agreeing within + 0.01 mm 
(paper amplitude) with the + 20-minute 
mean value. For weak, regular microseisms 
only +4 minutes are sufficient, for large, 
irregular microseisms + 10 minutes may be 


required for the same result. However, even | 


the measurement of + 7 minutes (14 measure- 
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ments) is quite lengthy. But it is possible to 
get mean values that agree within a few hun- 
dredths of a mm with the 20-minute mean 
value by shorter methods, e.g. by computing 
U (Arı + Ava + Ayr + Aa + Ag + A) 
(method IV). The following table contains a 
comparison of the results of the different 
methods in three cases, measured on the Upp- 
sala Wiechert records. As the periods were the 
same in each case, the paper amplitudes may 
be immediately compared. Columns (a) 
contain the results by the various methods, 
(b) give values obtained by multiplying (a) 
with certain factors so as to bring the amplitude 
for case A equal to 1.00. 


Method 


Case I II 


mm mmimm mmimm mmimm mm 
A [1.00 1.00/0.92 1.00]0.59 1I.00]/0.64 1.00 
B [0.40 0.4010.37 0.40|0.28 0.47|0.27 0.42 
C |0.75 0.75|0.67 0.73|0.43 0.73|0.42 0.66 


The values for methods II and especially II 
and IV (columns a) are naturally lower than 
for I. But the ratios between the amplitudes 
(columns b) are in perfect agreement in I and 
IL. In III case B has got a somewhat higher value 
in relation to the other cases. This is due to 
the fact that the microseisms are continuous in 
this case with only small variations, and then 
the activity is actually somewhat higher than 
shown by methods I or II. On the whole, it 
may be said that it is just a matter of conven- 
tion which method should be used. The 
method I used here is for our purpose as good 
as any other method and in addition it has the 
advantage of being a quick method. 


Discussion of the separate cases 


Two characteristic types of microseisms are 
shown in Fig. 1. Fig. 1 a shows regular, group 
microseisms, typical for Scoresby-Sund, Fig. 
1 b shows less regular, continuous microseisms, 
typical for Scandinavia and Iceland. 

The ground amplitudes and the corre- 
sponding periods are given in Figs. s—rr. In 
order to avoid confusion it was necessary to 
displace the zero points of the scales for the 
different stations. The scale used is always the 
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a. Regular, group microseisms (Scoresby-Sund, 
May 13—14, 1949). 


b. Less regular, continuous microseisms (Uppsala, 
October 3, 1949). 
Characteristic types of microseisms. 


ies 1. 


same. The scales as well as the amplitude 
curves have been marked by S = Scoresby- 
Sund, R = Reykjavik, B = Bergen, and U = 
Uppsala. 

There were some obvious inconsistencies or 
variations of the constants for the E component 
at Reykjavik. In case 1 An & Ar at Reykjavik, 
in cases 2—4, An < Ag, and in cases 5—7, 
AN > Ag. There is no explanation for this 
behaviour except instrumental. The computed 
Ag of Reykjavik have been multiplied by 
factors, in each case determined so as to give 
the same mean value of Ag as of An for each 


I12 


interval investigated. The Ag for Reykjavik 
given in Figs. s—11 have been obtained in 
this way. The factors are 0.60 in case 2; 0.53 
in 3; 0.39 in 4; 1.49 in 5; 1.18 in 6; and 1.92 
in case 7. Ay and Ag for Reykjavik are nat- 
urally not comparable if their magnitudes are 
concerned, but their variations are comparable. 

The periods given in Figs. $—11 are the 
arithmetic means of the periods on the two 
horizontal components. Exceptions occur 
where only one of the horizontal components 
has been in operation, when the period of this 
component only is given. 

The interruptions of the curves in Figs. 
s—II, occurring in a few cases, are due to 
instrumental failures, unless otherwise men- 
tioned under the separate cases. 

The observations of ocean surface waves 
(mean maximum height H, period T, and 
direction of propagation) for every third 
hour made at the weather ship “Polar Front” 
(61° N, 2° E) are given in Fig. 4 for all our 
cases. These observations have been published 
by the Norwegian Meteorological Office, 
Oslo. There is a very clear parallelism be- 
tween the periods and the heights of the ocean 
waves. 

The ratios between the microseismic periods 
(T,) and the simultaneous ocean wave periods 
(T) have been computed in all cases. For the 
microseisms the mean horizontal period is 
used. The following table gives the per- 
centage frequency distribution of these ratios, 
comprising all cases studied. 


Sta- IT 
tion mu | 
0.4/0.5/0.6/0.7{0.8/0.0/1.0/1.1/1.2/r.3/1.4/1.5]1.6 
US 52 3:08,29 RG Sn) rl 
Ba 744276170 |K78 [ro eye | Ser I 
R | 5124|28| 18| 9 8| 4 2 2 
S ı| 6] 30| ro] 27] 7 | 11| 3 | 1 | 4 


Total Arles Pal x) | | 1 1 | I 


The number of observations used in com- 
puting this table were 224 for U, 221 for B, 
223 for R, and 213 for S. The frequency 
maximum lies around 0.6 and is significantly 
higher than 0.5. This is not only statistically 
valid but is repeated in every case for every 
station. The higher ratios (> 1.0) occur for 
shorter ocean wave periods and are almost 
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exclusively due to cases 2 and 5. On the 
standing ocean wave theory for microseisms 
(Loncuet-Hiceins, 1950) the period ratio 
should be 0.5. It could be said that the ocean 
wave phenomena are too complex to be 
represented by only one period value. How- 
ever, the wave observations refer to the mean 
maximum waves, and these should show the 
best relation to the microseisms. It can also 
be argued that the ocean wave observations do 
not refer to the point of origin of the micro- 
seisms. However, this does not seem to give a 
complete and satisfactory explanation. They 
could be expected to be valid at least for 
microseismic storms in Scandinavia originating 
on the south-west coast of Norway. But also 
in these cases significantly higher ratios than 
0.5 are generally obtained. 

For reasons of space the following discus- 
sions of the separate cases will be made as 
brief as possible and only the essential points 
will be emphasized. For a more complete 
survey of the weather development the reader 
is referred to the official weather maps. In 
this study I used weather maps for every 
third hour. One typical weather map for each 
case is given in Figs. 2 and 3. All times given 
are GMT. 


Case 1: January 5—8, 1949 (Fig. 5) 
Weather AE ee) 


The weather development consists essentially 
of the motion of an intense cyclone from im- 
mediately south of Iceland (06! s/1) towards 
NE. At 06 8/1 its centre is over northern 
Scandinavia. 


Amplitudes 


R has a very large amplitude maximum. 
Both the increase and the decrease are very 
rapid, and much more rapid than at S. This 
indicates that only a very special situation is 
favourable for large microseisms at R, and 
when the cyclone moves quickly as in this 
case, the microseismic storm is of short dura- 
tion. This situation with a rapidly moving low 
(mean velocity so km/hour 06 5/1—06" 6/1 
and 30 km/hour 06" 6/1—06! 7/1) could be 
expected to be favourable for the genera- 
tion of standing surface waves on the ocean. 

R’s amplitudes begin to increase rapidly 
at 024 $/1 and reach a maximum around 
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c. May 14, 1949, 06h (case 3). 


Fig. 2. Weather maps for cases I—3. 


20b—22h 5/1. It seems to be a double maxi- 
mum (16h 5/1 and o2h—03" 6/1) on both 
components. The amplitudes at R begin to 
increase rapidly when the low is south of 
Iceland, but the maxima are not reached until 
the low is NE of Iceland and at a greater 
distance from R. The intensity of the low has 
increased slightly in the meantime. A possible 
reason for the rapid increase at o2h 5/1 are 
standing sea waves on the east side of the 
cyclone between 4/1 and 5/1; a rapid change 
of wind direction has taken place there. The 
increase cannot in this case be ascribed to 
polar air, but the double maximum may be 
due to the approach of polar air (a cold front 
of smaller extent passes over Iceland at 18" 
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b. April 20, 1949, 06h (case 2). 


5/1; after that time it is situated over the 
ocean immediately to the south of Iceland). 

The amplitudes at S begin to increase at 
o24 $/1, ie. simultaneously with R, but 
considerably slower. The maximum is reached 
about 06! 6/1, ie. later than at R. The mi- 
croseisms at S were too large to be measured 
around the maximum. We have in case I an 
intense cyclone of relatively limited horizontal 
extent. There is no coast effect at S; the wind 
is almost parallel to the east coast of Green- 
land. But at the maximum there is an intense 
inflow of polar air over the ocean between S 
and the cyclone centre. No definite conclusion 
about the importance of the polar air at S is 
possible only from this case. Considering the 
maximum at 06! 6/1 we have to take the fol- 
lowing facts into account: 

1. The distance from S to the cyclone 
centre is unchanged from $/1 to 6/1. 

2. The cyclone intensity has increased 
somewhat (960 mb at centre at 06" 6/1 against 
970 mb at o6h 5/1). This could possibly 
explain part of the increase, but this change 
alone could hardly have such a large effect. 

3. The essential difference between 6/1 
and 5/1 is that on 6/1 the polar air has got 
large influence at S; that is not at all the case 
OM es) i: 

The origin of the microseisms at S is likely 
to be found in the polar air over the open 
ocean, just to the east of S. 

Regarding the amplitudes at U and B we 
observe a rapid increase at B at 16" 5/1; it 
occurs at U at the same time (UE 17° 5/1), 
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b. March 17, 1950, 12h (case 5). 


‘cut 


A 


c. March 26, 1950, 12h (case 6). d. October 8, 1950, ooh (case 7). | 
Fig. 3. Weather maps for cases 4—7. 


though less marked. Our earlier experience more pronounced for B. It is to be observed . 
from the amplitudes at B and U leads us to that the maximum amplitudes at U (and B) 
suspect some effect localized to the coast are not reached until now, when the cyclone 
around B (see my paper 1951 b). The in- intensity is decreasing but when the coast 
creases coincide with the passages of fronts effect is fully developed: 


over the coast at B. The amplitudes are still CoS - 
unimportant, especially at U. There is a Distance | pressure | 
further increase at U at 081—00h 6/1 (a cold Time LU LB |of cyclone | Coast effect | | 
front passes Krakenes on the Norwegian km km (L) | 
coast at 065 6/1), and maximum is reached at 

h 
U about 12 7/1, when the coast effect has CHE le Cool ce ER 
been displaced to the central part of the o6h 6/1 |z,600 1,300] 960 __ increasing 
Norwegian coast. An > Ag for U and still 06h 7/1 |1,400 1,400 980  |considerable 
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Fig. 4. Observations of mean maximum height (H), periods (T), and direction of propagation (arrows) of 
ocean surface waves at 61°N, 2°E. 
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pe L RE ı 1 A à en ı f h 


207 249 040108 CRE ME RC NO NON CO PE ATI SEE DE Al Cdi 
sh of vA eh 1949 
Fig. 5. Case 1: January $—8, 1949. 
U = Uppsala, B = Bergen, R = Reykjavik, S = Scoresby-Sund. A = amplitudes, T = periods of microseisms. 
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The coast effect is unimportant up to 03" 6/1, 
increasing from 03" 6/1 to 03'—06 7/1, and 
maximum is reached about 06 7/1. There 
seems to be no doubt that the microseisms at 
U, B depend upon a coast effect (within the 
polar air). 

Towards the end of case 1 there is another 
storm at U, B, apparently closer to the sta- 
tions: 


Station and 


Increase starts at} Maximum at 


component 
SPE eyo Sie a fe 20h 7/1 o7h 8/1 
INSEE 23h 7/1 06h 8/1 
Eee asus 24h 7/1 ogh 8/1 


UN is completely unaffected. There is no 
trace of this storm at R, S. The reason seems 
to be two smaller cyclones over the North 
Sea—Skager Rack, which may possibly give 
rise to standing ocean waves. There is no 
effect on the B coast, but rather on the Swedish 
west coast. 


The ocean wave heights (Fig. 4) have a 
maximum from about 12? 6/1 to 06" 7/1, ice. 
later than the amplitude maxima at S and 
especially at R, but earlier than at U and B. 
The general run is similar to the microseisms 
Ae Se 


Periods 


The periods of the microseisms are re- 
markably constant at S (6.0—6.4 sec), whereas 
the periods at U, B vary considerably. Tu, Tp 
increase from about 4.8 sec to 7 sec from the 
beginning up to 198 7/1. What is the reason 
for this period increase? Various possibilities 
will be considered: 


1. The source moves away: 

a) The cyclone centre does not move away 
(see above). 

b) The effective coast is displaced from the 
B coast to the central part of the Norwegian 
coast. This may contribute to the period in- 
crease, but can hardly explain the increase 
after 06 7/1. 

c) The polar air is closest at 06" 7/1. 

2. The increase of intensity of the coast 
effect (not of the cyclone) seems to be of 
essential importance for the period increase. 


3. Displacement of the source in accordance 
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with Press and Ewine (1948 and later personal 
communication): 

a) The cyclone centre moves towards 
shallower water during 6/1—8/r. 

b) Effective coast: no explanation possible 
in this way. 

c) Polar air: no explanation possible. 

The rapid decrease of Ty, Tg from 19h 7/1 
(7 sec to 4 %—s sec up to 8/1) is remarkable. 
It is hardly due to a decrease of intensity of 
the source, as the amplitudes at B increase. 
The reason is instead a new source close by. 
It is significant that Ts decreases more rapidly 
than Ty. After the period minimum at 02"—o4» 
8/1 for Tg and at 06° 8/1 for Ty they increase 
again somewhat. 

We observe that there is no maximum of 
TR simultaneous with the maximum ampli- 
tudes at R. At R, S the periods seem to be 
more independent of the intensity of the 
source and vary mostly with the distance to 
the source. In every case Tr < Ts almost 
without exception. 


Types of microseisms 


In case ı the microseisms at U are typically 
continuous, slightly irregular, during 6/1—7/1; 
there is some indication of groups on E in the 
morning of 8/1. At B the microseisms are 
continuous up to about 02403) 8/1, after 
which there are clear group microseisms; on 
BZ the groups are clear already from about 
19? 7/1. There are no pronounced group 
microseisms at R, whereas S has in this case 
as almost always beautifully developed regular 
group microseisms. 


Case 2: April 17—23, 1949 (Fig. 6) 

Weather’ (Pie. 20h) 

The weather development is similar to the 
preceding case and is dominated by a cyclone 
which first (up to 19/4) moves rapidly from 
south towards north of Iceland and then stays 
practically in the same position at about the 
same distance from Iceland, Greenland, and 
Norway. 


Amplitudes 

At R we distinguish a series of microseismic 
storms: 

1. The amplitudes increase in the beginning 
to an unimportant maximum at OPARTS | Amie: 
when the cyclone has its centre on Iceland. 


/ 
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Fig. 6. Case 2: April 17—23, 1940. 


This corresponds to case 1 with a rapid cyclone 
approaching Iceland from the south. Standing 
sea waves may be expected due to rapid 
changes of wind direction. 

2. Most striking is the very rapid increase 
starting at 14" 22/4 and reaching a very large 
maximum around o8® 23/4. This storm ceased 
about 23b— 24" 23/4. It is very important that 
this storm occurs not until the cyclone has 
moved away from R and when its intensity 
is decreasing, but, on the other hand, just 
when the stream of polar air straight from the 
north towards Iceland is most fully developed. 
In other words, the cyclone has reached its 


critical position. However, the possible im- 
portance of standing sea waves to be expected 
to the SE of Iceland cannot be excluded. Still 
more so, as at S there is only a slight increase 
of the amplitudes at the time of the large 
storm at R. 

3. There are two smaller maxima at R at 
about 24" 19/4 and 10° 21/4, probably due to 
the approach of polar air. 

The amplitude variations at S are again 
slower than at R. A more rapid increase 
starts at 175 18/4 and maximum is reached 


approximately at 20% 20/4. If L denotes the | 


cyclone centre, we have: 
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Time ee SL er 

18h 18/4 700 970 
o6h 19/4 800 970 
18h 19/4 750 970 
06h 20/4 700 970 
18h 20/4 700 975 


The polar air off the S coast is best developed 

on 20/4. This is obviously the reason for the 
microseisms at S. It does not seem possible to 
Mascribe these microseisms to standing sea 
M waves. 
The amplitudes at U and B are on the whole 
Arather unimportant in this case. This is also 
to be expected from earlier experience for 
cyclones with centres in the Atlantic. They 
“never give rise to large microseisms in Scan- 
| dinavia. There are two well marked amplitude 
increases (I and II): 


| Increase starts at 


The following table gives the distances to 
the cyclone centre (L) and the pressure of the 
isobar closest to the centre, given on the maps. 


Distance L 
Time BL UL rere 
km km 

o6h 18/4 1,400 2,000 980 
06h 19/4 1,500 1,700 970 
06h 20/4 1,600 1,800 970 
06h 21/4 1,400 1,600 985 
06h 22/4 1,200 1,400 980 


It is obvious that the second storm (II above) 
| has no relation to the distance to the cyclone 
| centre, nor to the central pressure, but just 
occurs when the coast effect around B has 
| reached its maximum. It is to be observed 
| that during the second storm the cyclone 
centre is all the time over the Atlantic far 
outside the continental shelf. 
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119 


The ocean wave heights (Fig. 4) have only 
one maximum of relatively long duration. 
The variation is similar to the microseismic 
amplitude variation only at S, but not at the 
other three stations. 


Periods 


With regard to periods we observe that 
Tu varies on the whole in good parallelism 
with the amplitudes at U. That is not the 
case at B, probably due to relatively greater 
importance of near-by parts of the coast around 
B. There are minima of Tg at 03" 20/4 and 
12b— 13% 21/4. The last minimum is especially 
clearly marked and is also observable as a 
decrease of Ty. This period minimum is 
simultaneous with the onset of the storm on 
21/4, i.e. the same phenomenon as was found 
in case ı and in several other cases. The period 
minimum simultaneous with the onset of a 


Remarks 


ee Ui. .| 2th 18/4 and 054 19/4; 13h] poorly Cold front passes the B coast about 15h 18/4 
19/4 developed 
LS about the same time as U,| poorly 
but more gradual developed 
LU. .|ııh 21/4 (only E) 19h 21/4 This maximum is the largest for the whole in- 
Be: 11h 21/4 (especially N) ıgh 21/4 terval studied for UE, BN, and BE. The coast 


activity is limited to the B coast. Cold fronts 
are coming in over the B coast after 18h 21/4 


storm is due to a new source of microseisms 
situated close to the station and which begins 
to dominate the situation. This source is 
here located at the B coast. After the minimum 
Tp increases again and reaches a maximum at 
about the same time as the maximum ampli- 
tude. This shows that at the coast effect the 
period depends also on the amplitudes. 

Concerning Tr it is remarkable that it 
does not change during the large amplitude 
maximum towards the end of the interval. 
This fact is not in favour for an explanation 
by means of standing sea waves. 

Ts is as usual remarkably constant = 6 sec. 
But in the beginning of the interval Ts 
decreases from 6 sec to a minimum (5.2 sec) 
at 22h 18/4—or® 19/4 and then increases 
again, whereas TR at the same time increases 
from about 4 sec and reaches 5.2 sec at the 
time of minimum T's. It is essential to note 
that the Ts minimum coincides in time with 
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the more rapid amplitude increase, ie. a 
completely analogous behaviour to what we 
just found for B. The explanation is also for 
S very probably that a near-by source (polar 
air) begins to dominate the situation. During 
19/4 Ts is again increasing, while the source 
does not move away but the amplitudes 
increase. 


Types of microseisms 


The microseisms are continuous on the U 
records. There is some indication of group 
formation on the E component at U in the 
storm on 21/4. Also at B there are continuous 
microseisms up to the storm on 21/4, when 
there are clear group microseisms. The groups 
are especially clear on the Z component at 
B and last approximately from 13" 21/4 to 
035 22/4. The microseisms at R are typically 
continuous, especially 18/4—20/4. On the 
other hand, during the storm on 23/4 they 
are group microseisms. 


Case 3: May 13—15, 1949 (Fig. 7) 


Weather (Free) 


A cyclone moves rapidly in the direction 
from SW towards NE straight across Iceland 
or just to the SE of Iceland. 


Amplitudes 


The amplitudes are on the whole relatively 
unimportant, even at R; only S has a more 
important microseismic storm. 

The storm at S is of short duration. An 
extremely sharp increase begins at 17 13/5; 
see Fig. 1a. Maximum is reached about 
23>—24> 13/5. The decrease is slower than 
the increase, which is a general property for 
the microseismic storms at S. The distances 
to the cyclone centre (L) and central pressures 
are as follows: 


Distance 
Time SL Ê 

Lin mb 
06h 13/5 800 990 
ogh 13/5 800 990 
12h 13/5 950 990 
15h 13/5 1,000 985 
18h 13/5 1,000 990 
03h 14/5 1,350 990 
o6h 14/5 1,600 995 


MARKUS BÂTH 


The very rapid amplitude increase at S 
obviously cannot be ascribed to some mech- 
anism localized to the centre itself. But the 
polar air outside the coast of S begins to be 
well developed about 125 13/5. The maximum 
at S can hardly depend upon standing sea 
waves. There is no coast effect; the wind 
blows parallel to the coast. The short duration 
of the storm is due to a corresponding short 
duration of the conditions favourable for 
large microseisms (the polar air stream be- 
tween S and L). An anticyclone over Green- 
land spreads out over the adjacent ocean at S 
and contributes to the disappearance of the 
microseismic storm. 

At R there is a more unimportant amplitude 
maximum at about the same time as at S 
(oo"—o2R 14/5 at R), probably also depending 
on the polar air. It is interesting to see that the 
maximum does not occur earlier when the 
cyclone is nearest to R (065 13/5), but later 
when it has moved away towards NE and 
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Mreached a position most favourable for a 
“polar air effect over the ocean to the north 
‘of Iceland. There is another also unimportant 
“maximum at R about 23h 14/5, of which there 
M is no trace at the other stations. 


The amplitudes at U and B are again very 


4 small. The coast effect is limited to the northern 


part of the west coast of Norway. This explains 
why the variations are larger at U, and that 
the amplitudes at U are larger on the N 


component than on the E component (see 


my paper 1951 b). The U amplitudes start to 


Hincrease at 185 13/5 and have maxima at 


o7® 14/5. BN has a maximum at about the 


same time. There is a cold front along the 
{ central part of the Norwegian coast at 18 13/5, 
coincident with the amplitude increase at U. 
There is no effect at the B coast in this case. 


There is a very striking similarity between 


{the variation of ocean wave heights (Fig. 4) 
§ and the amplitudes at S, but not at the other 


stations. But the rapid increase of wave 


heights occurs about 14 hours before the 


rapid amplitude rise. The maximum wave 


% height (about 18 13/5) is reached s—6 hours 
| before the microseismic amplitude maximum 
‚at S. These time differences are due to the 
4 different localities. Whatever the mechanism 
‘is, there seems to be no doubt that the mi- 
À croseisms must be ascribed to ocean surface 


Waves. 


Periods 


Ts increases considerably from 16h 13/5 on, 
and it is probable that Ts has a minimum at 
the time of the approach of the polar air. Ts is 
very much lower in this case than in all other 
cases. This may be due to the season of the 
year (annual variation of periods; see my paper 
1949, pp. 23—24). Still we have Ts > Tr 
without exception as in all other cases studied. 
The variation of Ts resembles to a certain 
degree the variations of the amplitudes at 
S. The periods Ty, Ts follow each other on 


* the whole with an obvious variation parallel 


to the variation of the amplitudes; the maxima 


| of Ty, Tg coincide with the corresponding 


amplitude maxima. The variations of the 


1 periods are much more pronounced in rela- 


tion to the corresponding amplitude variations 
at U and at B than they are at S. Also at R 


| there is a period maximum at the time of the 
| first storm, but not at the second. 
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Types of microseisms 


There are no clear group microseisms at 
any time at B and U in this case. The mi- 
croseisms are continuous also at R. The mi- 
croseisms at S have a very regular group 
character (Fig. 1 a). 


Case 4: October 1—5, 1949 (Fig. 8) 


Weather (Fig. 3. a) 

Three lows follow each other in a west- 
easterly motion from the SE coast of Green- 
land to Russia (along the northern edge of an 
anticyclone over the British Isles and Central 
Europe). The lows move all the time very 
close to the latitude 65° N. Of the three lows 
the middlemost is of greatest importance with 
regard to microseisms. The third low of the 
series is only of limited extent. The lows will 
be denoted (a), (b), and (c) in the order of 


appearance. 


Amplitudes 


The amplitudes at B and U run parallel to 
each other, whereas the variations at R and S 
deviate completely from those of U and B as 
well as from each other. The following table 
gives the essential features of the variations 


at U and B: 


Station 
and ; 
Increase at | Maximum at Remarks 
compo 
nent 
I.UE| 21h2/10 |roh—ızh 3/10] Cold front 
BE| 18h 2/10 |1oh—12h 3/10} passes over 
BN| 21h 2/10 the B coast at 
18h 2/10 
II. UE] 055 5/10 — Cold front 


touches Nor- 
wegian coast 
at Krakenes 
at 03h 5/10, 
moving to- 
wards land 


UN|o6b—o7h 5/10 


It is significant that especially in I UE reacts 
much quicker than UN. This is typical for a 
coast effect localized to the coast around B 
(due west of U). The maxima at B and U in I 
are not reached while the cyclone is deepest 
but about half a day later. At that time the 
centre is already over land, while when the 
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sharp increase starts in I, the centre was far 
outside the continental edge. The increase II 
} is not so clear at Bas at U; it is due to cyclone 
{ (c). The second storm was finished at U (N) 
@ about oob—or* 6/10. 


The amplitudes at R have two maxima, one 
about 165 2/10 and another smaller one at 


1 24" 3/10. The E amplitudes at R begin to 


increase markedly already at 225 30/9; the 
N amplitudes have a marked increase about 
12 ı/ro. The first maximum is caused by 


| cyclone (a), the second by (c). At 15" 2/10 the 


distances from the centre of (a) to R, B, U are 


| approximately the same, but the amplitudes 


at B, U are considerably lower than at R. 
A cold front passes over Iceland at 155 2/10. 


| The second maximum at R (24h 3/10) seems 


remarkable; the cyclone (c) seems to be 


{ unimportant and is situated between Iceland 
| and Greenland. It is possible that at R we have 
| instead a maximum superimposed upon the 


general decrease. If this interpretation is correct, 
this maximum would be later, i.e. about 12h 
4/10, when the cyclone is closest to R. 


The time difference of the amplitude 


| maxima at R and U, B are for cyclone (a) 


about 20 hours, for cyclone (c) about 24 hours. 
The question arises if these time differences 


| are an effect only of the different positions of 


the stations or if other factors dominate. The 


4 last explanation seems to be most probable 


for the following reasons: 


1. The rapid increases at U, B could not 
be explained only by the very gradual approach 
of a cyclone, but are connected with the 
passages of cold fronts over the Norwegian 
coast. 


2. If the mechanisms responsible for the 
microseisms were exactly the same for U, B 
as for R, we could expect similar variations of 
the periods; but there is no similarity. 

3. The maximum amplitudes occur at U, B 
when the cyclone intensity is already de- 
creasing, but not so at R in this case. 

Obviously the coast effect at the Norwegian 
coast is observable only on the continental 
side and not on Iceland or Greenland. At R 
the oncoming of polar air from due north 
with long straight paths over the ocean to the 
north of Iceland seems to be of importance 
for large microseisms. Probably the source 
of the microseisms is located over the ocean 
Tellus V (1953), 2 
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but a contribution also from a coast effect 
along the northern coast of Iceland cannot be 
excluded. 


At S there is a maximum about 06" 3/10, 
due to cyclone (a). It is remarkable that there 
is no amplitude increase at U, B before the 
very rapid increases, whereas at S there is a 
gradual rise all the time. The maximum at S 
occurs about half a day after the time of the 
largest cyclone intensity, when the cyclone 
is moving away and its intensity decreasing, 
but just when the whole ocean between 
Greenland and Scandinavia to the north of 
Iceland is occupied by polar air. The mi- 
croseisms at S have obviously no close connec- 
tion with the cyclone centre itself. There is 
no coast effect either at S. There is another 
maximum at S (E) at 11° 5/10, due to cyclone 
(c). The situation is completely the same as at 
the larger maximum. 


The variation of ocean wave heights at 
61° N, 2° E (Fig. 4) is similar to the amplitude 
variation at Sand also at U and B. The maxi- 
mum height occurs about 06! 3/10, i.e. simul- 
taneous with the amplitude maximum at S and 
slightly earlier than at U and B. Also the small 
storm towards the end of the interval occurs 
on all the records mentioned. 


Periods 


Tyand T3 follow each other, and the periods 
vary in an obvious parallelism with the ampli- 
tudes; this is valid also for the smaller storm 
towards the end of the interval. Tp and 
especially Ts are remarkably constant, except 
for a decrease from 08" 4/10 onwards; this 
may be due to a new source, cyclone (c). 
Tr and Ts show no obvious connection with 
the amplitudes. There is no increase or any 
change at all of Tr when cyclone (a) passes 
over the largest ocean depths of the whole 
region to the NE of Iceland at 12%—15" 2/10. 


Types of microseisms 


At U there are in both storms (I and ILabove) 
a certain group formation but not very 
pronounced or regular (Fig. 1b). Also at B 
there are no regular group microseisms at 
any time in case 4. The microseisms at R are 
mostly continuous. 


3 15% 163 


Case 5: March 14—18, 1950 (Fig. 9) 


Weather (Fig. 3 b) 


On the whole the weather is dominated 
by a single, large low pressure area which is 
first lying practically immovable around 
50° N, 23° W and then (from the morning of 
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17/3) moves slowly towards NE between 
Iceland and N orway. 


Amplitudes 


The amplitudes. at U, B are quite unim- 
portant as always in cases with SW winds 
over Scandinavia. The main variations are 
summarized in the following table. 
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tation 
and Increase | Maximum 


ompo + at Remarks 


WI. UE} ıgh 15/3 — 
B | 17) 15/3 


Increase sharp on 
oıh 16/3 | UE; not so clear on 


Storm I seems to be a typical case of a coast 
ffect localized to the coast around B, to judge 
rom the behaviour of the amplitudes. How- 
ver, the weather maps used show only an 
#ncreased wind velocity on the B coast from 
(st to 188 15/3 and that an occluded front 
Masses the B coast about 215—22" 15/3. It is 
remarkable that the B amplitudes later de- 
crease again although the wind velocity con- 
Sinues to increase during 16/3 and then remains 
Mt a high value (s—6 Beaufort) up to at least 
8/3. According to the experience of Topp 
ind WIEGEL (1952) it is not always that a near- 
“oastal storm is obvious from the isobaric 
Übattern of a weather map alone. 

HN The amplitudes at R show a very large 
(maximum. A very rapid increase sets in at 

br5 17/3, and maximum is reached at 05 18/3. 
Mt is remarkable that the amplitude variations 
thre completely different at R and at U, B. 
here is no trace at U, B of the large maximum 
it R in spite of the fact that the cyclone is 
Situated between Iceland and Norway and 
“he distance from R to the centre is about the 
same as the distance from B. One may believe 
‘chat there is a microseismic barrier between 
“he cyclone and Scandinavia. However, the 
study of other cases, especially case 7, makes 
chis explanation improbable (see also the 
‘discussion in Conclusions). 
| It is important to note that both the in- 
‘crease and the maximum at R do not occur 
lantil the cyclone intensity is decreasing, 
‘whereas the distance from R to the centre is 
approximately unchanged all the time. Again 
Mt seems as if a certain position of the cyclone 
were of critical importance for large micro- 
fseisms at R, just as the case is for every 
station. The critical position for R is evi- 
dently over the ocean to the E-NE of Ice- 
land. But that is just the position when the 
‘stream of polar air from due north towards 
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Iceland is most developed. This is in com- 
plete accord with most of the other cases. It 
could be believed that the delay of the micro- 
seisms behind the largest cyclone intensity was 
only due to the time it takes for the ocean 
surface waves to reach their full development. 
This circumstance certainly contributes but if 
it were the only explanation, large microseisms 
should be obtained also for other wind direc- 
tions, if they have been blowing long enough 
over a sufficiently large body of water. How- 
ever, this is definitely not the case. For in- 
stance the microseisms in Scandinavia are never 
large when the winds are coming from SW 
or W. It is only in connection with polar air 
from NW—N that large microseisms are 
obtained. And this is a common property for 
all four stations investigated. 

The amplitudes at S vary on the whole, but 
not in detail, in parallelism with those at R. 
The maxima are approximately simultaneous. 
The variations at S are more gradual. The 
microseisms at S are most probably due to 
the polar air over the ocean to the west of the 
cyclone. 

A question of great importance is how the 
microseisms originate in the polar air. The 
two main theories will be considered. 

1. Standing sea waves (LONGUET-HIGGINs, 
1950). They could be expected at a coast but 
there seems to be no possibility for their 
existence over the open ocean in case 5. The 
wind direction has been practically unchanged 
all the time, as the cyclone has moved only 
very slowly. 

2. Organ-pipe theory (Press and Ewıng, 
1948, and later personal communication). If 
the source for the microseisms at R and S 
were the same, the periods at R and S should 
vary in parallelism, the only difference being 
due to possibly different distance to the source. 
But the periods Tr and Ts behave quite 
differently. 


Periods 


Ts has a pronounced minimum at 161—17h 
15/3, i.e. at the same time as the amplitudes 
begin to increase. This is a further indication 
that the source of this storm is on the B coast. 

The periods at U, B diminish from or? to 
toh 17/3, about 0.4 sec at U and about 0.8 
sec at B. There is a marked increase of 0.8 sec 
of Ty at or® 18/3. These cases are examples of 
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Fig. 10. Case 6: March 24—28, 1950. 


period variations without corresponding am- 
plitude variations. The reason is most probably 
that the position of the dominant source is 
variable without producing changes of ampli- 
tudes. For instance it is not excluded that the 
higher Ty from orb 18/3 on could be due to 
an origin from the cold air sector of the 
cyclone. It is known that U sometimes records 
microseisms with small amplitude (< 1 4) but 
large period (8—9 sec). In addition to the cases 
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mentioned in (1949) the following are typical: 
January 12, 1951; September 12, 1951; Decem- 
ber 22—23, 1951; January 7, 1952. These 
microseisms may be due to a source further 
out on the open ocean and not at the Nor- fl 
wegian coast (for a fuller discussion see my | 
paper 1949, pp. 140—141). The microseisms | 
from the Atlantic have always small amplitudes |. 
ac" I 

As in every other case Ts is larger than the H 
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beriods at the other three stations. In both this 
“fase and case 4, Tr is on the whole larger 
: han Tu and Des 


Types of microseisms 


The microseisms at U are continuous. In the 
orning of 18/3 at the time of the period 
crease, they become irregular, especially on 
he E component. There is no clear group 
Hormation at any time at B either. The same 
js true at R, also for the storm towards the end. 
At S the group microseisms (on 16/3) are not 
30 regular as in many other cases. 


Case 6: March 24—28, 1950 (Fig. 10) 


Weather (Fig. 3 c) 
|} During the interval investigated a cyclone 
‘moves slowly from south Greenland via Ice- 
and towards the region to the north of 
thorthern Norway. 


Amplitudes 


M The main amplitude variations at U, B are 
‘as follows. 


©) Station | Increase at | Remarks 


I. U |2ıh 25/3—17h 26/3 
B {around o2h 26/3 


les OF 10h 27/3 Cold front over cen- 
tral and northern 
part of Norwegian 
coast 
| III. UN 2th 27/3 


7 


" It is significant that in II the increase is much 
…learer at U than at B; this is characteristic for 
‘1 coast effect localized to the central and 
thorthern parts of the Norwegian coast (see my 
‘paper 1951 b). The same is true for storm I. 
"A cold front passes U about 20" 27/3; if this 
front has something to do with the approxi- 
‘nately simultaneous increase remains un- 
lertain, because if this were true, it is a rather 
linique case. The amplitudes at U, B have no 
‘maxima during the interval investigated; they 
lo not return to the original values, whereas 
‘he amplitudes at R, S do. The amplitudes 
ht U, B are on the whole relatively unimpor- 
sant. The wind velocity at and outside the 
Norwegian coast is also small. The cyclone is 
hot very intense, except for a shorter time 
on 27/3. 
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The amplitudes at R start to increase at 
20" 25/3 and reach their maximum at o9"—r0h 
26/3. It is obvious that the increase first starts 
when the cyclone is situated over the ocean to 
the NE of Iceland and not when it is nearest 
to Iceland. The reason cannot be varying 
cyclone intensity, as this is practically un- 
changed the whole time. It is clear that there 
is no close relation between the microseisms 
and the cyclone centre itself. The large 
microseisms at R are obtained when the cyclone 
has reached such a position that the polar air 
stream north of Iceland is most fully developed. 
From about 12h 26/3 the area of polar air moves 
away from Iceland and an anticyclone south 
of Iceland increases in importance; the mi- 
croseisms at R are then decreasing. 

The amplitudes at S are on the whole in good 
parallelism with those at R. There is a very 
marked increase at S at 161—18h 25/3, ive. 
somewhat earlier than at R, and the maximum 
is about 10" 26/3, i.e. about the same time as 
at R. It is evident that the marked increase at S 
occurs when the wind changes to a northerly 
direction and the polar air is coming down 
over the ocean to the east of S. The increase 
cannot be explained by the distance to the 
centre, which is about the same all the time 
before maximum, nor can it be explained by 
variations of cyclone intensity. There is no 
coast effect at S; the wind is parallel to the 
coast at S during the amplitude maximum. 

After the maxima the amplitudes at R, S 
decrease slowly, but at U, B the amplitudes 
are unchanged or increasing. The distance 
from the centre to U, B and R are about 
equal, while the distance to S is smaller. The 
stations U, B are not nearer to the cold air 
sector than are R, S. In spite of this the ampli- 
tudes decrease at R, S, but not at U, B. This 
is an indication of the importance of the 
Norwegian coast for Scandinavian micro- 
seisms. 

At the beginning of this case there is a 
cyclone around southeast Greenland with on- 
shore winds to the north of the centre. There 
is no microseismic effect at S, i.e. no coast 
effect. The on-shore winds are not polar air, 
but more stable, warm air. 


Periods 


There are clear minima of Ty at 23h 25/3— 
ot) 26/3 and at 19 27/3, in both cases at the 
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same time as the amplitudes increase or begin 
to increase (I and II above). Ty is larger than 
in the two preceding cases. 

Ts has a minimum about 185—204 26/3; a 
cold front passes B about 16 26/3 from NW 
to SE. After the minimum Ty increases all the 
time towards the end of the interval. The cold 
front at the time of minimum 73 had no 
observable effect on amplitudes at B or on 
amplitudes or periods at U (compare my paper 
1951 a). The connection between period 
minima and cold fronts is indisputable, and is 
a further support for the microseismic im- 
portance of polar air. 

TR increases up to the time of maximum 
amplitudes from 4 to 5 % sec and then remains 
approximately constant. There is some simi- 
larity between this variation and the variation 
of the periods of the ocean surface waves 
(Fig. 4). 

Ts decreases from the beginning up to 
18"—20" 25/3 and then runs parallel to Tr. The 
behaviour of Tr and Ts resembles very much 
case 2. Also for Ts the minimum is coincident 
with the largest amplitude increase. 


Types of microseisms 

The microseisms have a typically continuous 
character at U, B, and R. At S they have 
generally a regular group appearance. 


Case 7: October 7—11, 1950 (Fig. 11) 

Weather (Fig. 3 d) 

A whole series of cyclones moves rapidly 
towards NE between Iceland and the British 
Isles. From the microseismic point of view the 
weather development is more complicated 
than in all the preceding cases. For ease of 
reference the cyclones will be denoted (a), 
(b), (c), (d) in the order they appear. For more 


exact definition the following data are given: 


Cyclone Time Position of centre 
ee 
(a) 06h 7/10 TAN IT 
(b) 06h 8/10 SINE IE 

(c) 06h 9/10 58° N, 25° W 
(d) 06h 11/10 SZENEN: 
Amplitudes 


Due to the complicated weather develop- 
ment the amplitude variations are more 
difficult to explain than in the preceding cases. 
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The R amplitudes are at first decreasing 
from a maximum on 6/10—7/10 (RE maxi- 
mum about 142—24> 7/10). At o6! 6/ıo 
there is a cyclone (a) with a central pressure 
of 975 mb NE of Iceland and polar air is 
coming from the north over Iceland. This is 
the only obvious reason for this microseismic 
storm at R. The maximum does not occur 
earlier when the cyclone was much closer to 
Iceland and its intensity larger, but first when 
the critical position NE of Iceland has been 
reached. The next cyclone (b) has no greater 
effect at R, the probable reason being that the 
cyclone, when the critical position is reached, 
is only of limited extent and produces no 
important inflow of polar air. Polar air is 
coming in later, but then the whole system 
has moved away too much to be of greater 
significance, and also the cyclone (c) disturbs 
the conditions around Iceland. At R there is 
only a small maximum around 20" 9/10, 
superimposed on the general decrease. Cyclone’ 
(d) is the probable reason for an amplitude 
increase at R in the morning of 12/10 (cyclone 


effect). 


The amplitudes at S seem to be simpler to 
explain, as only cyclone (b) is of importance. 
A very rapid increase starts at 145—16> 8/10 
(before that the amplitudes were slightly 
decreasing from a maximum corresponding to 
cyclone (a); maximum is reached about 
07" 9/10. As in general at S the increase is 
much more rapid than the decrease. Some 
details of the development are given in the 
following table. 


Dis- 

; tance | L Distribution of | 

AIR { 
aay SL | mb polar air 


ooh 8/10 | 1,100] 955| no 

06h 8/10 | 1,100] 955] no 

12h 8/10 | 1,100 slight beginning over the 
ocean outside Greenland: 


ooh 9/10 | 1,050 | 975] well developed off the 
COASt Ot —S: 

06h 9/10 | 1,000 | 975 very well developed off S: 

12h 9/10 900 | 975] well developed; 

ooh 10/10} 1,000 | 980 decreasing 


| 

The microseismic storm at S cannot simply | 
be explained by variations of the distance or: 
the intensity of cyclone (b). Again clearly the > 
polar air is of decisive importance. But there: 
Tellus V (1953), 2 M 
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is no coast effect at S, as the wind is approxi- 
mately parallel to the coast. At least in win- 
tertime also ice along the Greenland coast 
prevents a coast effect. 

The variation of the amplitudes at U, B 
deviates completely from S. The cyclone (b) 
which was of great importance at S is of very 
little significance in Scandinavia. It gives SW 
winds over Scandinavia. The amplitudes at 
U, B indicate a clear coast effect localized to 
the B coast: 


1. The amplitudes at B are considerably 
higher than at U. 
2. Ag > AN without exception at U. 


3. The increase is earlier and more rapid for 
Ar than for Ay at U. 
There are two storms (I and I) at U and B. 


Station 
and 2 
Increase at Maximum at | Remarks 
compo- 
nent 
Ir ce 24h 7/10 05h 8/10 cyclone 
UN orth 8/10 (b); cold 
BE) beginning of fronts 
interval pass the 
(21h 7/10) 06h—o7h 8/10 | B coast 
BN 24h 7/10 about 
152776 
and 
24h 7/10; 
II. UEJabout 13h 1o/1olabout o6h 11/10 cyclone 
UN 16h 10/10 (c) 
BE 13h 10/10 o4h—o6h 11/10 
| BN not clear 


| 


At oo" 8/10 the cyclone centre (b) lies in the 
direction W 40° N from U, but the E com- 
ponent of the amplitudes is considerably larger 
than the N component. This happens always 
at U when cold fronts pass the B coast. This 
excludes the cyclone centre as source of the 
microseisms. The question remains if the 
source is located at the B coast or further out 
on the open ocean within the polar air. The 
last-mentioned possibility is excluded for the 
reason that the storms I and II are observed 
only in Scandinavia, and not at R, S. This is 
easily explained for a source at the B coast, 
as the distance to R is then too great. But it 
can hardly be explained by an origin within 
the polar air over the ocean, as the distance 
from R to the polar air sector is just about the 
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same as from U, B. On the other hand, for 
R, S the polar air over the open ocean, espe- 
cially with straight wind paths over a large 
body of water (large fetch), is of greatest 
microseismic importance. The special réle 
of the coast will be discussed in the Con- 
clusions below. 

At oo! 8/10 there is a cyclone (b) between 
Iceland and Norway; there is no microseismic 
storm at R, but there is one at U, B. It is 
interesting to compare this situation with 
case 5, where there was also a cyclone at 
about the same place between Iceland and 
Norway, but accompanied by a microseismic 
storm at R and not at U, B, i.e. just the reversed 
situation. The difference lies in the different 
form of the isobars. In case 5 there were SW 
winds over Scandinavia and polar air only 
over Iceland and over a large body of water 
to the north of Iceland. On the other hand, 
in case 7 polar air is coming in over the B 
coast, but there is no pronounced inflow of 
polar air north of Iceland, at any rate only a 
small fetch. That completely different micro- 
seismic situations are obtained for a cyclone 
with its centre in about the same place shows 
clearly that the source cannot be located at 
the cyclone centre. The different conditions 
could only be explained from differences in 
the distribution of different air masses. We also 
note that the idea of a source in the centre and 
a microseismic barrier, e.g. along the conti- 
nental edge outside N orway, is untenable. The 
cyclone centre (b) in case 7 is far outside the 
continental edge at oo! 8/10, and the idea 
mentioned could not explain why in one case 
we get a microseismic storm only in Scandi- 
navia, in another case only on Iceland. The 
conditions are simply explained by the distri- 
bution of the polar air without using the 
hypothetical explanation by means of a 
microseismic barrier, the existence of which 
is by no means proved. 

That Ag > Ay at U at 068 9/10 as all the 
time can only be explained by a coast effect 
at the B coast. At 065 9/10 the centre (b) lies 
in the direction W 60° N from U and the 
polar air over the ocean to the W of (b) is in | 
approximately the same direction. The coast | 
effect is limited to the B coast. | 

There is some similarity only between the : 
microseisms at S and the ocean wave heights ; 
(Fig. 4). The increase and the maximum of the : 
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waves occur earlier than of the microseisms 
@at S. 


Periods 


Ts is as always remarkably constant (= 6 
sec) and has no variation in parallelism with 
the amplitudes. The period variations are 
much larger at R, U, B. 

J Ts, Tu, especially Ty, show some remark- 
gable variations. Ty has one maximum at 
about o15 9/10 (= 6.4 sec) and another 
“maximum at o4 11/10 (= 5.6 sec) with a 
minimum between, about 065 10/10 (= 4.8 
{sec). The period variations are quite consider- 
‚able and no doubt significant. Tp deviates 
“from Ty in not showing the first maximum. 
There is a pronounced parallelism of Ty, Ts 
with the amplitudes only in the second maxi- 
mum. We have to observe that the amplitudes 
UE and UN behave rather differently, and 
that the parallelism between Ty and UN is 
“relatively good on the whole, but for UE it 
is good only for the last maximum. The 
amplitudes BN, BE, and UE are quite similar 
(due to the B coast), whereas UN deviates. 
A UN is probably more dependent on conditions 
© at the more northerly parts of the Norwegian 
M coast, i.e. a more distant source and higher 
‘periods. This distance effect is certainly the 

essential reason that Ty > Tg during the 

first part of the interval. A comparison of 
Ty and Tg for U clearly shows that during 

8/10 Ty > Tr, whereas for 9/10—11/10 they 
‚are approximately equal. The minimum 
| periods occur just when a new source, cyclone 

(c), comes to importance. The variations of 
Tu in case 7 cannot be explained by motion 
1 of the cyclone centre (b) or of the polar air 
| over varying ocean depths. 


| Types of microseisms 

(| The microseisms at U are irregular on 9/10, 
especially on the N component. The micro- 
. seisms are continuous at U, B and also at R. 


| Conclusions 
fH In the study of the special cases several 
questions have appeared which cannot be 
answered definitely due to the lack of sufficient 
| meteorological and especially oceanographic 
data. But in addition a number of well- 
established general conclusions can be drawn. 
, They will be summarized in this chapter. 
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1. Importance of polar air 


For every station investigated there is a 
certain critical position in which a cyclone of 
given intensity produces the maximum micro- 
seismic amplitudes. I think that quite some 
progress towards solving the microseismic 
problems could be made if these positions were 
determined for every station within a larger 
area and also the situations studied which 
produce rapid amplitude increases. For Upp- 
sala the maximum microseisms are obtained 
for a low with its centre in the region of 
northern Russia (see BATH 1949); for Bergen 
the critical position is outside the central part 
of the Norwegian coast (on-shore winds and 
waves at the coast around Bergen; see BATH 
1951 b). For both stations rapid amplitude 
increases are generally obtained when cold 
fronts pass the Norwegian coast in the direc- 
tion sea to land (see BATH 1951 a). In the 
present paper it has been shown that maximum 
amplitudes occur at Reykjavik when a low is 
over the ocean to the NE of Iceland and polar 
air is coming down from the north across 
Iceland along straight paths over a large body 
of water; and the critical position for Scoresby- 
Sund is over the ocean to the east of the sta- 
tion and when the area between Scoresby- 
Sund and the centre is occupied by polar air. 
All four stations have in common that the 
microseisms show no close relation to the 
cyclone centres themselves, but that they 
are at all stations closely related to the approach 
of polar, unstable air. This circumstance cannot 
simply be explained as a time lag between 
winds and ocean waves, but must be due to 
some special property of the polar air (see 
discussion of case 5). In general it seems diffi- 
cult to explain this as an effect of standing 
ocean waves, especially when the wind has 
been blowing in about the same direction for 
several days and not against a steep coast. 

Indications of some kind of cyclone effect 
with a closer connection to the cyclone itself 
has earlier been obtained for Uppsala in cases 
with a minimum coast effect (BATH, 1949). 
At Reykjavik a cyclone effect, possibly due 
to standing ocean waves, is obtained when a 
rapidly moving low is in the near vicinity of 
Iceland. 

In several cases there are greater similarities 
between the ocean wave observations at 
61° N, 2° E and the microseisms at Scoresby- 
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Sund than at any other of the stations investi- 
gated. The reason may be that the ocean wave 
observations correspond better to the condi- 
tions on the open ocean, whereas the micro- 
seisms at least in Scandinavia are dominated 
by a coast effect. There seems to be no doubt 
that the ocean surface waves are of essential 
importance, whatever the mechanism is. 


2. Coast effect 


Recognizing the common microseismic 
importance of the polar air for Scoresby-Sund, 
Reykjavik, Bergen, and Uppsala, the question 
arises what importance the Norwegian coast 
has in the case of Scandinavia, in other words 
what is the nature of the coast effect. There is 
no coast effect in Scoresby-Sund in the cases 
studied; at Reykjavik the coast effect may 
contribute, but for both stations the main 
source is probably located over the open 
ocean within the polar air. Some kind of 
coast effect on the Norwegian coast, especially 
within the polar air, is the main source of 
microseisms in Scandinavia. For a full dis- 
cussion of the evidences for this, see BATH 
(1951 c). 

The coast effect at a steep coast could be 
due to one or both of the following possi- 
bilities. 

Pressure fluctuations on the ocean bottom 
below standing sea waves, formed by reflec- 
tion from a steep coast (LONGUET-HIGGINs, 
1950). In this case standing sea waves can 
naturally be expected. This idea is supported 
by the behaviour of the periods (see below). 
The special importance of the polar wind 
would then be its property of producing 
higher ocean waves than within warm, stable 
air of the same wind velocity (see BATH, 1951 a, 
pp. 295—303). Concerning ocean wave phe- 
nomena of importance for formation of stand- 
ing waves at a coast, reference is made to an 
investigation by WILLIAMS and Isaacs ( 1952). 

The importance of surf at a steep coast still 
deserves further attention. 

The continental layers serve as a wave- 
guide for the microseismic waves, whereas 
the ocean bottom does not. It was shown by 
GUTENBERG (1932) that the microseisms origi- 
nating at the Norwegian coast are observed in 
Russia and far into Siberia, whereas in this 
paper it is shown that even at such relatively 
near places as Reykjavik and Bergen the 
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microseisms behave quite differently. Ewinc 
(1951) found from theoretical and empirical 
considerations that Rayleigh waves of the 
same period as the microseisms could not 
travel very far over an oceanic structure. 
He proposed the hypothesis that microseisms 
instead were body waves (P and SV), for 
which different absorption would not be 
expected. However, the empirical result, 
obtained here, does not seem to require the 
introduction of body waves, but could be 
explained by the hypothesis of Rayleigh waves 
or a combination of these and other surface 
waves. At any rate, if the microseisms were 
only body waves, I cannot see how the range 
could be so much larger over the continent 
than over the ocean. 

When this result had already been found 
for our region, I found that CARDER (1952) 
had got the same results for other parts of the 
Atlantic Ocean and for parts of the Pacific 
Ocean. The condition may be general and, 
as stated by CaRDER, microseisms which are 
well recorded on land would hardly be 
expected to originate from great distances at 
sea. It is not excluded that the continental 
channel is the same as for the short-period 
surface waves Lg and Rg, observed in earth- 
quake records when the path is purely con- 
tinental. See also a recent paper by GUTEN- 
BERG (1951). 


3. Periods of microseisms 


On the whole, it is more difficult to explain 
the observed period variations than the ampli- 
tude variations. But for every theory of 
microseisms it is essential that also the be- 
haviour of the periods can be explained. The 
period at Scoresby-Sund is remarkably con- 
stant and shows no obvious connections with 
the amplitude variations. The period varia- 
tions are larger at Reykjavik but also here there 
are in general no very obvious connections 
between the variations of period and ampli- 
tude. At Bergen and Uppsala, on the other 
hand, there are often large period variations, 
and the periods vary in obvious parallelism 
with the amplitudes. In addition the periods 
vary with the distance to the source of micro- 
seisms (for a full discussion of this problem 
for the Scandinavian region see BATH, 1952 b). 

It is characteristic with a period minimum 
coincident with a rapid amplitude increase 
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and the passage of a cold front over the Nor- 
wegian coast around Bergen. Similar observa- 
tions have been made in New Zealand by 
JONES (1949). There is a very striking resem- 

lance between this behaviour of the micro- 
seisms and the behaviour of periods and ampli- 
tudes in near-coastal storms, the latter dem- 
onstrated by Topp and Wircet (1952); see 
especially fig. 3 in that paper. They found 
that in near-coastal storms, the periods of the 
ocean waves dropped rapidly while breaker 
heights increased. There is no doubt that the 
microseisms at Bergen and Uppsala in these 
cases are due to ocean surface waves in near- 
coastal storms. Similar period minima are 
observed at Scoresby-Sund, apparently simul- 
taneous with the approach of polar air over 
the ocean outside Scoresby-Sund. 

Apparently due to varying distance to the 
source there are period variations sometimes 
observed at Bergen and Uppsala, with no 
simultaneous amplitude variations. 

If the standing ocean wave theory is valid, 
we should expect a parallelism between 
periods and amplitudes, for the organ-pipe 
theory not. As mentioned above, standing 
waves may be expected at the Norwegian 
coast, but often not over the open ocean in 
the cases studied. This means that standing 
ocean waves may be of importance for Scandi- 
navian microseisms, and some other mech- 
anism is the main reason for the microseisms 
on Iceland and Greenland. However, I have 
in no case succeeded in correlating the period 
variations at any place with varying depths of 
ocean and underlying sediments, as require 
by the organ-pipe theory. The ocean depths 
are well known (see Stocks, 1950), whereas 
there are no measurements at hand of sedi- 
ment thicknesses so far to the north as our 
region of the Atlantic Ocean. Some informa- 
tion for regions more to the south is obtained 
from results of the Albatross expedition, as 
well as from British and American investiga- 
tions. The standing wave theory requires a 
ratio of 0.5 between the periods of microseisms 
and of ocean waves. Using the ocean wave 
observations at 61° N, 2° E, we found this 
ratio to be significantly greater than 0.5 at all 
statiors investigated. 

The periods at Scoresby-Sund are practically 
with« ut exception larger than at the other 
three stations. 
mellus iV (1953),72 
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4. Microseismic barriers 


The fact that cyclones in the Atlantic are of 
relatively little microseismic importance in 
Scandinavia, whereas the effect on Iceland 
and Greenland may be large, could lead us 
to suspect a microseismic barrier west of 
Scandinavia. In a previous paper (BATH, 
1952 a) I mentioned that no barrier effect was 
observed at the continental edge. It is also 
clear that it is not a sufficient reason for large 
microseisms that a cyclone has passed the con- 
tinental edge. For instance a cyclone around 
north Russia has a far greater importance for 
Scandinavian microseisms than a cyclone 
situated over the North Sea, also within the 
continental edge, and moreover at a much 
shorter distance from Uppsala. Large increases 
occur in Scandinavia in several of the cases 
studied already when the cyclone is far outside 
the continental edge. The reversed behaviours 
of the microseisms in cases 5 and 7 (see discus- 
sion of these cases) also disproves the existence 
of a microseismic barrier. The microseismic 
amplitude variations at the different stations 
in our region can be explained by the distribu- 
tion and approach of polar air without using 
the artificial explanation by means of micro- 
seismic barriers. 

On the other hand, the apparently low 
ability of oceanic structures to transmit micro- 
seismic waves may be called a barrier effect. 
I have earlier offered a supplementary hy- 
pothesis (BATH, 1952 a) for the explanation 
of barriers. In addition to this hypothesis I 
would like to mention also the possibility 
that irregular distribution of microseisms 
(barrier effects) may be due to interchanging 
oceanic and continental structures. 


5. Types of microseisms 


The microseisms at Scoresby-Sund are in 
general much more regular than at the other 
three stations and show a clear group character. 
At Reykjavik as well as at Bergen and Uppsala 
the microseisms are usually continuous. Group 
microseisms are observed at Bergen and Upp- 
sala in some cases with a coast effect limited to 
the southwest coast of Norway (around 
Bergen). In general the microseisms are likely 
to be more regular for a limited source, whereas 
they become irregular when several sources 
are active simultaneously or when there is 
one source of large extension. 
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Experimental Determination of Electron Orbits in the Field 


of a Magnetic Dipole 


By ERNST-ÂKE BRUNBERG, The Royal Institute of Technology, Stockholm 


(Manuscript received February 18, 1952) 


Abstract 


An analogue method for determination ot electron orbits in a magnetic dipole field has been 
developed by Malmtors several years ago. The method is in principle a scale model experiment 
in a vacuum chamber where the results are obtained by measuring the deflection of an electron 
beam in the field of a “‘terrella’’. During the last years this method has been further developed, 
and particle orbits corresponding to the following momenta have now been measured: 


I) 2+ 10° — 10! eV/c 
Orbits arriving at magnetic latitude 58° in E, S, W and N directions at zenith angles of 0°, 
20°, 40°, 60°, 80°. Orbits arriving from zenith at different latitudes beginning with 38° and 


continuing at intervals of 2 degrees up to the pole. 


2) 101° — oo eV/c 
Orbits arriving at magnetic latitude 
0°, 10°, 20°, 30°, 40°, 50°, 60°, 70°, 


80°, and 90° 


from E, S, W and N directions at zenith angles of 
070327.0710096212:032, 7402.,432., 50 004, 722, Fans On 


The direction of orbits far away from the dipole is published as a function of momentum. 
The particles are assumed to move towards the earth and to be positively charged. 


This article includes only diagrams for the latitudes 90°—60°. The remaining ones will 


be published in the next issue of this journal. 


It is of great importance in cosmic physics 
to know the orbits of charged particles in a 
magnetic dipole field. This work aims to give 
a general survey of particle orbits in the entire 
range of momentum above 101° eV/c and thus 
continues the determination of orbits in the 
range 10°—101° eV/c, that have been carried 
out earlier in the same way by MALMFORS 
(1945). 

The problem of the paths of charged par- 
ticles in a magnetic dipole field has been tackled 
mathematically by STÖRMER (1904—1937), 
LEMAITRE and VALLARTA (1936) and many 
others (cf. HEISENBERG, 1943, and JANOSSY, 
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1948) in connection with the theories of the 
aurora and the cosmic radiation. However, 
a comprehensive survey of the most impor- 
tant part of the energy range is still lacking. 


The principle of the experiment 


The method is to shoot out an electron 
beam from the surface of a “terrella” and 
then measure the deflection of the beam in the 
magnetic dipole field. 

In this way it is possible to trace back a 
cosmic particle orbit and find its direction 
before being deflected in the earth’s magnetic 


field. 
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The experiment was carried out in a Cy- 
lindrical vacuum chamber made from brass, 
730 mm in length with a diameter of 400 mm. 
The ends of the chamber were covered with 
thick glass plates to look through. The terrella, 
with electron gun, was suspended inside the 
chamber and 130 mm from one end. As is 
further explained below, the electrons could 
be shot out in any direction and from any 
latitude of the terrella. At the other end of 
the vacuum chamber the direction of the 
electron beam was measured. In this end the 
dipole field is practically zero, and the direction 
of a particle orbit here is nearly the same as at 
an infinite distance. 

It is evident that there must exist a certain 
relation between the energies, dipole fields 
and geometric dimensions in cosmic ray case 
and for the model. In order to relate the 
model to this case the required conditions are 
set out below. MKSA units are used. 


Cosmic rays Model 
Radius of the earth..... R Re 
Magnetic dipole moment M Mn 
Energy of particles .. ...- V Ve 


Uniformity between the dimensions of the 
systems requires 


ee (1) 


The relativistic formula for the movement of a 
single-charged particle in a magnetic field is: 


Bo= ple (2) 

where p = - V(eV)? +2eVmc? (3) 
with B = magnetic field strength 

p = momentum of the particle 

o =radius of the curvature of the 

particle path 

eV = particle energy 

My = rest mass of the particle 

e = charge of the particle 


= velocity of light 


It is obvious that when considering cosmic 
particles a relativistic treatment of their motion 
is necessary. For model particles, however, 


eV < mec? 
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and thus equation (2) is reduced in the experi- 
ment to 


(B+ 0)m = 3-37 10° VVm (4) 


For the dipole field we have 


= le ee 
Rie 


em 


(5) 


R3 m 


After combining equations (1), (2), (3), (4), and 
(s) the following final expression 1s obtaine 
with constants inserted 


(R = 6.37" 108m, M = 3171-10 Aan) 


SR 
= : 10% ae 
p = 202- 101°. V/V, M, (6) 
For the model we have 
= 8.0 - 10? metres Vin = 150 volts. 


Rn 
My = 0—15 amp.metre? 

In the expression (6) p is in eV/c. This unit 
will be used in the diagrams (1—38) and if 
the particle energy is V electronvolts, p in 
these units is 


p=VVırzmelV  eVlc 


Apparatus 


In the general principles the apparatus is 
the same as used by Malmfors (see fig. 1). 
The construction of the dipole and the electron 
gun is, however, quite different. 


$ 
+ | 
* 


Fig. 1. View of the apparatus. The terrella can be seen 
in the nearest part of the vacuum chamber. 
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Generation of the magnetic dipole field 


Malmfors used a magnet cast in A/-Ni-Co 
steel in the form of a sphere, magnetised in a 
homogeneous magnetic field. However, it is 
difficult to cast such a sphere, and irregularities 
in the metal will result in deviations from the 
exact dipole field. Measurements on the field 
of a cast magnet showed that the axis of the 
dipole was situated about r mm from the 
centre of the sphere. 


The dipole field is now generated by a 
coil inside the hollow terrella; thereby, the 
above mentioned error is only 0.2 mm. The 
possibility of varying the strength of the mag- 
netic field is another advantage. Eq. (6) shows 
that the equivalent cosmic ray momentum is a 
function of the magnetic dipole moment of 
the model, and by varying the current through 
the coil, the desired momentum is obtained. 
This could also be carried out by varying 
the voltage of the electron gun, but it is 
preferable to have a constant voltage on the 
gun to facilitate the focusing of the electron 
beam. Finally, it is a great advantage to be 
able to adjust the beam without the presence 
of the magnetic dipole field. 


It can be shown (HARNWELL, 1938) that a 
spherical coil generates a magnetic dipole 
field when wound in such a manner that the 
current density along an arc in the meridian 
plane is proportional to the sine of the polar 
distance. 


A solid coil could be built up of many con- 
centric spheres and a magnetic dipole M ob- 
tained by superposition. It is easier, however, 
to use a number of coaxial cylindrical coils 
and a suitable design has been given by BROWN 
and SWEER (1945). 

Fig. 2 shows a crossection of the construction. 
Insulated copper wire (a) was wound on a 
brass core (b). The layers were separated from 
each other with paper and bakelite varnish. 
The whole coil was enclosed in a copper shell 
(c) and baked at a temperature of + 120 °C. 
Afterwards, the copper shell was evacuated 
to a pressure of 1 w. The evacuation was 
necessary to prevent possible bursting of the 
shell when placed in a vacuum chamber. 
The diameter of the wire was 0.6 mm and 
the entire resistance 50 2. 

The measurements of the generated magnetic 
field were carried out by means of a rotating 
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Fig. 2. Construction of the terrella. (Dimensions in 
millimeters.) 


fluxmeter coil of sufficiently small dimensions. 
The induced voltage was measured with a 
vacuum tube voltmeter. The variation of the 
field along the arc from pole to pole was 
compared with an ideal dipole field; the devia- 
tion was smaller than the probable average 
error of measurements (about 2%). The 
same range of deviation was obtained when 
measuring the field along a radius and com- 
paring it with the r-3-curve. 
The dipole was calculated to 


M = 10.73 amp.m?/amp. magnetizing as 
7 


It was necessary to compensate the earth’s 
magnetic and other disturbing fields. Good 
compensation was achieved by placing current 
coils on the walls, ceiling and floor of the 
room. The magnetic field thus left in the 
vacuum chamber was less than 0.03 gauss. 

Instruments were placed at a safe distance 
from the vacuum chamber or were magnet- 
ically screened. 


The electron gun 


The electron gun required an electrode 
system with small dimensions to give a well 
defined electron beam with a diameter of 
1—2 cm after more than half a metre travel. 
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Fig. 3. Construction of the electron gun. (Dimensions 
in millimeters.) 


Fig. 4. The electron gun in its mounting. 
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This short electrode system was necessary, 
otherwise the electrons would be appreciably 
deflected by the magnetic field before leaving 
the gun and lead to an error in the direction 
in which the electrons are shot out. The gun 
meets the specified requirements except in 
some few directions, perpendicular or nearly 
perpendicular to the magnetic lines of force. 
Especially in the lower momentum range, 
this deviation can be noticeable as discussed 
more in detail in the next section. 

The experiment was carried out in a va- 
cuum of 1075 mm Hg. 

Fig. 3 shows a drawing of the electron gun. 
The body is made of ceramic material (“Pyro- 
phyllite”) which can be machined, but hardens 
after firing. The electrode system is made of 
molybdenum plates. The cathode is a nickel 
wire with a diameter of 1 mm coated on the 
end surface. It can easily be replaced when the 
emission has been reduced through poisoning 
or other destruction of the carbonate film. 
The cathode is indirectly heated by a tungsten 
spiral, the last turns of which are reversed to 
avoid any disturbing magnetic field inside the 
electrode system. Fig. 4 gives a picture of the 
gun in its mounting. 

Fig. 5 is a diagram showing the principle 
of fixing the terrella and electron gun in the 
vacuum chamber. By means of the screw (S), 
the model earth with mounting and gun 
can be moved around the trunnions T—T. 
The gun in its mounting can be turned in the 
horizontal plane with respect to the terrella 
by means of vacuum-tight cones (C) and 
gears (G). For each turn the gun is moved 
through 8° in zenith angle by means of a 
cross and bevel gear drive (see fig. 6). 

Finally it is possible, from outside the 
vacuum chamber, to alter the latitude of 
the gun on the model by means of a friction 
drive between the model and the gun mount- 
ing. This arrangement is not shown on the 


figure. 


Determination of direction 


‚Fig. 7 shows the method for determining the 
direction of the electron beam as used by Malm- 
fors. For a given energy and direction of emis- 


Fig. 5. Construction of the system carrying the terrella 
and the electron gun. 
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sion, a certain path is obtained. The model 
(with gun) is turned until the beam strikes a 
plate with specially shaped teeth. The shadow 
of some tooth can be seen in the spot on the 
fluorescent screen. The quantities m, n and | 
give then the direction of the electron beam 
at a long distance from the dipole. The differ- 
ence between its direction here and at an 
infinite distance due to the deflection in the 
fringe of the dipole field is only one or two 
degrees, and is determined by an approximate 
method worked out by Malmfors (cf. page 
24). 

The results of measurements are given in 
the form of two angles (see fig. 8) Oy, the 
north latitude angle and Wg, a longitude angle, 
reckoned eastwards from the meridional plane 
through the electron gun. These angles cannot 
be directly read in the measurements but are 
calculated from the measured direction on the 
screen. 

The error in the initial direction of the 
electrons, which is due to the deflection of 
the electron beam inside the electron gun, 
has been estimated in the following manner. 

According to Malmfors, the deflection can 
be approximately determined from 


25 fit, D 


+ sin y (8) 


Op = —— 
: Soe If Fe VV 


where B = magnetic field strength 


257-100 VV, = Bo 

a, = distance between the cathode and the 

anode 

acceleration voltage 

=the angle between the electron gun 
and the magnetic field. 


Vin 


It is assumed that B is a constant inside the 
electrode system, and that the electric field 
strength is constant along the axis of the system 
between the cathode and the anode. 

As ad) = 0.46 cm and V,, = 150 volts we 
obtain in degrees w,° = 0.13 : B sin y. 

In the present electrode system, however, 
the electric field is weaker in the vicinity of 
the cathode and then increases towards the 
anode. Thus the error ought to be greater 
than stated above. + 

A graphic determination of the error inside 
the electrode system sustained these assump- 
tions. The error was estimated to be 
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Fig. 6. Method for adjusting the zenith- and azimuth 
angles of the gun. 


TFerrella 


Fluorescent screen 


Fig. 7. Determination of the direction of the beam. 


Fig. 8. The angles of WE and D}, in the system of 
coordinates of the terrella. 


@° =NT.3 * Op (9) 


When the initial direction of the electrons 
lies in the E—W plane at latitude À = 0°, 
the error means that the direction is more 
easterly than is given from the geometric 
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direction of the gun. The correction amounted 
to 10°—15° in the lowest energy range. 

As a control the curves have been compared 
with values calculated by STÖRMER (cf. Astro- 
physica Norvegica, Vol. II No.1). These orbits 
are given in the following table together 
with corresponding values corrected for 
deflection inside the gun. 


Stormer Zenith (Stôrmer) 
orbit No. Nie angle 0 Ve 
70 223 BLO we SOs TI 72° A= latitude 
10 © ° ° A+acA 
68 3.03 * 10 38 24 24 
[e] (e} fe} 
2 : = 40; 46 Fig. 9. Sketch showing the calculation of the final 
©. 6 o 8 : S : 
Le ‘ 10 a 5 ° I 7 direction of the orbits. 
71 2.23 * 10 109 107 


The difference between the measured values 
Wg and the very accurate calculations by Stör- 
mer does not exceed 3° which is satisfactory. 

In order to compensate the deflection in- 
side the electron gun the axis of the gun has 
been tilted so much that the electrons leave 
the gun in the desired direction. 

When it was impossible to adjust the zenith 
angle because of the construction of the gun’s 
mounting (only angles of 8°, 16°, 24°... 
were possible) two zenith angles were measured 
and the intermediate value interpolated. The 
values thus obtained correspond to the direc- 
tions (N, S, E or W) in the diagrams. 


Calculation of Oy and VY; 
The angles D and W% are obtained from the 


measured distances m, n and | on the screen 
and the angles v and A (fig. 9) using the 


following trigonometric relations: 


Sot = Î sin v +n cos v 
Som= (I cos v—n sin v) sin A—m cos A 


Son = — (I cos y—nsin v) cos A—msin A (10) 
Son Sox 
tg On = == Wer 


V Sor? + Som? 


The first 3 formulae show the transition 
from the system of coordinates at the screen 
to the system (OL, OM, ON) at the terrella. 

Over the whole range of energies 10,000 
points of measurements were required. In 
these circumstances even the calculations by 
equations (10) would take a long time and a 
simple calculating machine was built to reduce 


the work. The time for calculation of Som and 
Son, namely about 6 minutes manually, was 
reduced to about 1.5 minutes. The measure- 
ment and complete calculation required about 
6 minutes for each point. | 

Fig. 10 shows the principle used for the 
calculating machine. The accuracy of the 
machine was approximately 0.5 %. 


R-R, /0W 


Fig. 10. The principle of the calculating machine. The 
following expressions can be solved: 


Som = (cos v—n-sin v)-sin A—m cos A 
Son = (lcos v—n-sin v)-cos A—m sin A 
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Measurements 


Measurements have been taken for each 10° 
of latitude beginning at the equator. ®y and 
Wr have been drawn for each latitude and for 
paths starting in North, South, East and West 
directions as a function of momentum. Each 
latitude is represented by four diagrams: 


VE in the N—S plane 
On in the N—S » 
PE in the E—W » 
On in the E-W » 


In the diagram the dot-dashed curve marked 
9° corresponds to particles which arrive from 
zenith; curves marked 8° N, 16° N, 24°N, etc. 
correspond to particles arriving from north 
and make angles of 8°, 16°, 24°, etc. with the 
zenith. Curves marked 8° S, 16° S, 24° S, etc. 
correspond to particles arriving in the same 
way from South. 

The diagrams give the direction of the 
path at an infinite distance from the dipole 
as a function of momentum. The particles 
arriving at the earth are positively charged. 


Accuracy of measurements 


The results of the experiment are to be 
applied to research on cosmic rays in the 
magnetic field of the earth. The obtained 
average error in the experiment did not exceed 
2—3°. Experimental errors of this magnitude, 
however, can be neglected, because of the 
considerable deviation of the terrestrial mag- 
netic field from a dipole field. : 

Accuracy in determination of the angles VE 
and ®y depends mostly on errors in reading 
of m, n and v. According to equation (10) we 
have 


(5D)? + cos? DB - (öy)? = 


= (y (7) + on 


If/ = 100 mm, m = 2 mm, # = 2 mm, v = 2 


(11) 


O 


WE get 


V (OP)? + cos? D - (0P)?= 3°.5 


in the greater majority of cases. This gives an 
= fe} 
average error in ® and Y of 2—3°. 
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It should be pointed out that, according to 
equation (11), the uncertainty rises rapidly 
when D 90°. This is not a real error but 
depends only on the chosen system of coordi- 
nates. In these special cases the velocity vector 
is projected on a plane nearly perpendicular 
to itself; therefore, it is difficult to determine 
Wr. An approximate calculation shows, how- 
ever, that even when @y reaches values of 
7o—75°, the error still does not exceed 3°. 
Curves where the error is expected to exceed 
3° are dashed. As can be seen from the dia- 
grams, this happens when ®X is near 90°. 

For accuracy in determination of momen- 
tum p we have the following expression 


Ap = DA Vie h „. À Rm A Mn 
p 2 Va | |x I M, 
derived from equation (6). 
A Vin AM, A Re / 
If 7, = 1% AL a R = 
we obtain + 5.3, % 


The investigation of systematic errors was 
carried out by comparing measurements at 
one place in the northern hemisphere, it’s 
antipodal point in the southern hemisphere 
and two more points symmetrical to the 
equatorial plane. Deviation of these values 
was supposed to show eventual asymmetry of 
the apparatus, influence of extraneous magnetic 
fields, etc. These deviations were also < 3°. 

Besides these investigations, a control of 
the accuracy was made by comparison of the 
results with values calculated by STÖRMER (cf. 
Astrophysica Norvegica, Vol. II No. 1). The 
accuracy of his calculations was stated to be 


n n 
10 10 
5 5 
0 wen AV. 0 5 10 AB, 


Fig. 11. Random distribution of PE and D 
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about 0.2°. About 30 values were chosen 
from this work and special measurements were 
made with the apparatus. The x-axis in the 
diagrams of fig. 11 represents the angle error 
in degrees, and the y-axis the number of 
observations. The diagrams show that the 
average angle error is 2—2.5°, in rather good 
agreement with the expected range of error. 

In addition to this some Störmer orbits 
have been compared with values found by 
interpolation in the diagrams (1—38). 

The result is shown below. In most cases it 
was necessary to interpolate these values from 
several curves and diagrams; therefore, it is 
natural if the error is somewhat greater than 
found above. 


: VE | D 
a naar ue @ 
| tum angle E N 
mer) | mer) 
9| 1.59 > 10!|52°| 52° 1°| o |— 2°|— 71° 
I ” 76°| 79° |-19°| o |—22°| 2° 
63 fs © || 332 |—21°|| 0 |—2021 0 
19 | 6.21 - 1010|66° 63° 38° O 30° 5° 
64 » D || CO ari | Oy 143 0 
19 | 5.01 - 10!9165°| 70° 70°| o CANNES 
65 Y o | 73° |—27°| o |—27°| o 
66 | 3.96 1029] o | 70° |-13° | o |— 14° 1° 
31 u Sun 390110301 © 32-4 
45.12.352.20”157°| 40% 323160 314142 
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The threshold values of the forbidden cones 
could not be exactly determind because of 
the small dimensions of the apparatus. How- 
ever, the curves indicate a tendency toward 
the values calculated by LEMAITRE and VALLARTA 
for 0°, 20° and 30° latitude (cf. Phys. Rev. 
1936, Vol. 50., p. 500). 
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Abstract 


A method of calculating moving standard deviations to give an idea of long-periodic fluc- 
tuations in the interannual variability of summer temperature is described and discussed from a 
statistical point of view. The method was applied to summer temperature records from Is 
stations in Sweden and changes in the variability at those stations are treated. In addition the 
actual variability during the period 1920—39 is determined and compared to the interdiurnal 
variability of the same period. 

In the second part of this investigation are studied frequency fluctuations of different circula- 
tion patterns over southern Scandinavia based upon the catalogue of daily weather situations 
during the period 1881—1950 published by the German Weather Service of the U.S. Zone in 
1952. These fluctuations are compared with changes in the summer-temperature and its varia- 
bility in South Sweden during the same period. It is found that the recent rise of summer tem- 
perature has been connected with a definite decrease in frequency of north-meridional situations 
over Europe and with an increase of south and warm-meridional cases as well as zonal situations. 
High variability seems to be connected with meridional circulation, low variability with zonal 
circulation. In a special chapter the fluctuations in variability of the Stockholm summer tem- 


perature since 1756 are discussed. 


I. Introduction 


Investigations of changes in the temperature 
climate have always played a dominant role in 
studies of climatic fluctuations. Such investiga- 
tions usually deal with variations of mean 
temperature or of frequencies of temperature. 
It is generally taken for granted that fluctua- 
tions in climatic elements, as temperature, 
reflect in some way or another variations in 
the general circulation of the atmosphere. It 
is impossible, however, to conclude only from 
the temperature trend how the circulation has 
changed from one period to another. 

From glacial-meteorological investigations 
in northern Sweden, I have reached the con- 
clusion that not only the summer temperature 
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itself might be of importance in glacier 
retreat but also the variability from one sum- 
mer to another (WALLEN 1948). It also seemed 
to me that fluctuations in the variability of the 
temperature climate would probably even 
better reflect long-periodic fluctuations in the 
general circulation than do the long-period 
changes in temperature itself. In order to ia- 
vestigate this, J have taken up a study of the 
variability of the summer temperature in 
Sweden, which would in some respects 
complete the careful analysis of the fluctuatic ns 
in temperature itself carried out by Liy3- 
QUIST (1949). 

Variability in the temperature climate is 
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usually studied as the interdiurnal variability. 
Studies of long-periodic fluctuations in inter- 
diurnal variability would, however, require 
tremendous work lying beyond the scope of 
this investigation. Therefore, it has been nec- 
essary to adopt some other, less elaborate 
method. Some studies of the actual inter- 
diurnal variability has been made, however, 
and will be treated in a subsequent chapter. 
If the general circulation of the atmosphere 
is subject to quasi-periodic fluctuations of 
several years—which we certainly have reason 
to believe from earlier studies—it seems 
likely that the interannual variability of the 
temperature would very often change with 
these periods. For instance, if we assume that 
the general circulation oscillates between a 
zonal and more meridional type, as suggested 
by Witerr (1944), it may be accepted as a 
hypothesis that in Sweden the zonal type of 
circulation would imply a maritime climate 
of low variability and the meridional type a 
more continental climate of a larger variability. 


2. Statistical methods 


As the interannual variations in temperature 
are often rather great it is necessary even in 
studies of the temperature changes themselves 
to adopt some method of smoothing in order 
to separate more long-periodic variations. The 
usual method of smoothing applied in recent 
investigations has been that of overlapping 
means. 

It occurred that it might be possible to 
adopt some method of overlapping means 
also to the problem of second order: the fluc- 
tuations of the variability. In all series of 
temperature observations it is possible to 
calculate the p-year overlapping means and to 
determine, for each one of these overlapping 
periods, the standard deviation, giving a 
measure of the interannual variability of that 
same period. The graphical representation of 
these standard deviations would then give 
an idea of the long-range fluctuations of this 
variability. We may call the temperatures in a 
series of n summers (,, fo, fg... . fn, where I, 2 

n represent consecutive years. The p- 
years overlapping means (p< n) are 


p pet p+2 
at; Zr at; 
be = igh ae a ee 
1 » “Mo 7 Ma bise se 
p P 
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where 1, 2... indicate the order of the over- 
lapping periods, which are initiated by the 


temp. f1, fo etc. ' 
The standard deviations (c) for a single 
year of these periods are then 


p | p+I 3 
»2 (t; SE Em) > (t; tng) 
I 2 
a= \/ — no a 
2 p ; p 


There is no doubt but that this method of 
moving standard deviations gives some idea 
of the long-periodic fluctuations in the inter- 
annual variability of the climate, but there are 
the following complications involved with 
the method. 

We have presumed, and it has also evi- 
dently been proven by the method of over- 
lapping means, that the summer temperature 
is in itself subject to long-range fluctuations 
which must influence the year to year varia- 
bility and also the fluctuations demonstrated 
by the moving standard deviations. An in- 
creased variability may be caused both by a 
rise in the year to year temperature range 
and also to some extent by a more or less 
sudden transition from a period of low to 
one of high temperature or vice versa. In 
other words, the general trend of the original 
temperature curve influences the variability. 

In fig. ı the broken curve represents a part 
of an original temperature curve, the dotted - 


Tempt 


2 + 6 4 /0 la 4 16 72 20 22 & 26 2. 

Years 
Fig. 1. Principal sketch of a temperature curve with 
secular trend. Broken line: original curve; dotted line: 
overlapping mean curve; heavy straight line: probable 
trend of the original curve; light horizontal lines: mean 
for certain periods around which the standard devia- 

tion is calculated. 
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curve the overlapping means and the heavy 
straight line the probable trend of the original 
curve. The light horizontal lines represent the 
mean around which the standard deviation 
is calculated for some periods. If a trend 
does exist, the right method of studying the 
variability is, according to more exact sta- 
tistics, to calculate overlapping means of the 
standard deviations from the trend line and 
not from the horizontal lines representing the 
mean value as suggested in the above the- 
ory. We can express this mathematically as 
follows. 

Let o be the value of the standard deviation 
obtained by the method suggested above. 
Then as before mentioned 


Let o’ be the value of the standard deviation 
from the trend line T. This value is expressed by 


> (t — T)? 
ar= +——___- (2) 


p 


The relationship between the sums of the 
squares in these two expressions is quite 


generally 
I (i —tm)?= DL (te —T)? + kr (m)? (3) 


Where k expresses the slope of the trend line, 
T; are consecutive time values and r,, is the 
average of the time values. The last sum of 
squares can be rewritten 


Zen) 2-2) EE ( 


Inserting (1), (2) and (4) in (3) gives 


ee (5) 


12 


of? =o? kt, 


This expression gives an idea of the in- 
fluence upon o of a linear trend in the original 
temperature curve. Unfortunately it is not 
possible to determine in any exact way 
the trend of the original temperature curve. 
Thus it does not seem reasonable to try to 
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eliminate the exaggeration of the variability, 
introduced by long-range fluctuations of the 
temperature, as we do not even know if the 
trend always can be expressed as a linear 
one. However, as can be seen from fig. 1, 
the overlapping mean curve is a very good 
approximation to the real trend curve. When 
obvious trends do occur and can be estimated 
from the overlapping mean curve, the above 
formula can be used to estimate the influence 
of the trend upon the value of o. 

Now it is easily shown that the influence 
of linear trends of the order of magnitude 
occurring in the temperature observations 
used, is often very little. In order to do so it 
must first be decided which value of p should 
be used for the overlapping means. 

According to LILJEQUIST (1949) it is re- 
commended to use a period of years about 
half as long as the period of the climatic 
fluctuations to be studied. As we are interested 
in fluctuations greater than 5 to 10 years we 
should probably apply 5-year overlapping 
means. This is, however, not possible because 
it would give us to few years for a deter- 
mination of the standard deviations. If these 
deviations should have any meaning they 
must be calculated from a series of numbers 
including at least 20 values. Consequently it is 
necessary with this method to use at least 
20-years overlapping means. This implies 
then that we can not be sure of the phase 
fluctuations shorter than 20 years. 

If we adopt a value of 20 for p, equ.(s) 
becomes 


OÙ = 6733.3 k? (6) 


From fig. 1 and curves shown later on we 
find that the largest probable value of the 
temperature change during a 20-years period 
is around 2°C which implies a continuous 
increase or a decrease of the summer tempera- 
ture by o°.1/year. A reasonable mean value 
of o is 1°.00 C. Accepting these values, we 
obtain from equation(6) o’ = 0.82, implying 
that the variability obtained by the suggested 
method of moving standard deviations is in 
this case 15—20 %. When a high interannual 
variability is connected with a linear trend, 
the exaggeration is smaller. If we assume k = 
0.1, and o = I°.50, we get o’ = 1.30. 

The exaggeration of the approximate 
method is then quite small. However, if the 
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Fig. 2. Relation between o’ and o for some different 
values of k. 


year to year variability is small compared to 
the trend, the exaggeration is of larger signi- 
ficance. Assuming o = 0°.7 (a minimum 
value) and k = o°.ı/year we get o’ = 0.40. 

Fig. 2 shows the relation between o’ and o 
for some values of k. The rapid increase of the 
exaggeration with a decrease of o implies 
that the top values of fluctuations in the 
moving deviations are rather independent of 
the trend while the bottom values become too 
high. The fluctuations of the variability are 
thus always underestimated and not overesti- 
mated with the method suggested. 

In fig. 2 is also drawn a curve for k = 0.05 
showing that such a weak linear trend has no 
influence of importance if ¢ > 0.5. It will be 
seen later on that no values < 0.5 have been 
found. 

The above discussion has shown then that 
the first mentioned method of moving 
standard deviations can be applied if due 
attention is paid to the influence of trends in 
the original temperature curve. Therefore, 
in the following analysis this method has been 
used; but, special discussion has been taken 
up in every case where the overlapping mean 
curves show a trend of importance to be 
hidden in the original temperature curve. 


3. Variability of summer temperature in 
different parts of Sweden 
The actual variability 


20-years overlapping standard deviations of 
the summer temperature (June—Aug.) have 


Karesuando Tarnaby 
Riksgränsen Gäddede 
Haparanda Ostersund 
Jokkmokk Härnôsand 
Falun Göteborg 
Skara Växjö 
Linköping Lund 
Visby Stockholm 


been calculated for the following stations in 
Sweden in order to study changes from one 
place to another in the variability (tab. 1): 

The period of observation from these 
stations is rather short, starting for most of 
them around 1860 and for two of them? as 
late as around 1900. The observation period 
from Stockholm is considerably longer, start- 
ing in 1758. It will be dealt with in a special 
chapter below. All stations have reliable and 
homogeneous series of observations. 

Before entering into the question of long- 
range fluctuations of the variability, it seems 
advisable to discuss the actual variability of the 
summer temperature in different parts of 
Sweden. 

For that purpose we have selected the 20- 
year period of 1920—1939, which gives an idea 
of the order of magnitude of the variability 
during a recent period of high variability. The 
values from this period should, in fact, all be 
somewhat reduced, for a general trend of tempe- 
rature increase occurs during the period. How- 
ever, as seen from the Stockholm curve fig. 10 
for example, this reduction will be only around 
10 % as all standard deviation values reach 
around 1.30 and the trend slope never exceeds 
o.°1/year. The value is in this way changed 
from 1.32 to 1.20. In fig. 3 maps are given of 
the mean temperature and its standard devia- 
tion in the period of 1920—39 at the above- 
mentioned stations (cf. also tab. 1). In order 
to clarify even more the distribution of mean 
temperature and its standard deviation such 
values have been calculated also for the 
following stations: Riksgränsen, Umeä, Sten- 
sele, Sveg and Karlstad. The values have not 
been reduced for the trend in the temperature 
curve. 

The decrease of the summer temperature 
with increasing latitude and from east to west 
is very typical in northern Sweden. It is also 
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Fig. 3a. Geographical distribution of summer mean 
temperature (June—Aug.) in Sweden during the period 
1920—39. 


evident that there are only small variations in 
Svealand and Gôtaland, but nevertheless two 
maxima can be distinguished: one in central 
Svealand around the lakes of Mälaren and 
Hjälmaren and another on the west coast. 
On the whole the standard deviation in 
different parts of the country varies oppositely 
of the mean temperature. The absolute maxima 
of variability are found at the high mountain 
stations of Riksgränsen and Storlien, which 
show the lowest mean temperature. The 
lowest variability is found in the southern- 
most part of the country where the mean 
temperature is considerably higher. On the 
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Fig. 3b. 
deviation of summer temperature during the same 


period. 


Geographical 


other hand there is a variability maximum 
connected with the temperature maximum 
in the region around the middle Swedish 
lakes. 

There is nothing surprising in the relatively 
high values of variability in the mountain 
region. These values reflect the wellknown 
fact that the Scandinavian mountains are 
situated on a borderline zone between two 
different climatic areas: the Atlantic to the 
west is the summertime source of cool, 
maritime air-masses, and eastern Europe is 
the source of warm, continental air. The 
high variability shows the sensitiveness of the 
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Table 1. 20 years overlopping means and standard deviations 
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of summer temperature at selected Swedish stations. 
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1859—78 
1860—79 
1861—80 
1862—81 
1863—82 
1864—83 
1865—84 
1866—85 
1867—86 
1868—87 
1869—88 
1870—89 
187I—90 
1872—9I 
1873—92 
1874—93 
1875—94 
1876—95 
1877—96 
1878—97 
1879—98 
1880—99 
1881—oo 
1882—o1 
1833—02 
1884—03 
1885—04 
1886—05 
1887—06 
1888—07 
1889—08 
1890—09 
189I—10 
1892—11 
1893—12 
1894—13 
1895—14 
1896—15 
1897—16 
1898—17 
1899—18 
1900—19 
I9Q0I—20 
1902—21 
1903—22 
1904—23 
1905—24 
1906—25 
1907—26 
1908—27 
1909— 28 
1910—29 
I9II—30 
(cont.) 
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mountains for an interchange between mari- 
time and continental summers. 

It is also evident from the map that regions in 
the mountains of generally higher maritimity 
show in summer a greater variability than do 
more continental regions!. For instance, the 
pass across the Norwegian border from Stor- 
lien to Ostersund is in summer clearly de- 
monstrated as a region of higher temperature 
variability compared to the more shielded 
areas in Västerbotten. 

The general decrease of variability as one 
gets further south seems to be related to the 
shifting northwards in summer of the sub- 
tropical high-pressure belt. No doubt, the 
chances for more constantly warm summers 
must increase in areas more or less regularly 
under control of the subtropical high which 
extends over middle Europe in summer. 

In south Sweden, finally, it seems as if the 
west coast is subject to a higher variability 
than is the east coast for we find a value of 
1.18 in Gothenburg but only 1.05 in Visby. 
If we consider the fact that the trend towards a 
secular increase of temperature is larger in 
Visby than in Gothenburg, the standard 


1 This increase of variability with increasing ma- 
ritimity is rather unexpected and may be caused by 
influences of the topography. Sumnet (1953) has shown 
the variability to increase in North-America with the 
continentality. 


deviation from Visby is too high and the differ- 
ence between the east and the west coast 
even greater. This fact suggests the influence 
of the zonal westerly winds to be of special 
importance in the variability on the west coast. 

Summarizing, the variability of the summer 
temperature in different parts of Sweden can 
be expressed by the following means: 


Gôtaland = 1°.09 C Norrland coast and for- 
(6 stations) est land = 132€ 

(7 stations) 

Norrland mountain re- 
gion =r —saiG 

(5 stations) 


Svealand 7.2310 
(3 stations) 


Fluctuations of the variability 


It is already evident from a general inspection 
of the curves in fig. 5 that in the recent cen- 
tury the period of 1920—39 is one of relative- 
ly high variability. It is, therefore, interesting 
to compare the distribution of the variability 
of this recent period with an older one of 
low variability. In fig. 4 the temperature and 
its variability are plotted for the period of 
1880—99 at the same stations as for the recent 
period and isolines are drawn as in Ta 

It is seen that the distribution of variability 
is somewhat different from the one during the 
more recent period. A general decrease of 
variability from north to south is certainly 
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found even during 1880—99, but the region 
of low variability in south eastern Norrland 
is much more accentuated than in recent 
time. The smallest changes from the older to 
the more recent period, reaching not more 
than 20 %, are found in the far north and in 
the south. In southern Norrland and central 
Sweden, especially in south eastern Norrland, 
the change is considerable (from 1.20 to 0.80 
or 30 %); but, generally speaking, the mari- 
time regions seem to show the largest differ- 
ences. They are evidently more sensitive to 
changes of the variability than are the con- 
tinental areas. 

It should be emphasized that the general 
temperature distribution over the country 
does not show any significant changes from the 
period 1880—99 to 1920—39. This is seen from 
a comparison between fig. 3 a and fig. 4 a. 
This fact might suggest that the change of the 
variability is of larger significance in studies 
of the variations of the general circulation 
than is the mean temperature itself. An attempt 
to explain the changes in the variability with 
regard to fluctuations in the general circulation 
is made in chapter 5. 

From the comparison between the varia- 
bility of the two periods it is then clear that 
even inside a comparatively small land area 
as Sweden considerable differences in the 
fluctuations of variability do occur between 


Tellus V (1953), 2 


various parts of the region. This fact becomes 
even more evident from an inspection of the 
moving standard deviation curves for stations 
in different parts of the country (fig. Ge 

In northern Norrland the variability shows a 
maximum during a period from 1890 to 1910 
while there are very small indications of a 
maximum at that time in southern Norrland. 
In middle Sweden this maximum again shows 
up; it is, however, of much less intensity than 
in the northernmost part of the country. Only 
in Visby is it again very marked. Practically 
no maximum is found in southern Sweden 
at that time. 

We find a very marked minimum in the 
variability around 1915. This seems to hold 
good for all the country but the intensity is 
much larger in the northern than in the south- 
ern parts. This variability minimum is most 
accentuated in Karesuando and Haparanda 
and practically insignificant in Géteborg and 
Lund. 

During the following decades an increase 
connected with the rise of the summer mean 
temperature is rather clear over most of the 
country. This increase is—as was to be ex- 
pected—most accentuated in northern Sweden; 
at some stations, however, this was caused by 
the somewhat larger rise of temperature in 
this part of the country. The very small rise 
of the variability shown in the curves from the 


Fig. 4a. Geographical distribution of summer mean 
temperature in Sweden during the period 1880—99. 


southern parts—at some places only insig- 
nificant secondary maximas are found—is 
rather surprising. Although smaller, there is a 
very definite rise of the mean temperature 
also in this region. 

Finally, a decrease of the variability has 
occurred in all parts of the country durin 
the last decade showing up in a fall of the 
curves since around 1930. 

Looking at the country as a whole we can 
summarize as follows. The period around 
1880 was in all parts a period of low variability 
being then everywhere relieved by a very 
definite increase culminating in northern Norr- 
land around 1895 and showing in southern 
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Fig. 4b. Geographical distribution of the standard de- 
viation of summer temperature during the same perio d 


Sweden a maximum somewhat later, about 
1905. In southern Norrland and middle Sweden 
there is no culmination around the turn of 
the century but a continuous rise until the 
maximum of the period centered around 1930. 
In Norrland the maximum at the turn of the 
century was followed by a minimum around 
1915 showing up rather smoothed out even 
in the curves for southern Sweden but to a 
very small extent in southern Norrland. The 
rise from this minimum to the maximum of 
1930 is very much accentuated in northern 
and middle Sweden but is of no significance 
in south Sweden and in Gotland. 

From this summary it seems possible to 
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Fig. 5 a. 20-years overlapping means of summer temperature (right) and its standard deviation (left) 
at selected stations. 


conclude that the country can be separated 
into three regions showing definite differences 
in the change of the variability; namely, 
northern Norrland, southern Norrland and 
south Sweden. There are, however, certain 
reasons to conclude that differences exist 
even inside south Sweden. At least we might 
say that the curve of Visby definitely shows 
different features than do the others. Further- 
more, the curves of Gothenburg and Lund do 
not show any typical minimum about 1915 
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as do the other curves and likewise show no 
typical maximum at the turn of the century. 
There seems to exist a maritime western type 
and a maritime eastern type; this is by no 
means surprising if we consider the differ- 
ences in the influence of the normal general 
circulation in summer on the Swedish west 
coast and in Gotland. From this reasoning the 
decision was made to separate the stations in 
the following groups representing different 
types of fluctuations in the variability: 
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Fig. 5 b. 20-years overlapping means of summer temperature (right) and its standard deviation (left) 
at selected stations. 


1. Northern Norrland 2. Southern Norrland and northern Svealand 
Karesuando (Gaddede) 
Jokkmokk Ostersund 
(Tarnaby) Härnösand 
Haparanda Falun 
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3. Central middle and south Sweden 


Stockholm 
Linköping 
Skara 
Växjö 


4. Gotland 
Visby 


5. The west coast 


Göteborg 
Lund 


4. The interdiurnal summer temperature 
variability in Sweden 


In order to give a general idea of the order 
of magnitude of the interdiurnal summer 
variability in Sweden, the values given in 
table 2 have also been calculated. 

From the table it is evident that the inter- 
diurnal variability decreases during the sum- 
mer being everywhere largest in June and 
smallest in August. It is also clear that in general 
it decreases rapidly from north to south being 
in all months much larger in Gällivare than 
at the other stations. The order of magnitude 
of the interdiurnal variability is surprisingly 
similar to that of the interannual variability 
and although the calculations have been made 
for too few stations to give any definite re- 
sults there seems to be indications that the 
variation of the interdiurnal variability inside 
Sweden is rather similar to the interannual 
summer variability as given in fig. 3 b. This 
is indicated from a comparison between the 
values of col. 4 and col. 5 of table 2. 


Table 2. Interdiurnal variability of summer tem- 
perature 1920—1939 


Inter- 


June July Aug. er zu 

Summer 
allvare..... 1.90. 1.70 1.40 1.67, 1.34 
Härnösand.. 1.50 1.20 1.10 1.27 1.26 
Stockholm.. 1.60 1.40 1.00 1.33 1.32 
Göteborg... 1.20 1.10 0.90 1.07 1.18 
Malmö..... £4051 10 MNOO 1.132 1.00 
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5. Changes in summer circulation pattern 
over Scandinavia since 1880 and their rela- 
tion to temperature and its variability 


Changes in the circulation pattern. 


In order to investigate to what extent the 
fluctuations in temperature variability re- 
flect changes in the general circulation it is 
necessary to make in some way an analysis of 
such changes that have occurred over Europe 
during the last decades. There has been much 
discussion on the increased general circulation 
in middle latitudes during recent years but 
very few investigations have been made in 
order to find out in which way the general 
circulation patterns have changed. 

WiILLETT (1944) has introduced, in his 
extensive studies on the variations in the gene- 
ral circulation in all middle latitudes, the im- 
portant definitions: zonal and meridional, 
which characterize two principally different 
types of circulation: high-index and low- 
index circulation resp. These definitions have 
been used to a considerable extent by authors 
dealing with circulation patterns. Among these 
I wish to mention PETTERSSEN (1949), who 
has investigated the basic circulation changes 
of this century over the North Atlantic- 
European sector of the hemisphere. This in- 
vestigation, however, only deals with the 
winter. Nothing very definite has been 
published concerning the summer circulation. 
However, some indications regarding an in- 
creased importance of “blocking-situations” 
during the summers of the last decades, im- 
plying a special type of low-index, were 
found in a study of synoptic climatology by 
REX (1950, 1951). 

In 1952 the German weather service of 
the U.S. Zone published (Hess and Bre- 
ZOVSKY, 1952) an extensive catalogue of daily 
weather situations having occurred over 
Europe since 1881, principally based upon a 
classification of circulation patterns introduced 
by Baur. Of course this classification was made 
with particular reference to middle Europe 
and therefore is not especially well-suited for 
an analysis of secular variations in the circula- 
tion pattern over Scandinavia. Such an ana- 
lysis has been made, however, since it was 
expected that it would, in any event, give 
some indications as to how the circulation 
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patterns have changed at least in south and 
middle Sweden. 

Hess and Brezovsky have made their 
classification from soo mb upper-air charts, 
and introduced Willetts basic idea of zonal 
and meridional circulation types. Among 
eighteen different classes of weather situations, 
two are considered purely zonal in character, 
three are mixed situations and the rest are 
considered meridional. The authors have in 
every situation distinguished between cyclonic 
and anticyclonic situations and so obtained 
at least two sub-classes for each principal 
class. In characterizing the weather situation 
of each day since 1881, most of these principal 
and sub-classes were used. It is certainly im- 
possible to apply in Scandinavia all the charac- 
teristics of the different classes expressed by 
these authors; for instance, a situation being 
considered as anticyclonic in central Europe 
may very well be cyclonic over Scandinavia. 
On the whole, however, the general principle 
of zonal and meridional situations hold good 
also over central and southern Scandinavia. 
Among the mixed situations there are some 
which might be referred to as zonal in Scan- 
dinavia but as meridional over central Europe. 

In beginning, an attempt has been made to 
determine, from among the basic situations, 
which are zonal and which are meridional 
over Scandinavia. It was accepted that the 
following types were zonal (the designators 
and characteristics are the ones used by Hess 
and Brezovsky and the air-flow everywhere 
refers to southern Scandinavia). 


NW Cyclonic or anti-cyclonic north-westerly 
zonal flow. 
W Cyclonic or anti-cyclonic westerly flow. 


SW Cyclonic or anti-cyclonic south-westerly 
flow, considered as a mixed type by the 
authors. 

BM _ High-pressure ridge from the Azores to 


eastern Europe giving south-westerly to 
westerly zonal flow. 


The following situations were referred to 
as definitely meridional. 


HN _ High-pressure area over the North Sea; 
northerly flow. 

HF High-pressure area over Scandinavia; 
light variable flow. 

HNE High-pressure over the North Sea and 
Scandinavia; light variable flow. 


CG. VERLIEREN 


N Northerly meridional flow mostly caus- 
ed by blocking-highs over the eastern 
Atlantic. 

Trough over northern and middle 
Europe giving southerly meridional 
flow. 

TrW Trough over the North Sea and eastern 

Atlantic; southerly meridional flow. 


TM Cut-off low, central Europe; south- 
easterly meridional flow. 
TB Cut-off low, British Isles and eastern 


Atlantic; mainly southerly flow. 


S Blocking-ridge Scandinavia and the 
Baltic; southerly meridional flow. 

SE Blocking-ridge from south-east toward 
north-west; south-easterly flow in south- 
ern Scandinavia, generally north-westerly 
in the north. 


NE High-pressure ridge from British Isles 
over Scandinavia to Russia. Mainly 
southerly flow. 

Ww Blocking-high, central Russia; low 


eastern Atlantic. Southerly meridional 
flow. 


Besides there are two types of situations, 
which are, as type SE, difficult to classify 
definitely in the zonal or meridional groups: 


HM _ Expanded high-pressure over most of 
central Europe giving sometimes zonal, 
but more often north-westerly flow 
over Scandinavia. 

HB Blocking-high, British Isles giving gen- 
erally northerly meridional flow over 
Scandinavia, this flow transporting, 
however, more or less warm air masses 
around the high ridge extending north- 
wards. 


Both of these types have been referred to the 
meridional groups. 


The frequency of different situations during 
the summer months June—August was cal- 
culated for each year from 1881—1950 accord- 
ing to this classification. The frequency of 
all zonal and all meridional situations was then 
determined for each year, and finally 20-years 
overlapping means were calculated from these 
frequencies. The overlapping mean curves are 
shown in fig. 6 together with 20 years over- 
lapping means for Stockholm summer tem- 
perature from the same period. 
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Table 3. 20 years overlapping means of frequencies of certain summer circulation patterns 


Period | Zon | Mer | S | N w |s+w| HM | HB | Z+N | 
er eee EEE eee TE IE ee eee 
1881—00 48.4 43.0 7.8 1222 9.4 17.2 10.2 4.1 60.6 
1882—o1 47.3 44.7 7.9 122 10.5 18.4 9.8 4.4 59.5 
1883—02 46.9 45.1 7.7 13.0 102 17.9 Io.I 4.2 59.9 
1884—03 46.8 45.2 7.9 13.0 10.2 18.1 10.1 4.0 59.8 
1885—04 47.3 44.7 75 13.9 8.7 16.2 10.4 4.3 61.2 
1886—05 47.1 44.9 7.7 13.2 9.2 16.8 10.9 2 60.3 
1887—06 46.4 45.0 7.9 13.4 9.5 17.4 10.6 4.3 59.8 
1888—07 47.5 44.5 8.9 1277 9.6 18.5 9.5 3.9 60.2 
1889—08 47.1 44.9 9.3 12.6 9.9 19.2 9.3 3-9 59.7 
1890—09 47.2 44.8 9.9 13.0 9.0 18.9 OL 3.9 60.2 
1891—10 45.8 46.2 10.2 13.8 2 19.4 9.1 4.0 59.6 
1892—II 44.8 47.2 10.1 14.3 Or 19.2 ‘9.4 4.4 59.1 
1893—12 44.2 47.8 11.3 13.7, 9.7 21.0 9.3 3.8 57.9 
1894—13 44.5 47.5 11.0 ne 10.3 21.3 9.2 3.4 58.2 
1895—14 2.8 49.2 12.0 nee 10.9 22.9 9.8 38 56.1 
1896—15 42.2 49.8 11.9 14.1 10.8 22.7 9.9 3.2 56.3 
1897—16 2.4 49.6 12.0 13.9 10.6 22.6 9.9 3.4 56.3 
1898—17 41.4 50.6 124 Tig) II.0 23.4 10.3 3.3 55.1 
1899—18 42.1 49.9 12.3 13.7 10.9 Bane 2 3.4 55.8 
1900—19 44.1 47.9 12.3 13.8 10.2 22.5 8.9 27 57.9 
1901—20 44.0 48.0 12.4 14.0 10.9 23 8.2 2.6 58.0 
1902—21 45.1 46.9 II.9 14.4 10.2 22.1 8.0 BAR 59.5 
1903—22 46.1 45.9 2 14.0 9.9 22.0 7-4 2.4 60.1 
1904— 23 46.5 45-5 11.5 14.2 9.4 20.9 7.8 2.4 60.7 
1905—24 47.2 44.8 102492 13.4 9.7 279 723 2-1 60.6 
1906—25 47.1 44.9 T2} 13.8 9.6 Ph els DP 60.9 
1907—26 46.9 45.1 12.4 13.6 9.3 27 Tee 2.5 60.5 
1908— 27 46.2 45.8 12.6 1327. 9.7 2223 7.6 223 59.9 
1909—28 A722 44.8 7212 111,2 9.1 PB 8.0 2.4 60.4 
IQlo—29 47.0 45.0 11.6 12.8 8.6 20.2 9.1 2.6 59.8 
IQII—30 47.1 44.9 Tae 11.8 9.4 20.9 9.7 2 58.9 
I9I2—3I 48.0 44.0 11.8 aii ey 8.5 20.3 9.7 2.5 59.7 
1913—32 48.6 43.4 10.6 12.4 7.9 18.5 9.9 2.6 61.0 
1914—33 48.7 43.3 10.6 11.9 7.5 18.1 10.5 DT 60.6 
1915—34 49.4 42.6 10.0 1225 7.5 17.5 10.0 2.8 61.7 
1916—35 49.6 2.4 10.5 ce 8.1 18.6 10.0 2.8 60.7 
1917—36 49.8 22 10.0 II.0 7.8 17.8 10.5 2 60.8 
LOLS 377 50.6 41.4 9.1 10.9 8.0 17.1 10.0 3.3 075 
1919—38 49.6 2.4 9.4 10.6 8.7 18.1 10.4 32 60.2 
1920—39 48.4 43 6 10.3 10.3 9.5 19.8 10.2 52 58.7 
1921—40 48.6 43.4 10.6 9.7 9.4 20.0 10.1 3.6 58.3 
1922—41 48.5 43.5 11.4 9.0 8.9 20.3 10.7 3.5 57.5 
1923—42 46.7 45.3 11.8 8.8 8.9 2047 12.0 3.8 55.5 
1924—43 46.2 45.8 12.3 8.6 9.2 21.5 11.9 3.8 54.8 
1925—44 45-4 46.6 12.4 9.6 9.6 22.0 11.4 4.0 55.0 
1926—45 45.8 46.2 12.5 9.2 9.0 21.5 11.5 4.0 55.0 
1927—46 45.6 46.4 12.5 8.7 9.8 228 11,7 4.1 54.3 
1928—47 45.3 46.7 13.1 8.5 9.0 2241 12.4 4.1 53.8 
1929—48 44.9 47.1 13.5 DB 8.9 ZA 12.7 4.0 53.2 


Provided that the above classification can have occurred since 1880. In the mean, howe- 
be used to describe the circulation types also ver, the two types of circulation balance each 
over southern Scandinavia, it is clear that other surprisingly well (tab. 3). The mean fre- 
very obvious interchanges between periods quency of zonal summer situations is 47 or al- 
of dominating zonal and meridional circulation most exactly half of the total (92). In the mean 
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Fig. 6. 20-years overlapping mean curves of frequen- 


cies of zonal (zon) and meridional (mer) circulation 
patterns in summer together with 20-years overlapping 
means of Stockholm summer temperature. 


for the whole period 1881—1950, a zonal 
circulation of any kind is, in the summer, as 
common as a meridional one. 

Looking at the fluctuations we find that 
the summerly frequency of meridional situa- 
tions increased very markedly from 1890 
until the period centered around 1910, reaching 
a maximum of sı cases out of 92. From the 
period 1898—1919 it dropped to a minimum 
at the end of the 1920's showing in the 
mean a frequency of only 41 cases out of the 
total. From the period 1918—37 it has in- 
creased continuously again. Evidently the 
frequency of zonal cases shows the opposite 
variations, i.e. a minimum between 1905 and 
1910 and a maximum at the end of the 20’s. 


The connection of circulation changes with tem- 
perature 


Comparing now these fluctuations with 
those we have found in the temperature 
variability, it seems reasonable to consider 
only the stations in southern Sweden as the 
applied classification seems more likely to be 
valid in the southern parts than in the northern, 
Generally speaking the fluctuations in the 
variability from 1881—1950 as given in fig. 5 
seem to coincide very well with the changes 
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found in the frequency of meridional situa- 
tions. In details, however, the coincidence is 
not so good. A maximum in variability is 
reached between 1905 and 1910, then it 
decreases to a minimum which occurred, 
however, much earlier than that for the fre- 
quency of meridional situations, namely around 
191$ instead of around 1925. In Visby, how- 
ever, a secondary minimum in variability is 
found as late as 1927. An increase of the 
variability and of the meridional frequency 
is found everywhere after the minimum 
but the coincidence in time is not at all per- 
fect. 

The question is now if it is likely that a 
high temperature variability would correspond 
to a high frequency of meridional situations. 
Before going into this question we must 
investigate how the frequency fluctuations do 
correlate with the changes in the summer 
temperature itself. 

In order to do so we have studied the over- 
lapping means of the frequencies of some 
combined circulation types. It was assumed 
that it is possible to combine those meridional 
situations in the above scheme giving warm, 
southerly flow over Scandinavia on the one 
hand and the ones giving cold, northerly 
flow on the other. The following situations 
were referred to a general group of warm- 
meridional situations (w): 

HF, HNF and NE 

The following were classified as south- 
meridional (s) situations: 

Ww, SE, S, TrM, TrW, TM, TB 
and the following as north-meridional (n) 
NFHNeHB, HM: 

Assuming that in summer HM, in most 
cases, must give above normal temperatures 
(REx, 1951) in Scandinavia, we have then 
combined the HM frequency with the fre- 
quency of warm and south-meridional situa- 
tions and calculated the correlation coefficient 
between this combined frequency and the 
summer temperature in Stockholm based 
upon the period 1881—1950. The coefficient 
was found to ber = + 0.55 + 0.05 showing 
a rather significant correlation. 

The correlation coefficient between the 
same temperature and the frequency of zonal 
situations was found to be r = — 0.41 + 0.06 
showing that the zonal fow generally gives 
temperatures below normal in Stockholm. 
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Fig. 7. 20-years overlapping means of frequencies of south-meridional (S), warm-meridional (W), and north- 
meridional (N) situations together with the corresponding overlapping mean curve for the Stockholm summer 
temperature. Combined curves (Z + N) and (W + S) are also entered. 


Finally, the coefficient between a combination 
of zonal and north-meridional situations (ex- 
cept HM cases) was determined to be r = 
— 0.51 + 0.06. 

From these calculations of correlation coeffi- 
cients, it can be concluded that the summer 
temperature may be considered as being rather 
closely related to the frequency of these 
basic types of circulation. 

The fluctuations in frequency of south- 
meridional (s), warm-meridional (w) and 
north-meridional (n) situations since 1881 are 

‘shown by 20-year overlapping means in fig. 7 
a and b together with the corresponding 
temperature curve from Stockholm. Combined 
curves (w +s) and (z +n) are also entered in 
the diagrams. 

It is evident from these figures that there 
has been a continuous increase in number of 
south-meridional situations but on the whole a 
decrease of north-meridional cases. The warm- 
meridional situations do not show sucha definite 
trend but a certain increase is found there also. 
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Comparing now these curves with the tem- 
perature curves from Stockholm and other 
places in south Sweden we find quite an inti- 
mate relationship between the temperature 
and the frequency of warm and south-meri- 
dional situations. Looking at the curves from 
Göteborg and Växjö (fig. 5), the increase of 
the summer temperature from around 1890 
to the period centered in 1905—07 is quite 
obvious. The frequency of warm-meridional 
situations continuously increased during that 
period. At the same time, however, the 
north-meridional cases—though being fewer— 
are also increasing in frequency. Socks 
natural that the temperature rise can not be 
so large as when the warm-meridional situa- 
tions increase at the cost of north-meridional 
instead of the zonal cases as occurs towards 
the end of the investigated period. In fact 
no particular rise in the temperature curve 
from 1890—1907 is found, for instance, at 
Stockholm. 

The temperature practically everywhere in 
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south Sweden shows a decrease from the maxi- 
mum period between 1905 and 1910 and re- 
mains then quite low until it again starts to 
rise around 1925. From the period around 1907 
until that of 1927 the frequency of warm and 
south-meridional situations is continuously 
falling. The number of north-meridional 
cases is also decreasing but with a lower rate. 
The frequency of zonal cases rises considerably, 
and this period of low temperature may be 
characterized as one of typically zonal circula- 
tion. From the middle of the 1920’s an extra- 
ordinary increase of the summer temperature 
has set in all over northern and middle Europe. 
According to the circulation changes found 
in fig. 7 the south and warm-meridional 
situations again started to increase consider- 
ably, this time at the cost of both zonal and 
north-meridional circulation patterns. 

As mentioned the HM and HB cases could 
not be definitely referred either to the north- 
meridional or to the warm-meridional group. 
From the point of view of circulation patterns 
they are mainly north-meridional over Scan- 
dinavia but give quite often rise in summer to 
high temperatures. On the other hand a 
very small shift westwards of the axis of high 
pressure may cause a considerable cooling-off 
over Sweden and a complete change in the 
distribution of temperature and precipitation 
over Scandinavia. This fact has been shown 
by Rex (1951, 1952) in his study of blocking- 

ighs over northwestern Europe. The HM 
and HB cases in the German classification are 
namely quite often identical with the blocking- 


highs studied by Rex. 
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Rex in his investigation stressed, however, 
that the development of the temperature condi- 
tions in Scandinavia during the last decades 
is exactly that which could be expected if the 
frequency of westerly blocking-highs—centered 
not too for west—has increased in the same 
period of time. In order to clarify this point we 
have studied the frequency fluctuations also 
of the HB and HM cases. 

The average number of HM-cases per 
summer during the whole period 1881—1950 
is 10 and of HB cases between 3 and 4. These 
numbers are so small that it is impossible that 
any high correlation could be found between 
the summer temperature at a single station 
and the frequency of cases. One finds that the 
correlation between the summer temperature 
in Stockholm and the frequency of HM cases, 
for instance, is expressed by r = 0.14, which 
certainly is quite insignificant. Looking at the 
long-range fluctuations illustrated by 20-year 
overlapping means, the coincidence between 
the frequency trend and temperature on the 
whole in Sweden is, however, quite remark- 
able. The general trend for HM and HB is 
quite similar so we may combine the two 
cases to obtain the curve in fig. 8. The rela- 
tively high mean temperature around the 
turn of the century corresponds to a well- 
defined secondary maximum in the frequency 
of HM and HB; the cool summers period 
from around 1907 to 1920 corresponds to a 
very marked minimum in the same frequency; 
and the considerable warming up in the last 
two and a half decades corresponds to a very 
strong increase in the frequency of HM and 
HB. It should be noted that this increase in 
the summer HM-frequency has already been 
stressed by Brezovsky (1952) in a special 
study of the fluctuations in frequency of 
different situations. It seems unlikely that this 
coincidence should be just accidental, so we 
may at least conclude that the summer mean 
temperature in Sweden on the whole is fa- 


voured by a high frequency of blocking- 


highs over the British Isles extending over. 


central Europe. 

The increase of summer temperature in middle 
and northern Europe since around 1925 then 
seems to have been caused by a considerable in- 
crease of the frequency of south-meridional and 
warm-meridional situations and also of blocking- 
highs appearing rather far to the east. This fact 
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is in agreement with the finding of Eriksson 
(1943) that in July the southerly flow over 
Scandinavia and over the northern-most At- 
lantic has increased from the beginning of the 
century to the 1930's. Such a change towards 
a more warm-meridional circulation over 
Scandinavia has been shown to have occurr- 
ed also in winter both by Eriksson (1943) 
and PETTERSSEN (1949). It seems rather obvious 
that this is the fundamental change of the cir- 
culation over Europe having caused the consider- 
able warming up in northern latitudes which has 
so evidently shown up, for instance, in the re- 
treat of glaciers (AHLMANN 1948, WALLEN 1948, 
WILLETT 1950). 


Changes in circulation and the temperature 
variability. 


Comparing now more conclusively the 
fluctuations of the general circulation with 
changes in the variability of the summer tempe- 
rature, it is found from fig. 9 and a comparison 
between figs. 5 and 6—8 that there is on the 
whole a definite tendency to high variability 
during the period 1890 to 1910; this period 
was, especially in the latter part, characterized 
by a very high frequency of meridional situa- 
tions. The decrease in the variability around 
1910 coincides with a development towards a 
more zonal circulation, although the maximum 
of zonal situations occurs much later than does 
the minimum of the variability. It should be 
noted, however, that the increase of the 
variability, starting again around 1915, is to 
some extent caused by the general temperature 
trend towards a warmer summer climate, 
which is in turn a consequence of the decreasing 
frequency of north-meridional situations and 
increasing frequency of south-meridional and 
HM and HB situations as well. The high 
variability is then probably caused by a rapid 
interchange between zonal summers on the 
one side and summers with high frequency of 
south-meridional and HM, HB situations on 
the other. The maximum of the variability in 
the latter 1920's coincides with a high frequen- 
cy of south and warm-meridional circulation 
patterns and a rather high frequency of zonal 
situations. Later on, when the frequency of 
zonal cases goes down together with a con- 
tinued decrease of north-meridional cases and 
the south-meridional frequency remains high, 
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Fig. 9. Stockholm summer remperature and variability 

(V) curves (T881—1950) compared to the overlapping 

mean curve of frequencies of meridional (M) circula- 
tion patterns. 


the variability decreases again (compare also 
fig. 9). 

Summarizing the Auctuations in sum- 
mer circulation over north-western Europe 
since 1881 and comparing with fluctuations in 
temperature and its variability in southern 
Sweden we recognize the following periods: 


Around 1890 to around 1910: 

Period of summer temperature slightly 
above normal and high temperature varia- 
bility. The circulation was characterized 
by a high frequency of meridional circu- 
lation both from a northerly and southerly 
direction. The frequency of HM and HB 
cases was comparatively high. The zonal 
circulation was low and decreasing 
throughout the period. 


Around 1910 to around 1927: 

Period of summer temperature around 
or below normal, temperature and its 
variability increased during the latter part 
of the period. 

Circulation was first characterized by 
a rapid decrease in frequency of meridional 
situations and an increase of the zonal 
flow. South-meridional situations became 
more rare and HM and HB cases showed 
a definite minimum at the same time as 
the minimum in variability (around 1915). 

Towards the end of the period the 
frequency of south-meridional cases in- 
creased again. 


Around 1927 to 1950: 
Period of increasing summer-tempera- 
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ture and variability. The latter decreased 
at the end of the period. 

Circulation was characterized by a 
continuous increase of the frequency of 
south and warm-meridional situations and 
HM, HB cases. North-meridional and 
zonal circulation decreased rapidly. When, 
during the 40’s, the frequency of warm 
and south-meridional cases kept high 
together with that of HM and HB cases 
the variability decreased again. 


6. Analysis of the complete Stockholm 
temperature curve 


Temperature observations are available from 
Stockholm since 1756 so it is possible to 
make a particularly interesting study of the 
fluctuations of the variability. In fig. 1o is 
drawn the year to year changes in the summer 
temperature (June—Aug.) together with the 
20-years overlapping mean curve. Above 
these curves the 20-year overlapping standard 
deviation curve is given. 

From inspection it is even more evident 
than from the other curves that the long- 
range fluctuations in the variability really 
correspond to periods of a large year to year 
variability. A period of continuously rather 
high summer temperature in the 1750 and 60’s 
is followed by still warmer decades but of 
much higher variability. This suggests that the 
famous warm period at the end of the 18th 
century was characterized by a quite high 
variability (1.2—1.3). There is no doubt 
that the high variability is somewhat in- 
fluenced by the general increase of tempera- 
ture, but from the trend shown by the over- 
lapping mean curve, it is obvious that this 
influence can not be dominant. 

During the first decades of the roth century 
the summer temperature begins to fall and 
there occurs a slight decrease in the variability 
corresponding to a cooler and more monot- 
onous period of summers. Between 1810 and 
1825 there is another period of increasing 
variability combined with a much cooler 
climate than had the period at the end of the 
preceding century. After 1825 the summer 
temperature continues to drop while the 
variability shows a small decrease but still 
remains rather high. From the beginning of 
the 1840's a very typical increase in the mean 
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temperature is caused by a certain number 
of very warm summers interrupted by rather 
cold ones. This means that the variability 
again reaches very high values during the 
1850'S. 

A most remarkable drop in the summer 
temperature curve of Stockholm occurred at 
the beginning of the 1860's. For about Io years 
after 1860 the summers were extremely cold 
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and even later on remained cool for several dec- 
ades. Contemporarily the variability decreased 
quite rapidly around 1860 and remained con- 
tinuously low until the 1890’s. It should be 
noted that the influence of the temperature de- 
crease is to increase the variability so the 
values of variability obtained are in fact 
somewhat to high. Then again—as has been 
discussed above—the variability increased while 
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temperature was rather unchanged and remain- 
ed low until around 1925 when rather suddenly 
it started to increase very definitely. The varia- 
bility, having dropped again slightly around 
1915, started to rise very soon thereafter and 
reached a maximum for the period 1920— 
1939. Since then it has dropped again. 

Considering now the more general trend 
of the fluctuations in the Stockholm summer 
temperature and its variability, it is seen 
from fig. to that the summers generally 
became colder and colder from the middle of 
the 18th century until the 1860's and then, 
on the whole, remained cool until 1925. 
Since then they have become considerably 
warmer. During the same time the variability 
has changed as follows. 

During the last half of the 18th century the 
variability increased considerably and then 
remained high or rather increased (maxi- 
mum 1.60) for a century until it suddenly 
dropped to low values in connection with the 
deterioration of the summer climate around 
1860. During this century of high variability 
there were short periods of decrease, corre- 
lating fairly well with periods of cool sum- 
mers. On the whole, the variability during this 
first century rather definitely differed from 
that during the following one. After 1860 
the variability has been particularly low un- 
til it rose at the turn of the century. On 
the whole the variability during the century 
since 1860 was lower than during the foregoing 
one especially if we note that the rise after 
1925 was to some extent caused by the general 
temperature increase starting at that time. 

Even if the conclusions given in chapter s 
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regarding the relation between the variability 
and changes of the general circulation have 
to be proven further it is interesting to make 
an attempt to analyse the older part of the 
Stockholm temperature curve. It seems, at 
least, probable that the long period of high 
variability in the summer temperature, cha- 
racterizing the century from 1760 to 1860, 
was one of comparatively high frequency of 
meridional circulation. As the period from 
1860 to our days shows on the whole a lower 
variability it is likely to have been of a more 
zonal character than the first century. If it is 
allowed to draw conclusions for shorter periods 
we may suggest that the warm period 1770— 
1795 was one of a typical meridional cir- 
culation with dominating south to east merid- 
ional circulation. 1800—1820 shows a definite 
cooling off with lower variability indicating 
a period of more zonal character. The follow- 
ing decades show further cooling off until 
1840 but at least the 1820’s seems likely to 
have been of a more north-meridional char- 
acter. The famous 1840's and 1850's are 
again periods of high variability and increasing 
temperatures; they are therefore likely to 
have been meridional decades with dominat- 
ing southerly flow. From 1860 to the periods 
discussed in foregoing chapter, the circulation 
was probably again of a more zonal character. 
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On the Formation of Tropical Cyclones 


By S. SYONO, Geophysical Institute, Tokyo University 


(Manuscript received November 1951) 


Abstract 


In the present paper, an attempt is made to obtain dynamically the criterion for formation 
of tropical cyclones. Main factors important for typhogenesis are vertical thermal instability, 
horizontal dynamical stability, horizontal and vertical dimensions of perturbations. First the 
physical picture of a typhogenetical process is given. Secondly a method which enables us to 
treat the latent instability dynamically, a polytropic change method, is proposed. Thirdly using 
this method, the perturbation equations for the perturbations of typhoon-type are obtained. 
Solving these equations, the following criterion for the development of perturbations of ty- 
phoon-type is obtained 


2 develop, 


S 
= Sonne 


neutral, 


where g, S, Z. n and ry are gravity acceleration, vertical thermal instability, absolute vorticity, 
number of loops of perturbations in vertical direction and number which gives horizontal 
dimension of perturbations in 100 km. The left hand side has a similar form as Richardson 
number for the formation of turbulence. Fourthly, modes of perturbations are investigated. 


I. Introduction 


The problem of the formation of tropical 
cyclones is one of the oldest problems in me- 
teorology. Numerous diverse views concerning 
this problem have been proposed. However, 
until now we have not obtained a conclusive 
picture concerning the formation of tropical 
cyclones which has been accepted universally. 
One of the reasons is due to lack of appropri- 
ate data and another is due to shortage of dy- 
namical tools. Since the early 1940's, however, 
numerous data have been collected in the 
tropics mainly by the meteorologists in the 
United States and sufficient data for the re- 
examination of the view concerning the for- 
mation of tropical cyclones have been accu- 
mulated. On the other hand, since 1940, the 
“Vorticity dynamics” has been developed and 
is ready for the establishment of the mathe- 
matical theory of the problem. In 1947, 
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S. SAWYER proposed a hurricane theory. His the- 
ory threw light upon the mechanism of the for- 
mation of tropical cyclones. But his theory is 
not conclusive for three reasons. In his theory, 
the horizontal instability alone is treated dy- 
namically but no dynamical consideration of 
the vertical instability is given. Secondly, in 
his theory horizontal dynamical instability acts 
as a principal factor for the formation. But the 
assumption that the initial field is dynamically 
unstable in horizontal direction can not be 
accepted because such a field is soon broken 
and is not observed in the atmosphere. 
As many meteorologists, the writer con- 
siders that the vertical latent instability is one 
of the main factors for the formation of tropical 
cyclones. On this point, E. PALMÉN (1948) 
elucidated the importance of the vertical la- 
tent instability by using data obtained at Swan 
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Island. Thirdly he discussed the components of 
the equations of motion separately, this con- 
ventional method elucidates some partial fea- 
tures of the complicated phenomenon but it 
is not sufficient to obtain the criterion for 
the formation of such a complicated dynam- 
ical system. The equations of motion, of 
continuity, and of change of state should be 
treated as a coupled system. In 1950, H. REHL 
proposed a model of hurricane formation. 
However, for similar reasons as above his mod- 
el is also open to discussion from the dynam- 
ical point of view. 

In the present paper it is attempted to ob- 
tain dynamically the criterion of the formation 
of tropical cyclones. There exist two theories, 
i.e. one is a thermal theory and the other is 
a dynamical one. In the writer’s opinion, the 
two standpoints should be combined. In fact, 
in the present theory, the thermal instability 
and dynamical stability are equally important 
factors. 


2. Physical picture of the model 


Here we shall describe the model, on which 
the present theory is based. It is well-known 
that at the initial typhogenetical stage, there 
occur numerous scattered local convective 
showers and squalls, but no systematic wind. 
Then we may assume that there is no mean 
systematic converging current within the re- 
gion, where the typhoon is about to form. 
The following attempt is made to explain 
the initial pressure deepening. 

The pressure tendency at any level z is given 
by the well-known tendency equation 


Op _ ae 
3 =SQu-8 div (ov) dz, 


where p, 0, g, w, v, t and z are the pressure, 
the density, the gravity acceleration, the ver- 
tical velocity component, the horizontal veloc- 
ity, the time and the vertical coordinate re- 
spectively. 

First we consider the case, when mean mass 
convergence or divergence into the region 
does not exist, i. e. 


i div (ov) dz=o, 


% DIFOHNE 


then we get 


Here we can put 
gow=g(ew).+g(ew), 


where (ow), is the mean upward mass trans- 
port in the form of moist air across a unit area 


and (ow), is the mean downward mass trans- 
port in the form of rainfall across a unit area. 
At the present stage, since numerous local con- 
vection cells are surrounded by downward 
currents, the mean upward mass transport in 
the form of air is expected to be very small. 
The latter term is generally neglected, but it 
cannot be neglected, when the intensity of the 
rainfall is as large as in the case of the typhogene- 
tical stage. In fact the pressure tendency is 1 
mm/hour corresponding to an intensity of the 
rainfall of 10 mm/hour. The pressure at z rises 


or falls depending upon whether (gw), is larger 


or smaller than - (ow),. — (ew), decreases rapid- 
ly with height, so that the level, where (ew) 
vanishes, may not be so high. This level is a 
level of constant pressure and is referred to 
as N. P. V. L. Above N. P. V. L., the mean 
pressure rises and below N. P. V. L., the mean 
pressure falls. 

As soon as a horizontal pressure gradient 
is formed at any level, horizontal motion is 
excited. Then if the mean pressure within the 
region at any level above N. P. V. L. rises, 
mean mass divergence takes place from this 
region. (Fig. I, 2.) 

Now we consider the next stage. At this 
stage, there is systematic outflow aloft and 
inflow below. Then the net upward mass 
transport (ow), increases and the downward 
mass transport — (ow), may be neglected com- 


pared with (ow) The tendency equation is 
written as follows: 


H [ee] 
7 ca de a ow 
= 8 ui div (gv) dz-g div (gv) dz +g (ow), 
z H 


where H is the level of non-divergence which 
is referred to as N. M. D. L. Above N. M.D. LA 
the tendency equation becomes 
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Pressure 


A = 


Fig. 1. Initial stage of typhogenesis. 


Pressure 
falls 


7) a == 
= -g | div (ov) dz+g (ow). 


The pressure at any level rises or falls accord- 
ing to 


H oo 
-s | div (ev) dz+¢ (ow) = | dv (ov) dz. 
z H 


The instantaneous level of constant pressure 
z is given by the equation 


H 00 
— i div (ov) dz+g (ow) =g | a (ov) dz. 
z H 


Above N. P. V.L., the pressure gradient pro- 
duced is negative, and below N. P. V. L. the 
reverse holds. 

Next we study the vertical distribution of 
the mass divergence qualitatively. 

The radial component of the equations of 
motion of an air-ring is given by 


Ns), I Op 
a O0 


and 
M=r(ve+ro sin y), 


where r, @, v,, vo, fand M are the radial dis- 
tance of the air ring from the centre of the 
typhogenetical region, the mean latitude of 
the air ring, the radial and tangential velocity 
components, the Coriolis factor and the ab- 
solute circulation of the air ring. 

First we consider the case, where no system- 
atic rotational motion exists; then the sum 
of first two terms on the right hand side of 
the equation above vanish. Then the radial 
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acceleration is equal to the pressure gradient 
force. The wind in this case is called Eulerian 
wind. | 

The mass divergence and convergence are 
due to the Eulerian wind. 

If the pressure gradient generating process, 
i. e. convection, does not continue, the pressure 
gradient initially generated by convective 
showers is soon annihilated by the Eulerian 
wind. On the other hand, if convection con- 
tinues, the outward pressure gradient above 
N. P. V. L. is maintained and the outflow 
aloft continues. At lower levels, inflow due to 
the Eulerian wind continues as far as the in- 
ward pressure gradient is maintained due to 
the continuous outflow aloft. If no systematic 
rotational motion exists, the outflow aloft is 
soon balanced by the inflow below and the 
deepening of the surface system will cease. 
No deep low pressure of typhoon-type is 
generated, even though the convection is 
violent. This is the reason, why deep pressure 
lows are not produced in cases of ordinary 
thunderstorms and showers. (Fig. 2). 

Next we consider the case, where system- 
atic cyclonic motion exists. The systematic cy- 
clonic motion is generated by a systematic 
radial current according to the law of absolute 
circulation. In this case the sum of the first 
two terms on the right hand side of the above 
equation are positive, i.e. the resultant force of 
the centrifugal force and Coriolis force acting 


N.P.V Le 


N.M.D.L. 


Laer | \ 


Fig. 2. Second stage of typhogenesis. 
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upon the air-ring is directed outwards. Even if 
the air-ring is a constituent of the atmosphere 
initially at rest, the resultant force does not 
vanish except at the initial position of the air- 
ring. Below the level of no pressure gradient, 
i.e. N. P. G. L. the resultant force and the pres- 
sure gradient force counteract each other. On 
the other hand, above this level they act to 
drive the air-ring outwards. Even at N. P. G. L. 
the air-ring is driven outwards by the result- 
ant force. From this consideration, we may ex- 
pect that every air-ring is influenced by a 
strong outward resultant force at higher levels, 
but by a relatively weak force at lower 
levels and consequently; the outflow aloft is 
larger than the inflow below. Then a net 
mass deficit is caused in the region and deep- 
ening of the surface pressure takes place. 
This is the typhogenetical stage (Fig. 3). 


Cm 


A, 


TET EURE TRE CUITE HE 


Fig. 3. Third stage of typhogenesis. 
A: pressure gradient force. 
B: resultant force. 


Next we consider briefly the role of hori- 
zontal dynamical stability . We consider two 
concentric air-rings, and denote the inner one 
by (1) and the outer one by (2). Let us 
compare the radial accelerations acting upon 
them at the same radial distance. The radial 
accelerations are 


SAS YO NO 


dv, dv, i Li Viens (22) =, 
De M3) | en 


If the pressure gradient does not change with 
time, the second term on the right hand side 
vanishes and we get 


dv,\ > (dv, = 
(Fi), = (Ge) SM 
If the initial radial coordinate of the ring (1) 
is r, and that of (2) is r9+dr9, we may put 


OM? 
Ory 


where Z is initial absolute vorticity. Thus we 


get 
(52) = (52) = 
Z— O. 
dis) 7-2 GO) a, 

If the initial distribution of motion is dynam- 
ically stable in horizontal direction, i. e. 
Z > O, the inward acceleration of the second 
ring is smaller than that of the first ring. 
Then if the pressure gradient is kept constant, 
the inflow decreases with time. Thus the 
upward mass transport decreases. Therefore 
the horizontal dynamic stability acts against 
the typhogenesis. 

Further, M increases with the initial radius 
of the ring, thus a perturbation of larger 
scale in horizontal direction is expected to be 
more stable. 

Next we consider the role of the vertical 
latent instability. We have discussed the ty- 
phogenetical process with reference to hori- 
zontal motion. The discussion above is con- 
cerned with the removal of air from the typho- 
genetical region, but the cause of the process 
has not been discussed. In the author’s opinion, 
the cause of typhogenesis or the source of ener- 
gy for typhogenesis is the latent instability. As 
will be shown later, a high level anticyclone 
is developed as the result of typhogenetical 
circulation. From our analysis, however, we 
cannot say whether all high level anticyclones 
are the results of typhogenetical circulations or 
not. But even if a priori a high level anti- 
cyclone may initiate the typhogenetical proc- 
ess, it could not be more than a trigger action. 

E. Palmén showed that a necessary condi- 
tion for typhogenesis is a deep conditionally 
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unstable layer. The authors result supports 
this idea. The typhogenetical circulation is 
driven by the vertical instability, so that the 
deeper the instability layer becomes, the strong- 
er becomes the power of the typhogenetical 
engine. In fact, our solution gives the result 
that a deep perturbation is more unstable 
than a shallow one. At upper levels, air- 
rings are driven outwards by the resultant 
force of the pressure gradient force, the centri- 
fugal force and the Coriolis force. The air-rings 
which are driven outwards are those trans- 
ported from below by the vertical current cau- 
sed by the vertical instability. If the layer is 
vertically stable, the vertical mass transport is 
suppressed and the circulation will soon cease. 
In fact, we know that the typhoon weakens 
by the inflow of cold air. 


3. The perturbation equations 


We may assume that axial asymmetry does 
not act as the main factor for the formation 
of the typhoon. Using cylindrical coordinates 
(r, O, z), the equations of motion are 


DE Er Cent» 
de Ovo Ovo v,vo 

ci Es A aoe ee arated (2) 
dw aw aw eee (3) 
eon” er 0 dz À 


and the equation of continuity is 


do d(ov,) ov, d(ow) _ 
a ps (4) 


The equation of polytropic change is 


dp__ do 

FER (5) 
or 
OC I LA 2). 
a N ae do) 
where o is defined by 

o=kRT. (7) 


Here k is a polytropic constant; R the gas 

constant of the air; T the absolute temperature. 
As the initial state, we assume no vertical 

and radial components of the velocity and uni- 

form relative vorticity. Then we get 
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v,=W=0, (8) 
= do vo 
=; 0) 
and 
her: 
Dt (10) 
The equations of motion reduce to 
Fy ed 
ue r 0 or gar 
Te 
Le AGE (12) 


The symbol — denotes the initial values and 
the symbol ’ denotes perturbations. 

For the sake of convenience, we shall define 
some quantities 


2e ENT 2_ f2 
Be; HA. Ai r= Kr, (13) 
dve 12 = 
PE Ge (Cu 


where J, K and L are assumed to be constant. 
Now we consider the perturbations. Let 
V;=0,+0',, Vo=Vo+v'e, W=W+u", 
p=p+p,0=0+0" (14) 
The perturbation equations of (r), (2), (3), (4) 


and (6) are easily obtained by taking their 
variations. We get 
Ge gl, OV NET = OU 


+0 +0 +w +w to 
ot MORE Or Oz dz 


_2vev'9 I Op’ 1 op 


r 0 oO or er (15) 
vo y Ove = d'e eae do AE Ov (2) 
Da Be) 


tt ae - 
op) 
"292° (17) 
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00 I(0'v,) Pov) ev, ov’ d(v'w) 
a or u or De nr di 
alow’ 


Aw’) _ 8 
no (18) 
DE u Ge we Ls le 
et orale De D 
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u Kan Bei, PRE PES AE (19) 
Taking account of the assumptions for the ini- 
tial state, the equations above reduce to 


vo 
a + Zv’,+ Lrw’ =o, 


Ow if Op 1%, Tl Opa 
at =(-3 a De) (Eee) (22) 


90° Nor‘) , Bu’, Now) 


à Pr a (23) 
op’ Op _ de DECO 
at (3 ae (2-02 3 

do’ 
=0 + (24) 


Introducing m,, mo and m; defined by 
M,=0V',, Mo=0V'o, Mz=0w’, 


these equations become as follows: 


om, dp’ ; 

ome Zmo= Tee Kro’, (26) 
Ome 

+ Zm,+ Lrm,=o, (27) 
Om; DDR, 

a we (28) 

Jo’ Om, m, Om; 4 

TE ern © (29) 
op m. (dp  d0\ m. (op — a 
on (F-« )* 0 (3 % se) 

Jo’ 
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4. Polytropic change method 


It is well-known that the latent heat of wa- 
ter vapour contained in the air is one of the 
energy sources of the atmospheric motions. 
It is an important but difficult problem in 
dynamical meteorology, how to treat the 
condensation energy of water vapour in the 
dynamics of the atmospheric motions. We 
shall here use a conventional method for this 
purpose. 

First we consider the polytropic change of 
an air-parcel of unit mass. 

The equation connecting p and v in adiabatic 
processes 1s 

(31) 


where y and v are the ratio of specific 
heats and the specific volume respectively. 
Generalizing this relation to the case where y 
is not a constant but a parameter which is 
independent on the state variables and using 
k instead of y, we get 

(32) 


This is the equation connecting p and v in 
polytropic processes. In accordance with the 
potential temperature, the polytropic poten- 
tial temperature ©; is defined by 

k—1 


= T (£2) (33) 


This quantity was introduced by the writer 
in 1941. 

In order to make clear the physical meaning 
of k, we consider an example, i. e., an air-parcel 
of unit mass is ascending polytropically in the 
atmosphere. 

From the first law in thermodynamics, we 


get 


pv’ =const, 


pvk= const. 


dq=C,dT* + Ap*dv*, (34) 
where dq is quantity of heat supplied to the 
air parcel during the process and the asterisks* 
indicate the quantities of the ascending parcel. 
If the heat supplied per unit time is Q, dq= 
=Q dt. Using the equation of state 

PrE=RTLK (35) 


we get 


(36) 


If the parcel ascends with a constant velocity 
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dz > 
x. di= x. The hydrostatic equation for the 


surrounding atmosphere is 


dp= -godz, (37) 
so that 
dz d RT d 
fe ee oe a 
ac gow* gw* p (38) 


Introducing (38) into (36), we get 


Be Mra Tra RT N 
Ge So aang wee 
We can put p=p* and T is approximately 


equal to T* for the small interval of the height. 
Thus 


dT*_ AR( __Q \ dp* 

ips Cp on, p* (40) 
If Q and w* are independent of the state varia- 
bles, we may integrate this equation and get 


AR Q 
T*p* Cr (2 >= ae —CONSE (41) 


Eliminating T* by using (35), we get 
AR RT 
Cp (: an) lis 


Comparing (32) and (42), we get 


PEAR (<Q igs RQ 
k Cc; en MAC (c.+ (43) 


OT 


yk p* =const. (42) 


= RQ 
k=C, (c, +) (44) 
If the parcel ascends without addition of heat 
Q=o. Equation (44) then becomes 


» 2235 


= (45) 
If the parcel ascends with addition of heat 
Q= gw*/R, the equation (44) becomes 


Rat. (46) 


The former case means adiabatic change and 
the latter means isothermal change. 

When the parcel ascends moist adiabatically, 
the latent heat of water vapour is released. In 
this case Q depends on the state variables 
and the ascending velocity is not kept constant 
and the equation (40) is not integrable. But 
for a small interval of height, they may be put 
constant. Following a similar line of argument 
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as above, we get the polytropic constant for a 
finite interval. If the sign of Q/w* is not 
changed, we may define similarly, the poly- 
tropic constant for a finite interval of height 
by using the appropriate mean value of Q/w*. 
We shall call this process “pseudopolytropic 
change”. 

Now we return to the problem of the ty- 
phoon. In the main part of a typhoon, there 
exists an ascending current in which latent heat 
is released. On the other hand, if the circula- 
tion is closed, descending currents exist in the 
surrounding parts, in which the air is heated by 
dry adiabatic compression. Then k must be 
different in these two parts. Corresponding to 
these different values of k, the sign of the ver- 
tical stability changes, i. e., in the main part 
it is negative and in the surrounding part it is 
positive. Since the vertical stability in the 
surrounding parts acts against the vertical in- 
stability in the main part the circulation does 
not develop if the typhoon is treated as a 
closed system. However, the atmosphere is in- 
finite in the horizontal direction, the air trans- 
ported to the upper levels by the ascending 
current may be driven outwards. Then it is 
not reasonable to assume that the vertical 
stability prevents the development of a circula- 
tion. However, it is desirable to take into ac- 
count the spatial distribution of the vertical 
stability but in order to treat this strictly we 
encounter a mathematical difficulty. For the 
sake of simplicity, we assume that the verti- 
cal instability is constant throughout the space. 
Strictly speaking, this assumption corresponds 
to a fluid heated from below uniformly and 
does not represent the model depicted in sec- 
tion 2. The criterion obtained below under 
this assumption may be considered to give the 
most favorable conditions for typhogenesis. 
The exact theory may be laid for future in- 
vestigation. 

Using the equation of state, we get 


k—ı k—ı 
FN ante LEO (aay aa 
STR) Ae, 
k—I ul 
k nk 
=" a (47) 


Differentiating (47) logarithmically with re- 
spect to r, we get 
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1 OO, _1 RR) -s 8 
zu ame (48) 


and similarly 


1 204 _1 SB) -s 
ler: BE (49) 


Here S, and S, are polytropic stabilities in the 
directions r and z respectively. Through the 
assumption mentioned above they are con- 
stant. 


1 op k 3 
Gee Janis ah 
1 (dp 05 -£(L 25) s 
aes 3 o\p azo dz) %% 6) 
Then (30) becomes 
Po (S;m+Sem.)=0 À (52) 


5. The perturbations 


Now we assume m,, Mo, mz, p’ and 9’ are 
proportional to e*. Then 
dp’ 
or 


am,—- Zme= ——=—+ Kro’, 


ame+ Zm,+Lrm,=o0, 


, 


aM, = — = mvs 


1 Ar Mr Omg A 
ax rt) 
8 Dir coz 


(56) 


(57) 


The perturbations now considered develop 
during a period of several hours, so that « 
is of the order 10-3 ~ 105 sec—!. The vari- 
ous orders of magnitudes of other quantities 
in c.g.s. unit are summarized below: 


ap’ +o (S,m,+S;m;)=x00'. 


0=0(10~3), Z=0 (10-4), 

o=0 (10°), K=0 (109), L=o (10-10), 

S,=0 (10%), S,=o (107), (58) 
m,=0 (1), me=0 (1), mz=0 (1071) 

p’=0 (rot), o'=0 (10-5), 


SIISYON® 


am, & ame & Zm, x Zme = 0(1073), 
Lrm, = 0 (10%), 
’ D) / 
—=0 (1073), Po (10m2), 


dz 
Kro’=o (10%), go’ =o (1072), 
ao’ = 0 (108), xop’=o (10), 
ap’ = 0 (10), «op =0 (10), 
oS,m, = 0 (1), 0S;m;=0 (10). 


(59) 


This list gives the maximum values of the 
various quantities, and thus a comparison be- 
tween different terms in the perturbation 
equations must be made with caution. We 
get the following equations: 


am,- Zmo= - = + Kro’, (60) 
amo + Zm,+Lrmz=0, (61) 
= — 20’, (62) 
om, m, dm 
r fe ze 6 
Dre ED (63) 


ap’ +0S;m,;= 000’. (64) 


It is of some interest to elucidate here the 
physical meaning of neglected quantities. From 
the third equation «&o’ is neglected, which is the 
time variation in the density, From the fourth 
equation &S,m, is neglected; this means that the 
effect of the horizontal thermal stability is 
small compared with that of the vertical ther- 
mal stability. However, if the vertical stability 
is very small, the horizontal stability may be- 
come effective. From (60) and (61), we get 


7) , 
(02+ Z2) m,= 2% +aKre -LZr mz, (65) 


(«+ 22) mo=Z fe ZKro'- Larm,, 


Ir (66) 
or 
dp’ 
ME an aKro’-bLrm,, (67) 
op’ : 

mo=b Ei bKro’-aLrm,, (68) 

where 

rar: Z 

Tuer tar ge (69) 


From (62) and (64), we get 
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—p,, (70) 


where 


g=-7- (71) 


Hereafter the subscript z of S. is dropped. 
Introducing (70) into (67) and (68), we get 


„op aK +bLq | op’ 22 


th, = — an 5 D tp; (7) 
ep bK- alq, Op aLq , 
mo = b Dr + 7 oe = wise (73) 


Introducing (70) and (72) into (63), we get 
2 , 
ep = (£-4r) oe + Br “2 


en 


oz? “Or Or dz 
+ Ep'=0 


= 2 2 
= tee 4 c=* +Z 
oS gs gS 
1 (2(KS- LZ) | a+ Z?) 
S | £ o 


D= 


Now we assume 


p=R(ne (76) 


Then the differential SE for R is given by 
| dR | 
-(4 -ißB 10 ma)? a oe 
D 2 = 
(2 -or)r- 


a ft _ (KS-LZ_, (KS-LZ)* _ 
= r 20S gS (a? + Z?) 
_KS- i | dR o+Z7f ,, 
=> )r J dr Sg \P 5 
£ Es LZ 2gLZ 
a0 +2: o( 


(77) 


: tae 
In the second term of this equation, = is lar- 


ger than o (10~°) and 
KS-LZ (KS-LZ) =) - 
( 20S si S(a? + Z?) gS 
is smaller than o (10-9); thus we may neglect 
the latter compared with the former. We get 
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PEL RR=0, (79) 
where 
nn 
pe er (50) 


The equation (79) is the Bessel’s differential 
equation of zero order. The solution which is 
finite at r=o is J, (kr). 

As the equation (80) contains ß in the form 
B, it is invariant with respect to changes i in the 
sign 7 ß. The solution of (74) is given by 


p=e 3c" (A, cos Bz +B, sin Bz) Jy (kr) e*. (81) 
Introducing (81) into (70), we get 


tm=[las(£-2) nantes 


+ A,ß+B, (£-2) sin 


D 
xe 2C* 7, (kr) et, (82) 
where 
D'ART DE Tarai 
2C 20 @+Z2 5 2C 20 a2+Z2 


The boundary conditions at z=o (the earth 
surface) and z= H (the tropopause or the upper 
boundary of the unstable layer) are that the 
vertical velocity vanishes. Thus 


ID) 
nee) (6) 
ß,H=nn. (84) 
Then (82) becomes 
Dn 
m,=A, sin Baze 2Cn° Jo (kr) ent, (85) 
where 
al. (£- el Kr ae é 
A,=- +B, 
z e | A,B 6 2C, LE of, HI 
fp NOVA Ts 
(Er) le (9 
Then (81) becomes 
p= ae OS (Anz + En) 
Ze; 
n\, eats 25 PE 24 72 
Dy 
a Cy ls (kr) ent (87) 
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where 
2 LZERS 
20 a?+22 
fan LZ-Ks\® 
2 FR Zu 
i+ (Ls 1) 
COS &,= Bn (88) 


‚(8,228 
Vz * (£+ a? +2? 
6. Stability of the perturbations 


In the present section we shall study the 
stability of the perturbations. The stability cri- 
terion is obtained from. (80). 

Solving (80) with respect to «? + Z?, we get 


a2 + Z?=-y„(IFVI-6,) (89) 
eZ) 


Vn= 2 = 
RER, 
(8 a a) 


_ 290 (ko S+ KS+LZ) 
Wea ae 
9? 
CPS 
. 277 
? (ee) 
[07 
(40° Bi +g?) (KS-LZ)? 
g? (RE So+ KS+LZ)? 


Here we may neglect some small quantities. 


where 


On = 


(90) 


From (84) B,=n a where H=1.5x 105 cm, 


(the thickness of the troposphere), then we get 
Bn & 2n:10-6. g=0(10°), then 4ß2o?=n2x 
o(107) and g?=o (10°), thus g? may be neg- 


lected compared with 48202. k=o (10-7), K= 


K 
=0 (10-9), k®=0 (10-%), no) thus 


= may be neglected compared with k2S. 
Then (90) reduces to 
__g (koS+LZ) 
Vn 20 B2 > 
5 40°82 (RS-LZ) 
gt (RoS+LZz} 


(91) 


Ss. SYONO 


I. Barotropic atmosphere L=o 


When the atmosphere is barotropic, the 
vertical gradient of the velocity vanishes, then 
the vertical gradient of the vorticity also van- 
ishes. Thus the roots «, are given by 


Dal — Z2—y, (1+V1-6,), 
ao Az (I+V1-ö,), 
ada 22, (1-V1-6,), 


Ot = Ver (1 wi VI — bn), 
where 


(92) 


(93) 


Since the first root of Jp (kr) is 2.4, we may 
take 1 as the characteristic radius of the 
perturbations. ry is of the order 10° cm. 
We get k==4 x 107, where r, is given in 100 
km. K is ites than 1078, then 

On = 0 (1072). 
Then the roots «, become as follows: 
oft) = 4/ SRE 
An Er Er 
«=, 


SEE YA 


(94) 


Corresponding to the roots «(3 and «1, we 
get neutral perturbations of inertia oscillation 
type. When «, is nearly equal to + iZ, the 
equation (79) may not be considered as an 
approximate equation of (78), because the sec- 
ond term with a factor r in the coefficient of 


dR 
Cannot be neglected. In the present 


discussion, this case is of no interest. 
If-Z2-2y,<o, the former roots al) and 
an” give also neutral perturbations. If- Z2- 
—2Yn > 0, the first root «fl? gives unstable 
perturbations and the second root a?) gives 
damped perturbations. The perturbations cor- 
responding to the first root are of particular 
interest. The stability criterion is given by 


(95) 
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A N) 
Table I. Magnitude of $ number at T—280° defined by nr ÔT represent the instability in terms of the 


temperature difference between the moist adiabatic lapse rate and the actual lapse rate per 100 m. 


OT! = 

RSS 0.01 0.03 0.05 0.10 0.15 0.20 0.25 0.30 
TOR 3:5. 202 07.050.102 07.752702 sa Too S2SE702 7.0870. 08.750102 77.088708 
Se Oz 1.4 10% | 4.20 10% | 7.00 ıo®| 1.4 104 | 2.10 10% | 2.8 ro4 3.50 104| 4.20 104 
LORS 3.5 10% | 12.05 10% | 1.75 10% | 3.5 108 5.25 10° |. 7.0 10% M8.7s ro®| 1.05 104 
S-1o-* 1.4 10% | 4.20 ıol| 7.00 ıol| 1.4 ro? | 2.10 ı02| 2.8 102 3-5 0.702 || 208702 
or more unstable, when the S number increases. 
NS = BE _ 6 n2 2 unstable, (06) ne = the vertical is Reve S 
72 = L2 0 neutral. number becomes negative, then pertur ations 


The left hand side has the similar physical 
meaning as Richardson number. We call it 
“S number”. 

From the above criterion we may draw the 
following conclusions. 

(i) If the S number and the vertical dimen- 
sion of perturbations are given, the flatter per- 
turbations are most stable. 

(ii) If the S number and the horizontal di- 
mension of perturbations are given, the deeper 
perturbations are most unstable. 

(ii) If the horizontal and vertical dimensions 
of perturbations are given, perturbations are 


are always neutral. The magnitude of the criti- 
cal S number for various Z and ÖT are given 
in Table 1, where 6T means the temperature 
difference between the moist adiabatic lapse 
rate and the actual lapse rate. Table 2 gives the 
critical horizontal dimensions for various S 
numbers represented by ÔT. The perturba- 
tions larger than those given in Table 2 are 
stable. 


From Table 2 we can see that perturbations 
of typhoon type can develop for plausible 
values of vertical instability and that too strong 
shear prevents the development of perturba- 
tions. 


Table 2. The critical horizontal dimensions in 100 km for various z and values of the instability 
represented by 6T°. 


n |Z Ce 0.01 0.03 0.05 0.10 0.15 0.20 0.25 0.30 
Tor seGe = 
I 21.46 36.02 47.98 67-90 83.11 | 96.02 107.29 TES, 
5 4.29 eas 9:60 13-57 16.62 19.19 21.46 23.50 
I Io 2.135 3.72 4.80 6-79 GPA 9.60 10.73 INA 
20 1:06 1.86 2.40 3.40 4.16 4.80 Sa 5.74 
50 0.43 0.74 0.96 1.36 1.66 102 DAMES 2.35 
I 10.73 18.01 23.99 33-95 41.56 48.01 52.64 58.79 
5 DTS 3.72 4.80 6.79 8.31 9-60 10.73 AT 
2 Io 1.08 1.86 2.40 3-40 4.16 4.80 5-37 5.74 
20 0.54 0.93 1220 1.70 2.08 2.40 2.69 28 
50 0.22 0.37 0.48 0.68 0.83 0.96 1.07 1.18 
I 7.15 12-01 15.99 22.63 27.70 32.01 35.76 39.19 
5 1.43 2.48 3.20 4.52 5.54 6.40 GA 7.83 
3 Io 0.73 1224 1.60 2.26 2.97 3.20 3.58 3.82 
20 0.37 0.62 0.80 Ted 3 1.39 1.60 1.79 1.91 
50 0.14 0.25 0.32 0.45 0.55 0.64 0.72 0.78 
I een 9.01 12.00 16.98 20.78 24.01 26.82 29.39 
5 1.07 1.86 2.40 3.39 4.16 4.80 5.37 5.88 
4 TO 0.54 0.93 1.20 1.70 2.08 2.40 2.63 2.87 
20 027 0.47 0.60 SES 1.04 7.20 1.32 1.44 
50 O.II 0.19 0.24 0.34 0.42 0.48 0.54 0.59 
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oT? ro 0.5 | I | 2 | 3 | 4 | 5 | 6 7 | 8 | 9 Io 
n 
I 10 000 2 500 625 | 278 156 100 | 69 | SI 39 | 31 25 
2 2 500 rh el TER TEN TR | 13 | ED. 
; 1100 | 278 | 69 ie EG DE | lé aan nde 3 
4 750 a 39, Bez] al 6 IE RE 2 
I 20 000 5 000 | I 250 | 555 | 313 | 200 | 278 | 102 | 78 | 62 | 50 
2 $ 000 1 250 | | 139 | 78 Bra: 16 | 13 
0.1 3 2 220 555 | 139 | 620| 335 | = 22 | 15 | Il | 9 CARE 6. 
4 | 1250] 313 78 15 | V7) Taam EN A 
5 800 | 200 so | 22 | 13 | 8 | 6 4 3 | 3 | 2 
I 40 000 Io 000 2 500 | ig Geren | 625 | 400 | 278 | 204 | 169 | 124 | 100 
2 Io 000 2 500 625 || 278. | 1567| 1001 CORNE EME Re 
0.2 3 4 445 TUE arr nr ET | SR EST AE 17 14 II 
4 2 500 625 169 70 38 25 | 17 | 13 | Io | 8 6 
5 1 600 400 | 100 45 | AS A 6 | 8 6 5 4 
I 60 000 15 000 3750 I 666 | 937 | 600 | 417 | 306 | 234 | 285 | 150 
2 15 000 3 750 938. |BEs56 | 234 | 150. Rz Fr | 59 12407738 
3 6 665 1 666 | 417 185 | el] 34 267 | 127 17 
Oe PR 3 750 938 | 234 [ein] MON ee) | ae 9 
5 2 400 600 | 150 Er 38 RE: 17 12 9 7 6 
6 1 667 417 |) 104 46 DO we | 5 9 | 7 5 4 
II. Baroclinic atmosphere L =o the vertical stability. In our case S is negative, 


Nextuwescohsider the: efece of thesverti- 50 that the positive verte (000 


aha to prevent the development of the perturba- 

her order of le is“ 1079" ör smaller tions and vice versa. If S is very small, k2So 

LZ S 0 (1073) and k?So=o (10-12). Thus becomes comparable with LZ. If L=10 9 

o(LZ/k2So) is 107} or smaller. (10 m/sec/t 000 mxı00 km), Z is 1074, k= 

We get the following approximations: =10°7, o=10°, k2|Sjo becomes equal to LZ 

4K2p2 LZ when S=10-8. When T=280°, the difference 

== ae ( aoe ) EZ <\KS| between the lapse rate and the moist-adiabatic 

grk lapse rate is 0.03/100 m. If k=2: 10 8, the differ- 

0,20 SLZERS ence becomes 0.7/100 m. In the ordinary case, 

TER (98) L is much smaller than 1079, and S is compar- 

6 = 16 —— :-LZ=KS able with o.1/100 m. Then we may assume 
orks that k?So is larger than LZ. 


… BeL?Z? 2KS\. From this discussion, we can conclude that 
GS apa se D :LZ> |KS| the vertical shear has small effects upon typho- 
& genesis in normal cases. 


222 127282 

Since 4 Br <-e"{(107 6) and ee = 
g2kt g?ktS2 7. The rate of development of perturbations 
0 (10° 5), 6, may be neglected compared with 1. In the present section we shall study the 


The effect of the vertical shear is given by rate of the development of perturbations. For 
yn. As may easily be seen from (91), the effect the sake of simplicity, we shall consider a 
of the vertical shear is in the same direction as barotropic atmosphere. 
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The rate of development is given by a), 
The amplitudes of the perturbations are in- 
creased by a factor e, when t,=1/a, 

From (93) and (94), we get 

| a ee 
a A -Z +elS| 35. (99) 

In c.g.s. unit, we assume T=280°, k= 

2: 

a 1077, By=n* 2.0: 10-6 and g=980. Thus 
0 
we get 


(100) 


where ÔT is the temperature difference in °C 
per 100 m between the lapse rate and the 
moist-adiabatic lapse rate. In Table 3, the values 
of 5 ÔT/n?r? 10-6 in sec-! is given for various 
OT, n and r,. Since the rate of change, «1? 
given by (99) contains four parameters Z, dT, 
n and ry, it is not easy to give the values of 
t, for various parameters in a compact form, 
thus we confine ourselves to obtain the range 
of ty. 

Table 3 is divided into five parts with four 
lines: double lines, thick lines, thin lines and 
dotted lines. These lines separate the unstable 
perturbations from neutral perturbations for 
various values of Z, i.e., Z=2' 10~4sec~1, 1074 
sec-!, 5:10 5 sec“! and 2:1075 sec-"!. The 
regions to the left of these lines correspond to 
unstable perturbations. The entries of the table 
are 525 10~° in 10! sec. Thus if we take 
7 in 10-5 sec! and subtract Z? from the 
entries of the table and take the inverse square 
roots of the result, we get 1, in 10° sec. Putting 


N= Rate 7 100) 1019, we get 
BT, 


I 
un 105 sec= — x 2746". 
VN VN 


Table 4 gives 1, for various N. 


I 
ze h 
Table 4. tn = VN x 274 46%, 


oT 
N = TE rer 10-*} 1020, 


n’r; 


N| 0.01 | I | 100 | 1.000 | 10.000| 50.000 


2h46m | 53265 | 16mgs | 7m3s 


in | 2778 | 27h46m 
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Comparing Table 3 and Table 4, we find 
that perturbations of typhoon-scale may de- 
velop under very favorable conditions. Further- 
more, perturbations of the horizontal scale of 
a typhoon develop very slowly in the atmos- 
phere for n=z. 


8. Modes of the perturbations 


We shall here investigate the modes of the 
perturbations. 

The vertical velocity perturbations are giv- 
en by (85). 


Ante -—z 
= sin B,ze 2Cn. J,(krjetnt. (101) 


At lower levels within the typhoon, w is pos- 
itive and thus A, is positive. It follows that 
the vertical motion is directed upwards all 
through the troposphere when n=1. For n=2 
we get ascending motion in the lower part 
of the troposphere and descending motion 
in the upper part. For perturbations n> 2, 
ascending and descending currents alternate as 
shown in Fig. 4. 


tropopause 


Von 


n=3 


AIN WA 


Fig. 4. Vertical velocity distribution. 


From the investigation in 7, we see that the 
perturbations of higher values of n only can 
develop in the vicinity of the centre. Since the 
relative magnitudes of the perturbations of 
various order are determined by the initial per- 
turbations and since the magnitude of the type 
n=2 is larger than that of n=1, a descending 
current in the upper half of the troposphere 
may exist in the vicinity of the centre. 

Next we shall consider the pressure perturba- 
tions. They are given by (87) and (91), 


gSA, 
/ 2 
J pa [ S KS 
an, Br (+) 


a = Jo (kr) em, 


+ En) e 2Cn 


L 


DS 


cos (B,2+ 


(102) 
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where 
(103) 


Since S is negative and A, is positive, the 
constant coefficient of p’ is negative i. e. with 
an increase of the strength of the vertical cur- 
rent, the pressure at lower levels within the 
typhoon decreases. The vertical velocity per- 


, ERS EIN 
turbations vanish at z=- H (i=o, 1,...n). On 
n 
the other hand, the pressure perturbations 
; 1 
vanish at levels Hi lower than- H, because 
n 


En 18 positive. At these levels Hj, the sign of 
the pressure perturbations change. 

The heights H} of the nodes of the pressure 
perturbations are given by 


Br Hi+e,= ae eT 
or 
> _ ET K 
en Dos os Boy S 
| 2n mm nn \2002y, mn 
10 
_Hj2i-ı I, Ne | : 
ARS not 20° Hgk? 
For k= 1.17, T = 260°, R =2.87- 10°, o = 
=9: 108, we Be = 510°”. On the other 
2 2K 
hand, for K = 1079, Tes = SiO, top ati 


Thus for unstable perturbations, which have 
small values of n and ro, the latter term may 
be neglected in comparison with the former. 
Then (104) becomes 


km 


tropopause 


15 


10 


n=l n —2 n=4 


Fig. 5. Perturbed pressure distribution. 


n 2 TT 207 N 


Table 5 gives Hi (km) for various i and n. 
And Fig. 5 shows the change of the pressure 
perturbations with height. 


Table 5. H; (km). Heights of the nodes of 
pressure perturbations. 


4 I 2 3 4 5 
n 
I 6.38 
2 3.47 10.97 
3 2.37 Weil 12:37 
4 1.80 5.55 9.30 13-05 
5 1.46 4.46 7.46 10.46 13.46 
D 


n 

: A er + 
All perturbations contain a factor e 2Cn°. This 
factor shows that the perturbations decrease 


exponentially with height. Let us investigate - 


this factor. From (75), we get. 


D, gS [2KS | a+ Z| 
20 2 (a2 +22) S\ g oie 
KS g 
a2 Ze ua toe (106) 


For the same reason as in the case of H;, the 
first term may be neglected. Then the factor 
is the same for various values of n. Further- 


more, en $°10 7, and thus we get 
o 


sr ER + 7 
e 2Cn =e 20 =e 


(107) 


Table 6 gives the values of this factor for every 
km from z=o to z=15 km. 

Among other things, the vertical distribu- 
tion of p’ for n=T is interesting. In connection 
with the deepening of the low-level pressure 
system, the strength of the high-level high in- 
creases (H. RIEHL, 1948). 

Let us compute the ratio of the amplitudes 
of p’ for n=1 at z=o and z=H. From (102) 
and Table 6, we get 


De _cos (B.H+e;) te à 
Pr COS €; 3 = 
= -0.47. (108) 


We shall finally investigate the horizontal ve- 
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Table 6. Vertical diminishing factor. 
z (km) lo are ne er ete orc 
e-5-10 'z IT | 9:95 | 0.91 | 0.86 | 0.82 | 0.78 | 0.74 | 0.71 | 0.67 | 0.64 | 0.61 | o ‚58 | 0.55 | 0.52 | 0-50 | 0.47 
Bas perturbations. Introducing (102) into (72) PR. On _ mn \/ _Z2—9S Re . 
and (73), we get a2 + Z2 27n gSk? 5 
kA, 7 ZB? 
ov,= - gShA, cos (Anz + be an (111) 
\/ ; (2 =) a2+Z? Sk? 
2Vn ß} nr (de et & 5 : 
20 2Yn and introducing (111) into (109) and (r10,) 
an we get 
+En 2Cn™ kr) ent £ 
IX N) cos (Br&+E&n) x 
254; (= bLq D, u) © Be: 
Sr EE = n 
pe ARS: C2 CG 
ain \/ ass SEE en 
en OCT Er) EE 
Dn 
RE LT nn s 
r cos (B,z+e)e 2" Ji (kr) e®n = FE = 1 cos (Bnz+En) x 
r i 
GSBnAz aK+blq, sin (B,z+ 20k? \ Be an 402 
za nl 8 Me 
ae, \/ Bis 2002, xe 2 ae r) emt 
gun Kb PA t 
+ ) e 2Cn 95 (kr) en (109) un r sin (Br? +€n) 2 
= n ie 
E a gk? \/8: es 
OV'e— KS z s (Pnz+ ' 
Dn CA a z 
Bm cay een 
ars 
gSA, (ES D; 9) a p 40? 
x Cc ee " 
Vie ae er Ga) ew 
A n co 2yn Kp2 ZA r 
Te See DE = cos (Buz +8) X 
meostpreteye "7 (kr) ent 20%,k? \/6: pee) 
n 40? 
gSPnAz aK+bLg rsin (BP, z+ eye 
2 KS & ete ey (kr) em 
= BÈ + Vn K 3 ZA 
Dn Bn = rsın (Bnz+ En) 
LE n° iB en ie 
ae Jo (kr) er. (110) an gk? \/ 2 Moi 
Putting i 
D, 9S k? me ze 20° 7, (kr) e (173) 
L=o, 2C, 20° ne Be’ where 
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The second and third terms in (112) and (113) 
may be combined as follows: 


2 
ov, = nn cos (Bu2+ En) X 
B & 
Senne 

awz re: 

tie 
een Jere 
KB? Az : Te 


sin (B,z+28,) e 2° rJy(kr) ent, (115) 


kg 


= nAz 
OV'e—= p Pa 
Aw 2 eh ds 
n 40? 


cos (Br&+ En) x 


2 
ae! sin (nz + 2&) x 


REG 


we oe 
ee oo (ET) e116) 


The ratios of the first terms to the second 
terms in (115) and (116) are 


KB} Az 
keg 


EEE 
apa 
it er 
ib 2 FE kg 
Aw 
=0 (1075). 


Thus the second terms are negligibly small 
compared with the first terms. We get 


OV, = Pr Az = cos (Baz + En) x 

AVE 

er 
es eier) en, = fer) 

© 2A, 
Ve = = a cos (ByZ+En) x 

= / 2 Sas n 

kV pb ar 40? 


24 

xe 22° J, (kr) ent (118) 

From these results we see that the strengths 
of radial and rotational motions increase with 
the development of perturbations. The direc- 
tions of these currents are convergent and 
cyclonic in the low pressure regions and vice 
versa. The vertical distribution of v, and vo‘ 
is the same as that of p’ shown in Pigs. 


Ss. SYONO 


The speed of development of the cyclonic 
motion is proportional to the initial absolute 
vorticity and inversely proportional to the 
developing coefficient x. | 
a, decreases from a finite value to zero with 
the increase of the absolute vorticity. If %,=0, 


7 | 5 : 
the factor — becomes infinite but for this case 


le . 
the perturbation is of a finite intensity after 
an infinite time. 
Using the formula of the Bessel function 


5 (x J, (x)) = x Jo (x), the vorticity perturba- 
x 


tions are given by 


2 
c= 2.2 (me) = oe ee (Bn +&n) 
Br EN 
RoW Bat a 
Bore d GA (Rr) ane 
ror 
Pr Az 


ae cos (Bnz+&n) 


ay BS ‘ 


a ER 
‚Io (kr) en 2 em (arg) 


From (119), we see that at lower levels the 
cyclonic vorticity is concentrated in the in- 
nermost region r < r,, and outside this area a 
region of negative vorticity is produced. On 
the other hand at upper levels, the negative 
vorticity region is surrounded by a region 
of positive vorticity. 

Since the velocity components and the 
vorticity are inversely proportional to 9, their 
magnitudes increase with height, even if the 
pressure perturbations are constant. 

The monotonic change with height of the 
velocity components is given by the factor 

a 
ie 2, Relative magnitudes of this factor 
at every km are given in Table 7, where 
the value at the earth-surface is taken as unity 
and © is the density of the standard atmos- 
phere. 


Concluding remarks 


In the present paper we have studied the 
perturbations which develop as the result of 
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Table 7. Vertical diminishing factor of velocities. 


Re Er ee rr re ns nn es Se 


z (km) | o | I | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | nr | 12 | 13 | 14 | 15 
Qo,-£ 
7} 20 1.00 | 1.05 | 1.10 | 1.16 | 1.27 | 1.30 | 1.38 | 1.47 | 1.56 | 1.68 | 1.80 | 1.94 | 2-08 | 2.41 | 2.68 | 2.99 


vertical instability and we have obtained nu- 
merous plausible results. 
However, the following questions remain 
for future investigations. 
1. In the equation for the pressure perturbation 
(78) we neglected the second term in the 


aha BASE it 
coefficient of —-, this neglection is justified 


dr’ 

when &+Z?2 & 10-8 sec-2. But if « tends 
to +iZ, the neglection is not valid be- 
cause the neglected term contains a term 
which is inversely proportional to «?+ Z2. 
From our neglected equation we obtained 
perturbations of inertia oscillation-type. 
These perturbations are interesting, but we 
cannot use the solution without a detailed 
discussion of the exact equation. 

2. We have discussed the components of per- 
turbations for various n and k. However, 
in order to approach to the actual typhoon, 
we must combine these components by 
using Fourier series and Fourier-Bessel’s 
integral theorem. Corresponding to four 
roots a1), a2), «(3 and «4, four initial 
conditions for p’, v,, ve, w are given. If 
the solutions corresponding to a and a,” 
were correct, we could trace the develop- 
ment of the model-typhoon from any 
given initial state. 
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3. For the sake of simplicity, we assumed a 
fixed tropopause, but we can treat the 
actual atmosphere, by taking a two-layer 
atmosphere, where the upper layer is stable 
and the lower layer is latent unstable. 


4. For the sake of simplicity, we assumed that 
S is constant and the atmosphere initially 
unstable but in order to treat the actual 
process of typhogenesis, both assumptions 
must be abandoned, and at least the sign of 
S must be changed for the ascending and 
descending parts. At present, however, this 
problem is very difficult, because we cannot 
determine the regions of ascending and 
descending motion and they change during 
the development of the perturbations. 
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in a Blocking Situation 
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Abstract 


The synoptic situation under study is a blocking situation over Western Europe and the North- 
Eastern Atlantic. Barotropic tendency computations for the 500 mb surface are carried out and 
compared with the 12 hour observed height change. Tendency computations are also carried 
out for the mean height (averaged along the vertical with respect to pressure). The vertical 
motion compatible with the observed local vorticity changes is computed from the vorticity 
equation at several levels. The result is in fair agreement with the observed distribution of pre- 
cipitation and cloudiness. The computations indicate that horizontal divergence and horizontal 
vorticity advection are of about equal importance in the vorticity budget at soo mb, whereas 
the vertical vorticity transport seems to be less important. 


The equations governing atmospheric mo- 
tion are very complicated, and the amount of 
useful information about future development 
which we have been able to derive from the 
“exact” equations is certainly small. In recent 
years, however, one has attempted to overcome 
the difficulties by a simplified treatment of 
atmospheric motion, putting the emphasis on 
the largescale features and neglecting all de- 
tails. This line of attack has been advocated by 
RossBy, who in 1939 proposed to consider 
the atmospheric motion as horizontal and 
non-divergent, so that the motion is charac- 
terized by the absolute vertical vorticity bein 
conservative. CHARNEY (1949) has found that this 
property should be attributed to the motion at 
a certain “equivalent barotropic level” in the 
vicinity of the 500 mb level. This very simple 
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“barotropic” model has proved to be very 
fruitful; it leads to equations which we are 
able to handle and which explain a consider- 
able part of the observed height changes at the 
500 mb surface. On the other hand, it is clear 
that the barotropic model represents an over- 
simplification and has severe defects. How- 
ever, by correcting the model for the more 
important of these defects, we may hope to 
arrive at other models that give a more satis- 
factory description of the observed changes 
without adding too many complications. From 
this point of view, it seems to be important 
to study the errors due to the barotropic mo- 
del. This is not a very easy undertaking, for 
the discrepancies of the computed changes are 
caused not only by the shortcomings of the 
model itself, but are due to two other sources 
of errors as well, viz. errors in observation and 
analysis, and errors which are introduced in 
the process of numerical integration. 

For theoretical reasons, we should expect 
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the barotropic model to be most unsatisfac- 
tory in cases of strong development. Another 
case where the model seems inadequate is in 
certain quasi-stationary motions. Some exact 
solutions, representing stationary, horizontal 
and non-divergent motion on a rotating 
sphere are known. (Craic 1945, HÖöILAND 
1950). They have the character of waves of 
finite amplitude with troughs and wedges 
symmetrical with respect to the meridians. 
On the other hand, semi-stationary motions 
with a pronounced tilt of the trough and ridge 
lines relative to the meridians are often obser- 
ved. It is difficult to see how such motions 
could be explained from the conservation of 
vorticity. An example which frequently occurs 
in connection with a blocking situation is 
characterized by an almost zonal flow over 
the North Atlantic, bending sharply around a 
cyclone in the vicinity of the British Isles 
and continuing as a southerly or even south- 
easterly current towards Iceland and Eastern 
Greenland. Such a streamline pattern some- 
times persists for one or two weeks or more. 
If we assume that the air particles move hori- 
zontally through such a stationary field, we 
must conclude that their vorticity will increase 
during the transition from the zonal to the 
cyclonically curved flow. We should therefore 
expect important individual vorticity changes 
to take place in this region. 


Fig. 1. soo mb map at 1500 GMT, 24 May 1951 showing 
the gridpoints used in all calculations. Heights in deca- 
meters. 
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On the basis of such considerations, it was 
thought that a closer investigation of the vor- 
ticity budget in a blocking case might throw 
some light upon the nature of the deviations 
from the conservation of vorticity. Through 
such an investigation, some clarification of the 
dynamics of the blocking phenomenon might 
be realized. 

Two different courses of action seem to be 
possible. One might investigate the mean 
vorticity budget over a longer period by means 
of maps averaged with respect to time, or 
on the other hand, the analysis could be based 
on actual maps. In the former case, the eddy 
transport of vorticity would enter into the 
considerations as an important unknown fac- 
tor. Because of this factor the latter method 
was chosen. This method, in turn, also pre- 
sents complications; namely, that the vorticity 
changes connected with smaller scale wave 
perturbations that travel through the semi- 
stationary field will tend to confuse the long- 
trend effects. However, this difficulty is inevi- 
table since we cannot a priori be certain that 
these smaller scale motions are not part of 
the large scale mechanism. 

A blocking situation of the type discussed 
occured from the 19th to the 29th of May, 
1951. In its main features, the field of motion 
over Western Europe and the Eastern Atlantic 
was almost stationary during this period. The 
transition from zonal to cyclonic flow took 
place over the British Isles where the aerolo- 
gical network is good. The twelve hour period 
from 1500 GMT on the 24th of May to 0300 
GMT on the 25th was selected for the study. 
The 500 mb map at 1500 GMT is shown 
in Pig or. 

At first, a tendency computation was per- 
formed on the basis of the barotropic model, 
assuming the local vorticity changes at the 
500 mb level to be due to advection only, as 
expressed by the equation 
fawn (1) 
dt 
Here € means relative vertical vorticity, f time, 
v horizontal wind vector, V horizontal del- 
operator and n absolute vertical vorticity, 
which may be written as 


n=Ct+f (2) 


where f is the Coriolis parameter. 
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When the geostrophic wind approximation 
is introduced in (1), a Poisson equation for the 
height tendency results. This equation was 
treated as described in a report by the Starr 
MEMBERS, UNIVERSITY OF STOCKHOLM (1952). 
Employing a square grid on a conformal map, 
as shown in Fig. 1, the right-hand side was 
determined by means of the difference method 
for the 24th and 25th, and the height tendency 
was determined by relaxation. 


Using a method outlined by Bouin and 
CHARNEY (1951), the tendencies thus obtained 
at the two instants were multiplied by 6 hours 
and then added to obtain a 12-hour height 
change. The result is shown in Fig. 2, which 
was computed under the assumption that the 
height tendencies along the boundary of the 
rectangular area were zero. The correlation 
coefficient between the computed height 
changes and the observed changes in the same 
12-hour period for 49 grid points in the center 
of the area was 0.84. Fig. 3 shows the computed 
height changes when the actually observed 
changes were used as a boundary condition 
instead of zero values (compare UNIVERSITY OF 
STOCKHOLM 1952). Somewhat surprisingly, the 
correlation coefficient between observed and 
computed changes was only 0.70 in this case. 
This does not necessarily imply a marked 
depreciation in the results, but rather demon- 
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Fig. 3. Computed s00 mb height changes from 1500 
GMT, 24 May to 0300 GMT 25 May in decameters. 
Observed height changes used as a boundary condition. 


strates the inadequacy of judging these tests 
using correlation coefficients alone. The ob- 
served 12-hour height changes of the soo mb 
surface are shown in Fig. 4. 

The figures demonstrate that the computa- 
tions give rise and fall areas in approximately 
the correct positions. The most pronounced 
discrepancy, however, is that the fall area to 


Fig. 2. Computed s00 mb height change from 1500 
GMT, 24 May to 0300 GMT 25 May in decameters. 
Height changes at the boundary assumed to be zero. 
Mean soo mb pattern shown by thin solid lines. 


Fig. 4. Observed soo mb height change from 1500 
GMT, 24 May to 0300 GMT, 25 May in decameters. 
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the SE of the low center is calculated to be 
too weak. This error is undoubtedly due to 
the barotropic model; the assumption that 
local vorticity changes are due to vorticity 
advection only is too restrictive to give good 
results in this region. 

The “barotropic model” as designed by 
CHARNEY (1949) was based on the assumption 
that the vorticity equation averaged in th 


vertical direction may be written with suffi; 


clent accuracy 


oe — ' 


HV VU D-0 0 


The application of vorticity conservation to the 
motion at a certain “equivalent barotropic 
level”, which may be approximately represen- 
ted by the soo mb surface, however, rests on 
certain additional assumptions concerning the 
variation of wind with height. It is conceivable 
that the errors are mainly due to these assump- 
tions, and that (3) is a reasonably good approxi- 
mation. If the wind velocity is sub-divided 
into a mean wind v and a thermal wind vr, 


V=V + vr C=0+ lr (4) 


equation (3) may also be written 


Fig. 5. Mean map for 1500 GMT, 24 May 1951. Isopleths 
represent (2 — 700) in decameters. 
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Fig. 6. Computed 12-hour mean height changes from 
1500 GMT, 24 May to 0300 GMT, 25 May in deca- 
meters. Observed height changes used as a boundary 
condition. Mean Z pattern shown by thin solid lines. 


When geostrophic winds are used, it is readily 
seen that v is related to the mean height z 
in the same way as the wind at a particular 
level is related to the height of the correspon- 
ding isobaric surface. Therefore, eq. (5) can be 
used to determine the local change of the mean 
height in the same manner as the changes of 
the 500 mb surface were determined, provided 
that we are able to evaluate the right-hand side 
of the equation. An attempt was next made to 
test this equation. Maps of the mean height 
were constructed for the 24th and the 2sth of 
May by means of the following formula, the 
derivation of which is given in the Appendix: 


2 =) 0.075 27000 OS Ter 0 i 0.175 
2700 + 0.2 2509 + 0.15 Za00 + O.I 
Zgoo + 0.15 2199 + 0.05 Tiooo + 
+ 0.06 T700 — 0.24 T300 + 2.68 Tioo (6) 


The subscripts denote the pressures to which 
the heights and temperatures refer. In the 
computations, T'isoo and Ts) were replaced 
by constant mean values since their variation 
is seen to be of very little importance. In 
every gridpoint, the heights and temperatures 
at the different levels were read from the 
carefully analyzed constant pressure charts, 
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Fig. 7. Observed 12-hour changes in the mean height z 
from 1500 GMT, 24 May to 0300 GMT, 25 May in 
decameters. 


and the mean height was computed. Fig. 5 
shows the mean map for the 24th of May at 
1500 GMT. 

At first, the advection of mean absolute 
vorticity with the mean wind was determined 
from the two maps, and a 12-hour height 
change was computed, disregarding provi- 
sionally the last term of (5) and using the 
same relaxation technique as before. The re- 
sult is shown in Fig. 6, and the corresponding 
observed changes of the mean height are shown 
in Fig. 7. The correlation coefficient between 
these two change patterns was found to be 
0.76. As at the soo mb level, the computed 
falls to the SE of the center were too small. 

An attempt was then made to include the 
last term of (5). This may be done under the 
assumption that the thermal wind vr does 
not change direction with height and that the 
change in magnitude is similar in all places, 
so that we may write vr= Ap) VT 1000 Where 
VT 1000 is the surface value of the thermal wind. 
This is in fact the assumption on which the so- 
called 2%-dimensional model is based (com- 
pare ELIASSEN 1952). Equation (5) may then 
be written as 
at 
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The geostrophic value of the surface thermal 
wind can be determined from the map of 
Z—Z1000 (which is actually a map of the mean 
x aa) We are 


temperature since Z — Z1000 = a 
D 

thus able to evaluate this term under the 

assumptions mentioned, provided a good 

estimate can be given to the weighting factor 

A. 

The term was first computed assuming that 
A (p), is a linear function, so that A?= 1/3. 
The corresponding corrections to the mean 
height tendencies were again determined by 
relaxation. The computed 12-hour mean 
height changes, with this correction included, 
are shown in Fig. 8. The principal contribution 
of the thermal wind correction was found to 
be an increase in the intensity of the area of 
falls to the S and E of the center, which de- 
finitely improved the result. The inclusion of 
this correction increased the correlation 
coefficient between observed and computed 
mean height changes from 0.76 to 0.84. 

Actually, the function A(p) deviates con- 
siderably from linearity. It has been estimated 
from a series of ascent curves from Berlin by 
S. Smebye, and the values are given in an 
article by ELIASSEN (1952). If these values are 


used, then the weighting factor A2 is found 
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Fig. 8. Computed 12-hour mean height changes from 

1500 GMT, 24 May to 0300 GMT, 25 May including 

thermal wind correction with weighting factor A? = 1/3. 
Units are decameters. 
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Fig. 9. Computed mean height changes from 1500 

GMT, 24 May to 0300 GMT, 25 May including thermal 

wind correction with weighting factor A? = 1. Changes 
in decameters. 


to be closer to unity. When the computations 
of the 12-hour changes of the mean height 
were repeated with the values 2/3 and 1 for 
the weighting factor, the corresponding corre- 
lation coefficients between computed and ob- 
served changes were 0.86 and 0.89, respectively. 


Fig. 9 shows the computed changes for A? = 1. 
As far as can be judged from the correlation 
coefficients, this computation gives a_better 
final result than that based on the soo mb map. 
This is confirmed by comparing the computed 
changes shown in Fig. 9 with the observed 
changes shown in Fig. 7. Even the falls to the S 
and E of the center of low pressure show up 
in the computations with the correct strength. 
This improvement should not be thought 
of simply as a result of the inclusion of the 
thermal wind correction, for the simple baro- 
tropic equation crudely incorporates this same 
term when it is assumed to apply not at the level 
where v=?, but rather where v )=v? (com- 
pare CHARNEY, Fyortorr and von NEUMANN 
1950). However, when eq. (1) is used, the 
thermal wind term is assumed to be a constant 
fraction of the total advection of relative 
vorticity, whereas it is computed indepen- 
dently using eq. (5). 

Making use of the vorticity equation in a 
diagnostic sense, an attempt was next made to 
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determine the three-dimensional field of 
motion responsible for the observed vorticity 
changes. When pressure is used as a vertical 
coordinate and friction is neglected, the com- 
plete vorticity equation may be written (see eq. 
(s) in ELIASSEN 1952) 


dés CNE OOM al 
imine ul aa ip 127 Vo 


(8) 


Here p denotes pressure, « specific volume, 
and w the individual rate of change of pressure, 
which is a measure of vertical motion. More- 
over, \7 means here differentiation along an 
isobaric surface, and 0/dt rate of change in an 
isobaric surface. Obviously, the local change 
of vorticity is due partly to advection and 
partly to non-advective changes represented by 
the last three terms, which denote vertical trans- 
port of vorticity, the effect of horizontal diver- 
gence and transformation of horizontal vor- 
ticity into vertical vorticity, respectively. The 
expression for the geostrophic wind has been 
introduced into the last term. These non- 
advective vorticity changes are seen to be 
closely related to the field of w, i.e. to the 
vertical motion. In fact, the equation may be 
used to determine æ when the local and 
advective changes are evaluated from the 
synoptic maps. Thus the three non-advective 
effects may be calculated separately. 

The idea to compute the vertical motion 
from the vorticity equation is not new; similar 
calculations were made by DEDEBANT and 
WEHRLE (1933) and by SAWYER (1949). 

When & is considered as the unknown 
function, eq. (8) is a linear, first-order diffe- 
rential equation. It determines the change of 
æ along the characteristic lines, which are 
seen to be the absolute vortex lines. Knowing 
the distribution of ® over a certain surface, 
i.e. the earth’s surface, the three dimensional 
field of vertical motion may be determined by 
integration along the characteristics. 

In a barotropic atmosphere \/a =o, and the 
absolute vortex lines are vertical. (Strictly 
speaking, they are parallel to the earth’s axis; 
the discrepancy, which is unimportant, is 
introduced by the neglect of the horizontal 
component of the earth’s rotation in the deri- 
vation of (8)). The atmosphere is, however, 
baroclinic in most regions, and the vortex 
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lines are therefore sloping. Generally, the 
inclination between the vortex lines and the 
vertical is of minor importance. Strongly 
baroclinic regions, for example in frontal zones 
where the absolute vortex lines are roughly 
parallel to the frontal slope, are an exception. 
In the situation considered, the slope of the 
vortex lines was found in most places to be 
very steep, so that their horizontal projection 
was much smaller than the grid-size used. It 
was therefore decided to neglect the slope 
altogether and to consider the characteristics 
as vertical. This corresponds to neglecting the 
I 


rer a : Vo of eq. (8), which simplifies 


the computations considerably. Upon doing 
so, eq. (8) may be written in the form 


où 
ie 2) DAT (9) 
DENE 0 
which may immediately be integrated to 
obtain 
2 Ie 
6 
CRCRS ER 
N D=p1 1 P—Po 7 P 


pi 


(10) 


Letting po denote the pressure at the ground 
where the value of w is known from the 
kinematical boundary condition, this equation 
may be used to determine w at the arbitrary 
level py. 

In the computations, the integral in eq. (10) 
was evaluated by means of the trapezoidal 
rule from values of es + v : V7 obtained from 
analyzed maps of the 1,000, 850, 700, $00 and 
300 mb surfaces. For instance, the integral 
for the layer from 1,000 mb to 850 mb was 
approximated by 
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and similarly for the other layers. From these 
expressions, the values of ® were determined 
at 850, 700, soo and 300 mb in all gridpoints, 
assuming w =o at 1,000 mb. This is justified 
because the 1,000 mb surface is close to the 
ground, and the value of w at the ground 
(which is assumed to be level) is much smaller 
than in the free atmosphere. 

The local vorticity changes at the various 
layers were approximated by the 12-hour 
change from 1500 GMT the 24th to 0300 
GMT the 25th of May. This change may be 


taken to represent the value of En at 2100 


GMT on the 24th of May (in the middle of 
the interval). In order to make the various 
terms of the equation applicable at the same 
instant, the vorticity advection must also be 
determined for the middle instant of the 12 
hour interval. Since there are unfortunately 
no upper air maps available at 2100 GMT, 
the computation of the vorticity advection 
was based on mean maps formed by averaging 

the maps for 1500 GMT and 0300 GMT. 
The values of -» computed in this way are 
shown in Figs. 10a. to 10 d. To a certain 
extent, the results can be corroborated by 
comparison with Fig. 11, the surface map for 
0000 GMT the 25th of May, which also shows 
the areas where precipitation was reported 
during the 12 hour period from 1800 GMT 
the 24th to 0600 GMT the 2sth. The precipita- 
tion areas over land are accurate; however, 
those over water are based on sparse data. 
The correspondence between the computed 
vertical motion and the precipitation areas 
is quite good, and on the whole, the computed 
vertical motion agrees with what one would 
expect in such a synoptic situation. The warm 
front south of Greenland, for instance, is very 
discernible in the 850 mb vertical motion. 
Weak descending motion in the warm sector, 
stronger ascending motion ahead of the front, 
and descending motion again far ahead of the 
front (especially at higher levels) all agree 
with the usual picture. The downward mo 
tion to the west of the principal low center 
and over Norway are also in agreement with 
reports of clear or partly cloudy skies. The 
old occlusion across Scotland is likewise il- 
lustrated by weak ascending motion associated 
with reports of light rainfall ahead of it. The 
correspondence over western France and the 
Tellus V (1953), 2 
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Fig. 10. Twelve hour mean vertical velocity — m in tens of millibars/6 hours from ısoo GMT, 24 May to 
0300 GMT, 25 May at (a) 850 mb, (b) 700 mb, (c) soo mb and (d) 300 mb. Corresponding mean constant 
pressure heights shown by thin solid lines. 


Bay of Biscay is somewhat poorer. As a whole, 
the results illustrate the feasibility of computing 
vertical velocities through the use of the 
vorticity equation. The computations were 
here based on observed local vorticity changes; 
however, if such changes could have been 
correctly predicted at several levels, the com- 
putations of vertical motion could have been 
based on them as well. Therefore, if we some 
time in the future will be able to compute 
satisfactory prognostic maps for several iso- 
baric surfaces, we will also be able to determine 
the corresponding vertical motion accurately 
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enough to be directly applicable in weather 
forecasting. 

In the 2 %-dimensional model, w is assumed 
to have a very simple distribution along the 
vertical, characterized by a maximum value 
at a certain middle pressure level, and zero 
values at the top of the atmosphere and at the 
ground. The computed @ values show very 
little tendency towards such a simple distribu- 
tion. However, one cannot draw too definite 
conclusions from the computations regarding 
the vertical distribution of m, because the 
errors in the computed values must be expec- 
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ted to increase rapidly with height due to 
doubtful analysis at higher levels. The com- 
puted w values have in some places a maximum 
as low as 600—700 mb; in other places, no 
maximum is reached even at 300 mb. The 
latter is true in several grid points in the 
area to the SE of the low center. As far as 
this result can be trusted, it indicates that the 
level of non-divergence is here actually located 
above 300 mb, which may explain the poor 
results of the barotropic tendency computa- 
tions in this region. 


From the computed values of w, the indi- 
vidual terms of the vorticity equation were 
evaluated. As one would expect, these terms 
show quite ragged distributions. In order to 
reduce these random errors, the values were 
smoothed by averaging over a square con- 
taining nine grid points; the averaged values 
may be interpreted as circulation changes 
instead of vorticity changes. The result for 
the soo mb surface is shown in figures 12 a 
to 12.d. 


A comparison between these figures shows 
that horizontal divergence and the vorticity 
advection seem to be of equal importance in 
accounting for the local vorticity changes, 
whereas the vertical vorticity transport seems 
to be of lesser importance. In fact, the mean 
absolute values in (6 hours)-!/12 hours were 


0.221 for D 135 for — v : V/ 7, 0.044 for 
od 
7 
—@ Wy and 0.175 for 7 = 


In judging these results, however, one 
should have in mind that the errors resulting 
from imperfect observations and analysis, as 
well as from the geostrophic wind approxima- 
tion will tend to increase the absolute magni- 
tudes of the computed values of w and the 
divergence term. The mean absolute values 
dr 
op 
fore likely to be overestimates. 

To the W of the low center, the advection 
seems to account for the greater part of the 
local vorticity decrease observed in this region, 
but the local vorticity increase observed SE 
and E of the center must be almost entirely 
ascribed to the effect of horizontal conver- 
gence. Therefore as far as the calculations can 


do . 
for — — and 5,7 given above, are there- 
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Fig. 11. Surface map for 0000 GMT, 25 May showing 
reported precipitation areas during the period from 1800 
GMT, 24 May to 0600 GMT, 25 May. 


be trusted, the indications are that the vor- 
ticity increase necessary to keep the field of 
motion semi-stationary is produced by this 
horizontal convergence. 


Appendix. 
The mean height Z is defined as 


Po 
= I 
= —|zd 
z af p 


where py is put equal to 1,000 mb. Integrating 
by parts, and using the hydrostatic equation, 
this equation may be written 


(12) 


Po 
= ROUES 
== Z1000 + Fa T dp (13) 


Thus z may be determined simply by graphi- 
cal integration of the ascent curve, plotted 
in a pT-diagram. 

However, if the values of z and T at several 
isobaric levels are already known, we may toa 
certain extent utilize this information, and 
thereby save work. 

Suppose that we know from the aerological 
report that the isobaric surface Po has the height 
% and the temperature Ty; the isobaric surface 
pı has the height z, and temperature T,, and 
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Dr’ (b) — V 7; (c) —@ > and (d) 7 7 at 500 mb obtained by averaging 
over nine gridpoints. Units (6 hours)~*/12 hours. Mean 500 mb pattern shown by thin solid lines. 


Fig. 12. Smoothed values of (a) 


so on. The integral in (13) is then broken up Z (po) = Zo 
into a sum of integrals over the layers between 
these levels. Consider one of these integrals, z (p)=2ı 
over the layer from p, to pı- Within this , Role 
layer, we approximate z by a polynomial in = (Po) = mer (15) 
p of the third degree, $ Po 
ET. 
z(p) = A+ Bp + Cp’ + Dp? (14) 2° (pi)= en 


where A, B, C and D are constants. This func- These equations determine the values of A, B, 
tion must satisfy the four conditions C and D. When the values of the coefficients 
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thus obtained are introduced into (14), and 
the integration is carried out, we obtain 


Po 
I 
Jen p0) + 5 Ble? — pr) = 
pı 


I 
+ ac (Po? — pr?) + 7 D (po! — pr’) = 


Pay A FE 
> Poe Pr) 
TER Er dU 
en 6 
le) (9 


This formula was applied to the layers 1,000— 
850 mb, 8s0—700 mb, 700—500 mb, $00—300 
mb, 300—200 mb and 200— 100 mb. Clearly, 
the same formula does not work for the 
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uppermost layer, 100—o mb. For this layer it 
was instead assumed that the temperature is 
constant and equal to Too; this gives 


100 


f? dp = (210 + : Tun) 100 mb (17) 


0 


When the integrals are added, we finally ob- 


tain 
Z= 0.075 Z1900 + 0-15 Z850 + 0.1752700 + 
0.2 Z590 + 0.15 Z390 + 0.10 Zggo + O.15 Z100 + 
+ - (0.00188 Ti999 +0.00208 T75) — 
— 0.00833 T399 +0.09167 Tjoo) 


(18) 


which is identical with (6). 
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Abstract 


A general discussion of multiple parameter models of the atmosphere for forecasting pur- 
poses is given and various ways of arriving at such models are presented. For some two-parameter 
models it is possible to derive a simple graphical method of integration which may be used 
before high-speed electronic computors are available. A numerical example is given. Finally 
the accuracy of a two-parameter representation is investigated more closely and a comparison 
is made with the accuracy of the radiosonde data. The conclusion is reached that a two- 
parameter model is adequate for representing the field of pressure-heights with an accuracy 
which is only slightly less than the accuracy of the radiosonde observations themselves. Hence 
models of three or more parameters are of interest merely to the extent that they are 
capable of describing in a more satisfactory manner the processes of the atmosphere. 


I. Introduction 


The barotropic model of the atmosphere has 
proven to give a very useful first approxima- 
tion to the actual motion of the atmosphere 
for computational purposes. The few complete 
numerical forecasts that have been published 
up to now show that the success of the baro- 
tropic model is comparable with the success of 
standard methods for forecasting at use in the 
weather services. However, the computations 
are extremely voluminous as is seen for 
example from the article by CHARNEY, FJoR- 
Torr and VON NEUMANN (1950). Therefore 
the method seems not applicable in prcatice 
until we have large-capacity high-speed elec- 
tronic computors at our disposal. 

One of the main reasons for the large 
amount of computations required is the necess- 
ity to fulfil the Courant-Friedrich-Lewy’s cri- 
terion for computational stability. This defines 
the maximum time interval At over which 
extrapolation in time may be carried out, 
when the grid size As is given, In order to have 


1 This work has been supported by grants from the 
Knut and Alice Wallenberg Foundation, Stockholm 
and the U. S. Weather Bureau, Washington D. C. 
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a fairly good resolution the grid-size must 
not be chosen too large and At turns out to be 
1—2 hours. This seems to be a surprisingly 
short time period to a practical meteorologist, 
who is accustomed to extrapolate over periods 
of 24 hours with considerable success. It is 
quite obvious that these time extrapolations are 
done in different ways. In the former case a 
strict Eularian method is adopted and no direct 
use is made of the fact that the absolute vor- 
ticity is conserved. In the practical extrapola- 
tion, on the other hand, certain quasi- conser- 
vative fields are moved over the map and the 
local effects are computed by integrating in 
time at various points on the map. The ques- 
tion arises if a similar method of integration 
can be advised for the barotropic model. This 
has lately been done by Fyorrorr (1952). By 
his method one person can make a complete 
24-hour forecast for an area of about 107 km? 
in three hours using nothing but graphical 
methods of addition and subtraction. The ex- 
trapolation in time is carried out in one 
time-step through the introduction of a certain 
mean stream field which varies slowly and in 
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which the vorticity field may be advected. One 
might expect, however, that even if the vor- 
ticity field is advected by the actual wind this 
may be done over time periods of 6—12 
hours before it is necessary to compute a new 
stream field. The forecast equations of the 
barotropic model may be written. 


À jm (1) 
n =? et (2) 
on 
n(t + At) =n(t— At) 4 an 2At 65) 


where z is the height of the pressure surface 
defining the motion, 7 is the absolute vorticity, 
g is the acceleration of gravity, fis the Coriolis 
parameter and J (1, «:) = dx,/0x da,/Py — 
— 0%,/0x 9a, /dy. In order to obtain 2n/df we 
must know the two fields of and z. However, 
the relative changes with time of z are con- 
siderably slower than corresponding changes 
of 7 and it may therefore be possible to use 
(1) and (3) alternatively over a few time steps 
without computing new values z from (2) at 
every instant. This procedure in which the 
vorticity field is advected in the initial 
wind field over several elementary time 
steps may be utilized for computations with 
a 2-parameter model. In this case it does 
not seem be possible to derive a slowly vary- 
ing mean field in which some invariant 
quantity is advected, as was done by Fjor- 
TOFT (1952), but we may thus try to use the 
actual wind field itself and extrapolate over 
6—12 hours before deriving the new stream 
field. Such a method will, of course, not be 
of immediate practical value in the way 
Fjortoft’s method for the barotropic model is, 
since we cannot make a 24 hour forecast in 
one step but have to use two or three, and 
each time step will mean 8—10 hours’ work 
for one person. Still, it may be of some interest 
for preliminary computations before high- 
speed electronic computors are avarilable and 
a few examples will be given here. 


1 The reader is referred to an article by CHARNEY 
and Phirrips (1953) for a detailed discussion of these 
problems. 
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Before we discuss this method in detail it 
is of some interest to approach the multiple- 
parameter models of the atmosphere for 
computational purposes in a somewhat more 
general way than has been done previously 
(PHILLIPS 1951, ELIASSEN 1952, EADY 1952, 
ARNASON 1952). 


2. Multiple-parameter models 


In all models of the atmosphere advanced 
so far for computational purposes, one makes 
use of the geostrophic approximation in order 
to express the wind field in terms of the pres- 
sure field in the vorticity equation. A set of 
dynamically consistent equations incorpora- 
ting the geostrophic approximation was derived 
by CHARNEY (1948) which implies essentially 
an elimination of oscillations of the inertia- 
gravity type that develop in connection with 
the adjustment of the wind and pressure fields 
to each other. This is probably a good approxi- 
mation when we deal with large-scale quasi- 
barotropic phenomena, but becomes more 
questionable when we try to incorporate 
baroclinic processes which are of a smaller 
scale. We shall, however, in this analysis adopt 
this approximation and leave the question of 
the goodness of the geostrophic approxima- 
tion for further investigations We shall here 
also neglect effects of friction and non-adia- 
batic processes. 

The multiple-parameter models now es- 
sentially represent attempts to simplify the 
complicated three-dimensional atmospheric 
structures and behavior patterns by intro- 
ducing simple assumptions about the distri- 
bution of the dependent variables with res- 
pect to the vertical coordinate. The name of 
the model indicates the number of param- 
eters used in order to describe the field of 
motion. This representation is used partly 
to simplify the computations and partly to 
get more accurate expressions for the deriva- 
tives along the vertical than finite diffe- 
rences will give, in particular if only a few 
gridpoints are used. As our vertical coordinate 
we may choose z (height), p (pressure) or © 
(potential temperature) or any function of 
these three. The dynamic and adiabatic equa- 
tions become considerably simpler if pod 
are chosen. In the former case dz/dt and w (= 
dp/dt) will be the unknown quantities, which 
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are determined from the vorticity equation 
and the adiabatic equation, in the latter case 
Op/ot (where y = c, T + g2) is our only de- 
pendent variable and may be determined from 
the vorticity equation. The adiabatic equa- 
tion is automatically fulfilled by our choice of 
coordinate system. It is well-known, how- 
ever, that the boundary conditions at the 
ground become considerably more compli- 
cated in that case. We shall use pas the vertical 
coordinate in the analysis below. 

The two basic equations that govern the 
motion are the vorticity equation and the 
adiabatic equation 


(5) 


v is the horizontal velocity, and k denotes the 
vertical unit vector. The continuity equation 
has already been used to eliminate the diver- 
gence term from the vorticity equation. The 
boundary conditions at the top and at the 
bottom of the atmosphere are (o denotes 


density) 


@=0 at p=o (6) 
d a 
re Op pi (7) 


Equation (s) may be transformed with the 
aid of the hydrostatic equation and we obtain 


where, of course, 9 In O/dp also is a unique 
function of z. 

In order to solve the two equations (4) and 
(8) with the boundary conditions (6) and (7) 
we may now follow the procedure used by 
ELIASSEN (1952) and Eapy (1952) and assume 
certain parametric representations of z and 
w and thus also of dz/dt: 
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ap + > Zn (x, 91)" An(p) 
(9) 


(10) 


ne) 


5 Om (XP) “Be (D) 


m=I 


œ (es y; p> t) = 


Zn and @,, are functions to be determined from 
the field of motion, while À, ( p) and By, (p) 
are power expressions of p given once and 
for all. H(p) is the height of the various 
pressure surfaces in the standard atmosphere. 
Both Eapy and ELIASSEN use N = 2 and 
M = 1 (however Eapy has z as the vertical 
coordinate instead of p). The boundary con- 
ditions are simplified to © = o at p = po and 
p =o in Ettassen’s model and Eapy puts 
w= 0 at the top and the bottom of the 
atsmosphere. The simplest function B,, (p) 
satisfying this boundary condition is a parabolic 
function of the vertical coordinate. 

Another possibility is to solve for w from 
equation (8) and introduce in (4) which then 
contains only one unknown function dz/dt. 
The boundary condition (6) must be trans- 
formed in a similar way in order to eli- 
minate w. It is now sufficient to introduce 
one assumption such as (9), which has to 
satisfy the vorticity equation and the boundary 
conditions. It is of course also possible to 
eliminate Jz/dt, but in that case we have to 
use both assumptions (9) and (ro). We shall 
proceed as EADY and ELIASSEN. 

It is quite obvious that the approximate 
expressions for dz/dt and w given by (9) and 
(10) cannot satisfy the equations at all levels. 
At each level where we require them to be 
satisfied we obtain one relation containing 
dz,/ dt and w,, but N + M represents the maxi- 
mum number of conditions we may derive 
from the two equations (9) and (ro) and the 
two boundary conditions without overdeter- 
mining the problem. It is also seen from equa- 
tion (4) that the expressions we finally obtain 
for dz/t cannot give an exact solution. The 
nonlinear terms in the equations will contain 
products of the functions A, (p) and B, (p). 
Thus dy/dt and also dz/dt as well as © cannot 
be represented by the functions A,(p) and 
B,, (p) originally chosen. Higher powers of p 
will appear but are automatically excluded by 
the method of representation. However, this is 
true for any numerical solution of a no1.-linear 
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equation of this type. The very fact that the 
function z according to the differential equa- 
tions very soon will contain all possible powers 
represents one of the major reasons why this 
equation cannot be solved by known analyti- 
cal methods at present. 

From a physical point of view it is not 
unreasonable to disregard these higher powers 
of p. It means that we assume that the energy 
originally contained in the large perturbations 
dissipates via the small ones but that no ener- 
gy transfer takes place in the opposite di- 
rection. Of course the same thing is true in 
the horizontal dimensions where we dis- 
regard the small perturbations which cannot 
be caught by the grid we use. Even if these 
harmonics were important for the develop- 
ment of atmospheric disturbances of the scale 
we observe in the atmosphere, it is very 
unlikely that for example the geostrophic 
approximation would be valid in this case. 
This assumption does not permit us to consider 
the detailed structure of the motion. 

It was indicated above that one way to 
obtain a sufficient number of equations for the 
determination of dz/dt and w is to assume 
that the two equations (4) and (8) as well as 
the boundary conditions are satisfied at a pro- 
per number of levels. However, this is not the 
most rational way, in particular not when 
M and N are small. Instead we may proceed 
in the same way as Eapy and devide the 
interval P = Po to p=o into a sufficient 
number of intervals and require that the two 
basic equations are satisfied if integrated over 
these intervals. Two conditions are again 
obtained by satisfying the two boundary 
conditions. Another method was employed by 
ELIASSEN who required that certain weighted 
pressure means of the vorticity equation and 
adiabatic equation should be satisfied if inte- 
grated over the whole interval p = po to 
p = 0. The question of how to determine the 
appropriate equations is important particularly 
when the number of parameters is small. 

In the case of the barotropic or one-parametric 
model it is quite obvious that the single equation 
needed has to be obtained from the dynamics, 
in other words from the vorticity equation. 
This is done by satisfying the equation at one 
level approximately in the middle of the 
atmosphere (the equivalent barotropic level) 
or by demanding that the vorticity equation 
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be satisfied after integration over the vertical. 
This model does not permit any conclusions 
concerning conditions at the upper and lower 
boundaires. 

The next model is one where two parameters 
are used to describe the field of motion, N = 2. 
M can be chosen arbitrarily to some extent, 
but little is gained if M > 2. (ELIASSEN, and 
similarly Eapy, puts M = 1. The 2 + M equa- 
tions needed to determine the time-derivatives 
of the parameters may again be derived in 
various ways. ELIASSEN (and also Eapy) has 
chosen his B, (p) so that the two boundary 
conditions at the top and the bottom of the 
atmosphere are automatically fulfilled and 
all three equations needed (N + M = 3 in 
this case) are derived from integral require- 
ments regarding the vorticity equation and 
the adiabatic equation as was indicated in a 
previous paragraph. In this way the model is 
fitted to the observed structure of the atmos- 
phere. However, it can be shown (see sec 3) 
that these integral conditions are equivalent to 
the requirement that the equations be valid 
at certain levels. These levels may not be 
chosen arbitrarily, but are determined by the 
character of the functions A, (p), A, (p) and 
B, (p). The two levels, where the vorticity 
equation is satisfied, may then be called the 
equivalent-baroclinic levels for a two-parameter 
model of the atmosphere. 

It is true that the two-parameter models 
described above permit unstable waves to 
develop, which have many characteristics in 
common with actually observed disturbances, 
but some essential features seem not to be 
accounted for. The assumption introduced 
both by Eapy and Eırassen regarding the 
distribution of the vertical velocity means 
that the level of non-divergence is a horizontal 
or quasihorizontal surface which does not vary 
with time. It is easily seen from equation (8) 
that in order to remove this restriction we 
must use a more detailed field of z. If £ denotes 
any of the variables x, 7. t we obtain from (9) 


+5 © 
% 4) * 


If only one function dA,/dp £ o as Eliassen 
assumes, nothing is gained by introducing 
more than one function B,,(p). If on the other 
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hand both A, (p) and A, (p) are mutually in- 
dependent functions which vary with p we 
may incorporate a more realistic field of ver- 
tical motion. Since at least one of these func- 
tions contains second powers of p the static 
stability may vary horizontally in this model, 
even though only two parameters are needed 
in order to describe the field of z. Now being 
able to introduce two parameters for describ- 
ing the field of vertical motion we need a 
fourth equation to determine the problem. 
This fourth equation may be obtained in many 
different ways but it seems most natural to 
the author to try to allow for a more refined 
lower boundary condition by assuming that 
By (po) # By (po) # 0 and by using equation 
(7) to obtain a relation between w, and w. 
If we extend the model in this way it does 
not seem possible to derive a similar corre- 
spondance between the finite difference method 
and the parameter representation similar to 
the one indicated above. This is intimately 
related to the fact that the level of non-diver- 
gence now is not necessairly an isobaric sur- 
face. However, further discussion will have to 
be deferred until more experience has been 
gained through actual computations. 

As yet very few attempts have been made 
to use more than two parameters for de- 
scribing the field of motion. ARNASON (1952, 
1953) has discussed such models and it is 
clear from the presentation above that such 
models can be derived in many different ways. 
However, the work involved increases very 
rapidly, mainly depending upon the fact that 
the non-linear terms of the type v : \77 will 
be split up in N? terms when the expressions 
(9) and (10) are introduced. It is therefore 
questionable whether this is the most practical 
way to proceed if we want to incorporate 
additional details. It may be better to use finite 
difference methods also along the vertical 
axis, an approximation that becomes the better 
the more layers are used. 

If using a parametric representation of the 
kind discussed here it is quite important to 
choose the parameters in such a way they are 
not correlated to each other since otherwise 
the full descriptive capacity of the multiple 
parameter representation is not utilized. This 
criticism applies to both ELIASSEN’s and EADY'S 
2-parameter models as to ARNASON S multiple- 
parameter model. In principle this is the reason 
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why a more detailed field of motion may be 
described by means of only two parameters than 
is the case in the models mentioned above. 


3. A simplified graphical method for 
forecasting with a two-parametric 
model 


The following analysis is based on the model 
suggested by ELIASSEN, but may also be 
applied to Eapy’s model. It seems not possible, 
however, to extend it to more detailed models. 
We shall neglect the vertical advection of 
vorticity as well as the term k : Vo x dv/op 
and finally assume that ¢ < f in the divergence 
term of the vorticity equation. The definitions 
of h and A (p) are arbitrary to a certain degree 
and will be stated at a later stage. The two 
forecast equations derived by ELIASSEN (1952) 
are 


à 4 
= ve 
Ir + 
FRÈRES VAG ef) EV EN Cr bv 
2 oh _ 
Ver—r ($F PV ony Hi} Em 0: 7 (13) 
where 


BANS. pt) = Z (x,y, t) + Hp) + 
HARM BOs une 


F= 7 Vzxk vr=2Yhxk ( (4 
c= 5 Vie u a: 


and r are constants determined from 
Am (p) and By (p) and h, is a constant = nor- 
mal value of h. These two equations were 
obtained from (4) by assuming 


Po Yo 
of V(x, y, pt) dp= 0; Jp: V(x, y, p,t)dp=0o 
(15) 


Except for the constants these equations are 
identical with those derived by Eapy using a 
method similar to the one outlined in the 
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previous section. The differences in the 
constants are insignificant from a principal 
point of view. 

In order to be able to use a graphical com- 
puting method similar to the one advanced by 
Fyortorr for the barotropic case we must 
transform (12) and (13) into two conservation 
theorems (cf CHARNEY, PHILLIPS, 1953). In fact 


it is possible to find two quantities each of 


which are invariant at one level. ELIASSEN 
pointed out that if two levels p, and p, are 
chosen so that A (pı) and A (pz) are the two 


roots of the equation 


[A(p)—bA(p)—a =0 


(16) 


we may write 


at 
where 
Jı = [v + Alp) vr] : VIE +f + Alpı) orl 
Jo = [¥ + A(p.) vr] : VIE + f + Alpe) Er] 


Multiplying (18) with A (p,) and A (p,) and 


adding to (17) gives us respectively 


di = dl 
HN) rf (F +v-Vh) <0 


where the individual derivatives d,/dt and d,/dt 
refer to the levels p, and p, respectively. Ac- 
cording to the definitions vr : \7h = o and 
thus (20) and (21) may be transformed into 


7 i +f)—Alp) rf Lo =) 
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Introduce the finite difference expression for 
&; according to 


Dore 


ce fed (zi (23) 


—z;) (i =1, 2) 


where m is the scale factor for the map pro- 
jection used, d is the gridsize used in computing 
2;, which represents the mean value of z in the 
four neighbouring points. We now define h to 
be the difference in height between the two 
levels derived above. Thus also A (p) is de- 
termined. Introducing a function J (g) as 
Fyortorr (1952) does we obtain 


Cae 
ag 


OA 3 4Je@] Koveya 


21 — K(p) (@: — 22) ] = -v° 


where 


(25) 


To bring K; (p) inside the V-operator with a 
minimum of error we introduce a quantity G 
being a constant and defined as the mean value 
of z, — 7, over the area. Thus we obtain 


Ki (9) vi V (21 — 22) = Ki (p}vi- 
-V(ti=2,—G)= Vi V (Ki (9) : 


G}) G) -7 Ki (9) 
(26) 


"12123 


(21 — 2; 


where the second term may be neglected. We 
finally obtain 


5 [Zi — 23 + J(p) — Ki(p) (21 —z.—G] = 
VE — 2 + Je) — Ki(@) - 
, (2, a G)| (i 1; 2) (27) 


These two equations are well suited for 
graphical solution. 

We shall now compare these equations 
with those obtained by using finite differences 
in order to evaluate the derivatives along the 
vertical (cf. CHARNEY, PHILLIPS, 1953). We then 
choose three levels at equal distances: level a: 
250 mb; level b: soo mb; level c: 750 mb. At 
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level a and c we get by applying the vorticity 
equation 


da Wb 

eG. +f) = ae (28) 
d. Wp 

di (©. ee Ap (29) 


(4p = 250 mb in this case). The finite differ- 
ence form of the adiabatic equation at level b 
gives 


0 
dt (&— a) Have "V(2— Za) = u 


(30) 


s being the vertical stability. If we introduce 
certain mean values for s, and o, which are in- 
dependant of time and put vj = 1/, (v, + v.) we 
can express &, in terms of quantities at the two 
levels a and c and introduce in (28) and (29). 
The two equations obtained in this way will 
be of exactly the same form as (22), except for 
the numerical constants. However, in deriving 
(22) to be valid at p = p, and p = p, we have 
considered the normal stratification of the 
atmosphere, while the three levels chosen in 
the finite difference procedure were chosen 
more or less arbitrarily. By defining A (p) 
from actual soundings we can derive the two 
levels where the vorticity equation preferably 
should be applied in finite difference form. 

The interpretation of these two levels may 
be expressed as follows: If the wind field is 
assumed to vary with height in the manner 
expressed by (14) there are two levels (but 
only two, depending upon the non-linear 
variation of the vorticity advection) with the 
aid of which the mean vorticity advection of 
the atmosphere and the differential vorticity 
advection within this atmosphere may be 
expressed. 


4. Discussion of numerical values and 
accuracy of representation 


a. Determination of the levels py and p, 


ELIASsEN (1952) has given approximate 
expressions for the empirical constants that 
enter the equations and that determine the 
two levels p, and ps. As far as the author 
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knows one month’s radiosonde observations 
were used and A (p) was determined with the 
aid of equation (14) using the method of 
least squares. It should first be pointed out 
that this function may be determined in 
many different ways and ultimately it should 
be chosen to give the best possible forecasts. 
However, the following remarks will be given 
already here. 

Equation (14) is used in various ways in the 
equations (12) and (13) namely in order to 
compute the density or temperature, the 
geostrophic wind and the geostrophic vor- 
ticity. In order to obtain these quantities one 
or several differentiations of z have to be 
performed. If now A(p) is determined by 
fitting the observed z-values according to (14) 
we cannot be certain that this will be the 
best function À (p) to represent the wind field 
or the vorticity field with the aid of two param- 
eters. We have determined A (p), first by 
fitting equation (14) directly to the observed 
values of z, secondly by computing a series 
of hodographs at different points on the map 
using in each case three radiosonde observa- 
tions, then determining the variation of the 
thermal wind with height and finally de- 
termining the function A (p) which will give 
the smallest error of 


Vz — \/Z = A(p)- Wh (31) 


In both cases z was first determined from 


which integral was evaluated by finite differ- 
ences using Ap = 100 mb. As h we have here 
chosen the observed difference of z between 
200 mb and 1,000 mb to represent the ther- 
mal structure of the whole troposphere. 
Admittedly a fairly small amount of data has 
been used for these determinations: about so 
radiosonde observations in a winter situation 
and about the same number in a summer 
situation. In both cases the material covers 
the whole North American continent. The 
function A (p) has been computed for win- 
ter and summer separately and is given in 
table 1 column 2 and 3 respectively and 
graphically in fig. 1. The method first indi- 


cated above has been used to obtain these 
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Fig. 1. À (p) computed for a winter and a summer situ- 
ation. A, (p) and A, (p) for the summer situation. 


values. Practically the same results are obtained 
if equation (31) is used. Some differences be- 
tween summer and winter are observed as is 
to be expected, however, mainly at high 
altitudes. The variations between tropical and 
polar regions probably are larger, but are 
difficult to take into consideration since the 
same value of A(p) must be used over the 
whole area in the equations (12) and (13). 
The values given here differ considerably 

from those given by ELIASSEN (19 52). Because 
of this the values of a and b also become differ- 
ent. We obtain a = 0.10, which corresponds to 
a = 0.30—0.35 with Eliassen’s definition of h. 
b seems not to be significantly different from 
zero. From equation (16) and (14) we then 
obtain 

| pr & 350 mb 

| pe & 825 mb 


which thus are the approximate pressures at 
the equivalent-baroclinic levels of this two- 
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b. Accuracy of a two-parameter representation 


There are two different factors that have to 
be considered when going from simple to 
more complicated models of the atmosphere 
for forecasting purposes. Primarily, how would 
it be possible to treat the hydrodynamic 
equations more exactly without increasing the 
computational work too much, secondly, how 
much details may we incorporate in our mo- 
dels without going beyond the accuracy of the 
input data, the temperature soundings? It is 
therefore of particular interest to investigate 
the accuracy of the representation given by 
equation (14) or more generally by assuming 
N = 2 in equation (9). A few estimates have 
therefore been made using the summer case. 
Column 4 of table 1 gives the standard devia- 
tion (o,) of the value of 2,4; (x, y, p,t) — H (p) 
and in column 5 the standard deviation (o.) 
is given of the errors of z (x, y, pst} ae 
using the representation proposed by Eliassen. 
The standard deviation of the corresponding 
temperature error (07,) at various levels is 
given in column 6. It is clearly seen that con- 
ditions near the ground are not very well 
described, while the temperature field in the 
rest of the troposphere is represented with an 
accuracy of about 2.5° C. Between 200 and 
300 mb the accuracy is again slightly less, 
obviously depending upon the variation of 
the height of the tropopause. It should be 
remembered that the function A (p) used here 
has been determined to fit this particular case, 
and thus these figures probably become 
slightly larger if the application is made to 


parameter model of the atmosphere. another weather situation. However, the 
Table 1. 
ee ee ee 
cm REND) Ma EE [vs BE ee 
A(p) A (p) O1 O2 ü 3 5 O4 
mb winter [summer | 10 | 10 = A,(p) A,(p) Io Ts Os Io 
case case feet | feet C feet C° | °C/roo feet 
| 
1000 — 0.53 | — 0.47 18 Io 0.II — 0.47 8 I 
900 — 0.4I | — 0.38 17 6 hee 0.27 — 0.36 4 05 0.50 2 
800 — 0.30 | — 0.27 23 4 3.0 0.49 — 0.25 4 0.30 3 
22 — 0.19 | — 0.16 31 be 0.70 — 0.16 5 0.25 4 
00 O0 050.03 2 i 0.96 —= ©.06 6 ; 0.10 5 
1.5 : | 
> > a 2 7 5, 1421 0.05 7 a 0.10 ÿ 
4 .24 27 : 9 ae 1.55 0.18 8 EG 0.15 9 
300 0.41 0.47 2 Io 43 1.94 0.33 8 : 0.15 12 
200 0.47 0.53 88 Io 3.5 1.99 0.56 8 3-5 O.15 15 
100 0:28 | - 0.07. | 39 21 50 = 
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values of A (p) seem to vary only a small 
amount from one case to another as can been 
seen from a comparison between the winter 
and summer cases (col. 2 and 3). Normally 
we therefore may expect these values to be 
perhaps 10—20 % larger, at 200 and 100 mb 
maybe still somewhat more. 

In section 2 it was very briefly indicated 
how a somewhat more adjustable 2-parameter 
model might be obtained by trying to de- 
termine two empirical functions A, (p) and 
A, (p) instead of just one as in Eliassen’s mo- 
del. Before developing the equations that 
would govern such a model it is of some in- 
terest to see if such an attempt would give a 
more accurate representation of reality. In 
the columns 7 and 8 of table 1 these two 
functions are given as obtained in this par- 
ticular case (cf fig. 1). As our parameter z, 


we have chosen 2 — H (p) and for the other 
parameter 2, 


Za (x, y, t) = [zo (x, y, 200, t) — H (200) — 
— A, (200) : 2, (x, y, t)] — [200s (x, y, 1000, t)— 


— H (1000) — A, (1000) : z, (x, y, f)] 


In the columns 9, 10 and 11 the standard de- 
viations of the errors in this representation of 
z (x, y, p,t),{o3], and the corresponding standard 

eviations of the errors of the temperature 
(ors) and static stability (o,) are given. A 
slight improvement is obtained if compared 
with the representation proposed by Eliassen, 
but it is surprisingly small. The error of the 
determination of stability is quite large in the 
layers close to the surface of the earth but varies 
between 0.10 and 0.15° C/100 m in the middle 
and upper part of the troposphere. However, 
the errors of this representation still show 
certain systematic characteristics that probably 
are real. This cannot be judged, however, 
until we have made a comparison with the 
accuracy of the various radiosondes in use at 
present. | | 

À comparison between six different kinds 
of radiosondes has recently been made in 
Payerne, Switzerland (cf. e.g. NYBERG, 1952). 
From this investigation we have determined 
the standard deviation of the error in com- 
puting the height of various pressure levels 
and these values have been entered in column 
12 of table 1. Above 600 mb the errors in the ob- 
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servations are as large as or larger than the errors 
introduced by a fwo-parameter representation. 

In relatity the situation is somwehat better 
than these figures indicate. In USA for example 
the same kind of radiosonde is used over the 
whole continent which reduces the errors, 
in particular when computing gradients of 
z. However, these figures clearly indicate that 
the exactness of the dynamic equations must not be 
carried too far but must also be determined from 
the accuracy of the observations. 

The fact that the errors in the representation 
of z remain comparatively small up to about 
200 mb depends upon our choice of the second 
parameter as the difference between 200 mb 
and 1,000 mb. In order to increase the accuracy 
at lower levels, where the observations justify 
it, this second parameter should be chosen as 
the difference between for example 300 mb 
and 1000 mb or may be even 400 mb and 
1000 mb. 


5. An example of a baroclinic forecast 


We shall finally give an example of fore- 
casting with the two-parameter model using 
the graphical method indicated in section 3. 
We have chosen a case over USA in which a 
inarked baroclinic development took place. The 
300 mb and 700 mb maps have been used since 
the heights of the 350 and 825 mb levels were 
not available. The variable = has in this case 
been computed from Z = 0.4 [z(x, y, 300, f) 
+ z (x, y, 700, t)], while h has been assumed 
to be h = (z (x, y, 300, t) — z (z, y, 700, f)). 
As a comparison a barotropic forecast with 
Fjortoft's method has been made using the 
soo mb map. The initial maps as well as the 
observed and forecasted 12-hour changes are 
shown in figs. 2—10. The effects of the 
mountains have been included in a way that 
was described by Born, CHARNEY, (1951). 
In tables 2 and 3 some figures are given to 
express the success of the forecasts. 


Table 2. Correlation coefficients between observed 
and computed 12-hour changes 


Obs. soo-mb change — Comp. barotropic 
ELVA SE nee nr ann a dass eeepc 0.72 
Obs. 300 + 700-mb change — Comp. baro- 
| CUDIGES TE SZ Chan ee Sharer oer cie RON 
) Obs. 300 — 700-mb change — Comp. baro- 
ChiniG 2 SN CHA Eee 0.72 


Fig. 2. Contour-field of the 700 mb surface Feb. 6, 1951, 


15 GMT. Lines are drawn for every 100 feet. 


BE RMS 18) VEIN 


Fig. 3. Contour-field of the 300 mb surface Feb. 6, 1951, 
15 GMT. Lines are drawn for every 200 feet. 


Fig. 4. Computed 12-hour change of the 700 mb surface 
from Feb. 6, 15 GMT—Feb. 7, 03 GMT 1951. Unit 
100 feet. 


Table 3. Standard deviations of observed and 
computed 12-hour changes (100 feet). 
NOLE so DEchan es Ee ren er ye 2 
\ Barotropic forecast (using 500 mb) ........... 1.9 
{ Obs. chance ot Boole 700 mb 4.6 
\UCompschansegore ee en a 4.4 


| Obs change of 300 — 700 mb 
\ Comp. change of z,— 2, 


A slight improvement is obtained if com- 


Fig. 5. Observed 12-hour change of the 700 mb surface 
from Feb. 6, 15 GMT—Feb. 7, 03 GMT 1951. Unit 
100 feet. 


pared with the barotropic forecast and in 
particular the intensification of the trough 
over the Middle West is better described by 
the 2-parameter model. When inspecting the 
maps more closely there are two facts that 
should be brought out. 

1. The computed centre of maximum fall is 
systematically displaced to the north of the 
observed centre, an observation that has also 
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Fig. 6. Computed 12-hour change of the 300 mb surface 
from Feb. 6, 15 GMT—Feb. 7, 03 GMT 1951. Unit 
100 feet. 


Fig. 8. Computed 12-hour change of the soo mb surface 
from Feb. 6, 15 GMT—Feb. 7, 03 GMT 1951. Unit 
100 feet. 


been made in Princeton (CHARNEY, PHILLIPS, 
1953). | 

2. The observed changes are considerably 
more concentrated into sharp maxima than 
is computed. This seems to be an indication 
that the influence between different points is 
too large in the way the forecast is made. It 
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Fig. 7. Observed 12-hour change of the 300 mb surface 
from Feb. 6, 15 GMI—Feb. 7, 03 GMT 1951. Unit 
100 feet. 


Fig. 9. Observed 12-hour change of the 500 mb surface 
from Feb. 6, 15 GMT-—-Feb. 7, 03 GMT 1951. Unit 
100 feet. 


may also to some extent depend upon using 
too large a grid-size, which will smooth the 
field. Some smoothing is also introduced by 
the graphical method. 

However, the final judgement cannot be 
made until a series of forecasts has been 
made, preferably by using a machine. 
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Abstract 


In a case of strong development the tendencies at the 500 and 1000 mb surfaces are 
computed with a two parametric model. It is possible to take the thermal structure of 
the atmosphere into account in a crude manner by expressing the vertical temperature 
distribution as a function of the deviation of the 1,000—500 mb thickness from a standard 


value. 


The height of the soo mb surface and the 1,000—s00 mb thickness are used as 


variables. The tendencies are computed for two maps 12 hrs apart, the intermediate tend- 
encies interpolated for every 2 hrs and the 12 hr change of the 1000 and soo mb surfaces 
are then computed. The tendencies computed represent an improvement over those 
computed with the barotropic model, but the 12 hr changes are still only about 2/; of the 


observed changes. 


The following tendency computations were 
started during a stay at “The Institute of 
Meteorology, University of Stockholm” in 
spring 1952 and finished in Oslo. 

The different atmospheric models have been 
treated by Dr. ELIASSEN 1952. The model I 
have used in the following computations is 
given there. I will therefore only repeat the 
main assumptions and the equations which 
form the basis for the computations. 

We try to represent the temperature field 
as a sum of a “normal” vertical temperature 
distribution and a deviation from this normal. 


T= Ty + per (p) (h — h„) or spec. vol. 


Œ— 0, +g Ay’ (p) (h— h,) where Ai’, (p) is 
a function of p only and:h —h, a function of 
x, y, t, chosen to give the deviation of the 
1,000—500 mb thickness from the “normal”. 
The ’ means derivation with respect to p. The 
assumption leads to the following vorticity 
equation. 
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+ AYh) VHS VI ZH AV . 


+ 2 dA'© V?h— 
2 | 
i à [€ + V:Z+AV°h o)|- O 
dp 2 / 
where A=A,—A’, Z,= = 
We can now in different ways from this 
equation form two equations which can be 


used to find Z, and h,. We have, however, to 
assume something about the vertical velocity w. 


We put = W (x y t)-B(p) and assume 
adiabatic motion. Then 

O, Ome CO; 

ray a NF O ut) ra) =O 


220 
or 
fp 7% Op _ ay 20 
= + V _ = Gv) ra) O, p dp 
; ee Op... 
We ignore the variation of @ with x, y, f, 


and the deviation of « from the chosen normal. 
We then get 
gAırh, | gAlVh y. pop 


Ky an oO 


hy + rk x VZ-Vh=h,+ 


IR ge poe _%n 
rie Ve 


Here the left hand side is a function of x, y, t, 
only, and we must have 
A 
An OF 


[u 
s 


KE ak 
i 


If we now apply the equation (1) to the 1,000 
mb level and also take the mean value of 
equation (1), we get the following two equa- 
tions. 


V? Z, +k x VZ-9(f +492) + 


+ Ak x Vie VE +2 A @ 7h = 0 
(2) 


VZ + Ao V'h+kx (WZ + 


+ Ae Th) Of 772 f 


£ AGE h) — 


— (£ EA, Vth) Op =O (3) 


[03 
le] 


In the following we shall use the equations (2) 
and (2)—(3). If we now introduce the ex- 
pression for W in these two last equations, 
we get the following two equations which 
form the basis for the computations. 


STGIUIRIDESESINDEIBEYZE 


WA h, u E (£ + v:Z) AL 


a) Veh | Leo 


& 2 ‘By (f 27% ( f 
EM (2 + v2) + By’ + 


A =) vith | ok NAS AE 


+k x View| S48 (OZ + Av | — 


A® 
ee: Vie VE Vth= 0: a 


Vy, VZ.v(r+#v2) = 
+ Ark x Ve DEV h + 2A B(h + 


+kXÉVZ-Vh = 0 (5) 


Here Z can be expressed by Z,, the height 
of the 500 mb surface. We now can solve the 
equations if we can find values of A and B. B 
is, as we have scen, given when A is given. 
To get a first estimate of the values to choose 
for A, I used the ascents for about a month 
from an arbitrary station, Berlin, and calculated 
T— I 
h—h, 
all ascents. I could then find a mean value 
of A,’. The values of A,’ thus found started 
from near zero at the ground and increased 
gradually to a maximum at about 300 mb. 
From there it dropped to zero between 225 
and 250 mb. Higher up the values were 
negative. From A,’ we may get A. With 
some small adjustments these values of A were 
used in the computations. Since 


V B= Vi Ze (Ai (s00) — A,) Wh 


- T, and h, were the mean values for 


we must have A, (500) — A, (1,000) = —1. 
The function A (p) must also give a reasonable 
variation of wind with height. The functions 
A (p) and B (p) are shown in the figures 4 and 
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Fig. 1. Observed change in soo mb contour height 
03-15 G.M.T. 4-12-1951. Lines are drawn for every 
80 m with some intermediate 40 m lines stippled in. 


s in Dr. Eliassens paper. The constants which 
enter into the equations are. 


A 08 A? = 0.6, JT ll AAN Nome. 
BE O27. a Om 


The unit for B is mb/decameter. A is a pure 
number. 

I have used the same area and gridsize in my 
computations as those that were used in the 
barotropic tendency-computations in Stock- 
holm. The distance between two points in the 
grid is about 300 km. From the map we must 
read Z;, the height of the soo mb surface 
and h, the thickness 1,000—5$00 mb. Then the 
following terms, or combinations of terms 
must be computed for every point in the grid 


kx Va V(f+Évez) kr V2. 
4 82 7. se Sh ya hh 
iv ni à 
kxVZ-Vh kxVZ VE Veh 


We solve first the equation for h,. It has the 
form V7? h, — ah, + b = o where a is positive 
almost everywhere, and there are no difficulties 
in finding h, by relaxation. The values we find 
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Fig. 2. Observed change in surface pressure 03-15 
G.M.T. Lines are drawn for every 10 mb with some 
intermediate 5 mb lines stippled in. 


for h, may now be used in equation (5) which 
is also solved by relaxation. 

For a numerical test of the model I picked 
a situation with fairly strong development. 
Whesday chosen was Deesuyoscaulnsche 
figures 1 and 2 we see the observed 12 hrs 
changes at soo mb and at the surface from 
03 to 15 G. M.T. The soo mb and 1,000 mb 
tendencies were computed at the beginning 
and end of the period. 

A major difficulty is to verify the results. 
The computed tendencies are instantaneous 
tendencies, and they can only be compared 
with a difference between two maps 12 hrs 
apart. The best procedure would of course be 
to make a series of short time steps. The amount 
of work involved in such a scheeme is, how- 
ever, very large. 

If we simply take the mean of the two 
tendencies, we can not expect to get a good 
result in a case like the one I have treated with 
rapid movement of strong tendency centres. 
We would get a picture with two maxima of 
falls and an area with smaller falls between 
them. To get a better verification I have tried 
to interpolate the intermediate tendencies for 
every 2 hrs. 05-07-09-11-13 G.M.T. The 
isallobaric centres were assumed to move with 
constant speed and to change intensity linearly. 
In the same way the position of the o lines were 


Fig. 3. Interpolated 12 hrs change of s00 mb contour 
height 03-15 G.M.T. Lines are drawn for every 40 m. 


found by letting each point on them move 
parallell to the centres with constant speed. 
The rest of the intermediate tendencies were 
then interpolated on the maps. The change 
from 03—15 is then found by letting the first 
tendency work 1 hr, the next 5 2 hrs and the 
last ı hr. The 12-hour changes obtained in 
this way are shown in the figures 3 and 4. In 
these two figures the lines are drawn for every 
40 m. 

These two maps may now be compared 
with the two observed maps in fig. I and 2. 
We see that the position of the o lines and 
the centres of rise an fall are in fairly good 
agreement with the observed positions. Both 
rises and falls are, however, too small. Except 
for the falls over Balicum on the 1,000 mb 
map, they seem to be about 2/3 of the observed 
values. The two soo mb tendencies computed 
with the barotropic model are both much 
smaller than the two computed with the two — 
parametric model. 12 hrs changes computed 
from them seems then be much too small. 

It is possible that the two-parametric model 
might give better agreement with another 
choise of A. In the equation for h, the terms 
containing 


k x vie v(r+8vz) 
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Fig. 4. Interpolated change of 1,000 mb contour height 
03-15 G.M.T. Lines are drawn for every 4o m. 


k x VA:VEVch kx<V2-Yh 


are usually the most important. The constants 
in these terms depend on A and B. It is the 
values near the ground which come in as a 
consequence of the way I have chosen to apply 
the equation. That is perhaps not a very good 
choise, as the values near the ground may be 
difficult to estimate. 

In the equation for Z,, it is, apart from the 
barotropic term, usually the terms with the 
constants A® and A’B which are most im- 
portant. B is proportional to A’ as we have 
seen. The windshear along the vertical is 
determined by A’. Since A is zero at the mean 
level, big values of A’ will also give big 
values of 42. Another possibility is to choose 
B without letting the adiabatic equation be 
fullfilled for every value of p. 

It may be of some interest to study the ver- 
tical motion which results from the adiabatic 
equation with the assumptions made earlier. 


We hveo = W (x y 2) -B(p) = (he = 
7k x V/ Z;:\/h)- B(p). Here all quantities 


are known and we can calculate w. The vertical 
motion at the beginning and the end of the 
period is shown in the figures 5 and 6. The 
surface fronts are shown as dashed lines. The 


Tellus V (1953), 2 


COMPUTATIONS WITH A TWO-PARAMETRIC MODEL 


Fig. 5. Computed vertical motion 03 G.M.T. Units 

are given in the text. A stands for ascending motion, 

D for descending motion below 250 mb. Surface fronts 
are stippled in. 


model gives the same direction of the verti- 
cal motion up to 250 mb and above this 
level it is in the opposite direction. The 
figures will therefore give the horizontal dis- 
tribution of the vertical velocity at all levels. 
The unit used on the map will, however, 
change from level to level and is given in the 
table below. 


800 mb I cm/sec 
700 1;7 
500 6 
300 2,8 
250 O 
225 — 3 
100 =, 
75 0 


Let us look at the conditions below 250 
mb. In the fig. 5 we see an area of ascending 


Tellus V (1953), 2 


77 
= 
Neh} 


Fig. 6. Computed vertical motion 15 G.M.T. 


air above the tip of the deepening cyclone. 
Ahead of it and behind it there are areas of 
descending motion. There is also an area of 
fairly strong descending motion behind an old 
depression off Norway. Above and ahead of 
the warm sector near Denmark we see an 
area of ascending air. A considerable pressure 
fall occurred over Balticum. In the fig. 6 we 
find the same general picture. In the tropo- 
sphere the ascending motion over the deepen- 
ing cyclone is considerably stronger, and the 
descending motion behind it is also quite 
strong. We have still the area of descending 
motion ahead of the deepening cyclone. 
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